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On the Theory of Diffraction 


By W. FRANZ 


University of Minster 
Communicated by D. Gabor; MS. received 28th November 1949 


ABSTRACT. A new method is described for calculating the diffraction of an acoustical 
or electromagnetic wave by successive approximations. Kirchhoff’s theory is a special 
case of the first step in the new method, which, unlike Kirchhoff’s, is not restricted to black 
screens, and applies to long waves as well as short ones. 


§1. BASIC IDEAS 

HE diffraction theory of Kirchhoff is, though very useful, not entirely 

satisfactory in its theoretical foundations. The main defect is that there 

exists no higher approximation, thus it is not possible to estimate or to 
amend the errors. Moreover it applies only for diffraction by a black screen, 
while diffraction by an object of arbitrary reflectivity is also of great interest. 
The deficiencies and difficulties of the Kirchhoff theory can be avoided by a new 
conception of the fundamental ideas, resulting in a consistent method of successive 
approximation to be described below. 

Let us take as zero order approximation a primary wave up, satisfying the wave 
equation for empty space, ignoring all diffracting objects. ‘The next step is to 
take account of these neglected objects, and this may be done in a purely experi- 
mental way by some plausible modifications of uy. We may, for example, cancel 
Uy in the geometrical shadow of an object, and add a reflected wave roughly 
according to geometrical optics. ‘The resulting modified function up’ now, of 
course, contradicts the wave equation in some parts of space, for instance, at the 
boundaries of the geometrical shadow. ‘Thus we have next to restore the wave 
equation by adding a first order correction du, in order to adapt the first order 
function u, =u, + du, to the wave equation, which determines du, unambiguously 
if only outgoing waves are taken into account. Now again dw, is calculated 
without regard to diffracting objects, therefore wu, has to be modified yielding uy’, 
which again violates the wave equation and therefore demands a second order 
correction 8u,, and soon. The whole procedure consists, as we see, of two kinds 
of steps: ‘intermediate steps’ u,—>wu,’ taking into account the influence of 
diffracting objects, and ‘principal steps’ u,’—u,,,, adapting uw to the wave 
equation. The first order of this method applied to a black screen yields 
Kirchhoft’s theory. 
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SP IDIRVONKEINNIL) SAM ales) Oo APPROXIMATION 


Consider an acoustical wave uw, satisfying 


1 0? 
(\- a5n)¥=0 e ereneuene (1) 


or an electromagnetic wave, given by the electrical field strength E and magnetic 
intensity H, which satisfy Maxwell’s equations 


culH=<— ; divi =O 1 
oH («= 3) eiotenae (1 a) 
curl E= —p—> ; divH=0; 


To perform the (v+1)th step of our approximation, we start from the modified 
function of the vth approximation u,’, E,’, H,’. The corrections dw, ,,, dE, .1, 
5H,.,, must satisfy the inhomogeneous equations 


io: = 
(a = 2 sa) Olea ae oe ee Em oined (2) 


curl 3H, )=¢ 5 8E, 1-5 edivok, 4 —=p,; 


0 : 
curl OE, athe, oho =j 3) Wel pee 


Here the right-hand sides 


~ po kG, 
=p = a ro 5) ew ee caer (3) 
: E,’ 
j,= —curlH,’ +« = ; py = —ediv E,’; 
oH’ See eh (3 a) 
j,* =curl EY +p ae p,* = —p div H,,’ 


play the part of fictitious sources of radiation; in the electromagnetic case they 
represent electric and magnetic charges and currents. The solution of equations 
(2) and (2a) 1s well known; if we mark the coordinates of the observer by a suffix 
P, the coordinates of the radiating point by a suffix J, we have 


1 / [e(3)] 
Ge) 
sul = Boe Sd se eee (4) 
! LW) [I] 7 y 
SE®), =— ig ytd Ua A ( )I 
»-_ Loe Ge eeeca ees to*(d 
®t Ve p 
OH 4m curl, | rps an 47 Ot TPS oe 7 ary 
eee (4a) 


The squared ’brackets indicate retardation, meaning that within the brackets the 


time t has to be replaced by t—rps/c. rp, is the distance between P and J, and 
dz, the differential of space. 
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§3. INTERMEDIATE STEPS: SURFACES OF DISCONTINUITY 


The intermediate steps are rather arbitrary, and it requires skill to choose 
suitable modifications of w, fit for the needs and conditions of a given case. 
We have to keep this freedom of the intermediate steps in order to obtain the best 
possible service from our method. ‘This means that our approximation method 
does not generally converge or diverge; if we perform the intermediate steps 
shrewdly enough we can ensure convergence (e.g. by just guessing the correct 
wave function); if on the other hand we choose u,/ unwisely, we can, of course, 
prevent convergence. 

The most important way of proceeding is to restrict the modifications of u, 
to surface discontinuities. ‘That means: we follow our free propagating wave 
u, till it meets the surface of the first obstacle, here we stop and consider from 
equation (4) how to continue the path of the radiation. We take, so to speak, a 
pair of scissors and clip our radiation field along the critical surfaces—these ‘cuts’ 
according to equation (3) represent a surface density of radiation sources, exactly 
as Huyghens’ principle demands. For illustration we see in Figure 1 a reflecting 
screen with an aperture. ‘The arrows indicate the direction of propagation of 
the wave, and the broken lines show where the scissors have to work: we clip 
along the opening, and also along the back surface of the screen immediately 
after reflection (the ‘reflection zones’ of Figure 1 are meant to be infinitely 


Figure 1. Remecas screen with opening. 


Reflection Zone 


narrow). Let us now calculate the secondary radiation du,,,. Of course the 
derivatives which we want for equation (3) are not defined for our discontinuous 
functions, but we may approximate these by continuous ones and proceed to 
the limit after integrating equation (4). In that way we get by Gauss’ generalized 
theorem in the electromagnetic case 


fi. ar = [do x HY; le. di =e [do hee 


ocean (5) 
[i*¢r=— |doxE,’; [p.* dr =p [do . Hy’. | 


The vector surface element do points, as shown in Figure 1, towards the clipped 
part of the space. For E,’, H,,’ we have to take E,, H, along the aperture, and the 
field strength of the reflected wave along the reflection zone. As (5) 1s valid 
within infinitely small regions of the surface, where rp; is constant, we may 
insert a factor 1/rp, and retard each single integrand of (5): 


(lB ap, = fda, x HD , [OD ape (day EOD, | 


Vpy Ypy J Lps | 
j ; é FJ : H,'(J 
(GAO a, = — | day x . Die | ce, ee 
py . rps Tp J TPs 
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These formulae enable us to apply (4a) to surface discontinuities. In the scalar 
case we have to be more careful because of the second derivative contained in 
(3), which makes it wrong just to put 1/r and retardation into integrals analogous 
to (5). We get the correct result by twice applying Gauss’ theorem to the exact 
function (4): ; 

OD fea = div p | es do; — = [esc dG 70 Sie (7) 


If we insert (6) into (4a), we obtain the corresponding formula of the 
electromagnetic theory correctly for closed surfaces. However, we meet with 
difficulties in the case of surfaces which are not closed surfaces, as here the surface 
currents demand line charges along the boundaries, which are not considered 
in (6). The simplest way to overcome these difficulties is to eliminate the p’s 
from (4a) and to express everything by the j’s, which can easily be done after 
differentiating with respect to time; since we are sure that (4a) satisfies the 
Maxwell equations (1a) identically in j and j*, we take from (6) only the first 
right-hand terms of (4a) and obtain the rest by application of (1a): 


0 valve: [E,’] 1 [Eo a 

5 5B a= G5 out! | do x ; + g-curleurl {de x rat . 
sitet (7a 

a a Pa Lea ee 

5, Ha = 7 5 curl | do x = moe Fa Curlew | dex oul 


(7) as well as (7 a) satisfies the field equations (1) and (1 a) identically, for arbitrary 
choice of u,’, grad u,’, E,’ and H,’ as we see at once: every integrand is a function of 
t—r/c, divided by r, and therefore a solution of (1), consequently curl curl in (7a) 
can be replaced by grad div—(1/c?)(@?/dt?). (Note that curl curl = grad div — A.) 

The equations (7) and (7a) are just the formulae of Kirchhoff’s theory. 
Equation (7a) was not hitherto known in this simple and useful form. Some 
early attempts to arrive at these formulae were refuted by Kottler (1923), as 
violating Maxwell’s equations. Kottler thought it necessary to introduce a 
correction by applying equation (7) to each component of E and H and adding 
certain line integrals. “The formulae (7a) were first derived in 1947 by the author * 
(Franz 1949) for monochromatic waves, but the above derivation for arbitrary 
time dependence is much simpler and clearer. Another generalization—arbitrarily 
variable anisotropic « and —is derived in another note (Franz 1948). 

If the radiation is monochromatic with time-dependence exp(—iwt), d/dt 


means a factor —iw and retardation a factor exp (¢krp;), where k=w/c, so we get 
instead of (7) and (7a) 


ieee exp ik 
Bea = 5 tm (I) grad, ( pt “2 ) caUCANE ree (3) day, 


py Tps 


SEeee pout [dex po ee eee 
An ‘ r Atiwe 


exp tkr | 


if 


curl curl | do x 


exp kr 


i 
8H, = curl | do x HL’ os 


1 

—— curl curl | do x EY’ 
Aniwp a 
Some obvious suffixes J and P are omitted. 


* See note at end of paper. 


’ 
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Though (8) is identical with the Kirchhoff integral and (8 a) is its consistent 
electromagnetic analogue, our altered conception of the fundamental ideas 
yields a series of essential improvements. First, of course, we have an infinite 
Sequence of approximations instead of a single not very clearly established 
approximate formula. But even in the first approximation we get important 
advantages from our new aspect: we can apply (8) and (8a) as well to the space 
before as behind the screen, to screens of arbitrary reflectivity or even transparency 


and to long waves as well as to short ones—which is all outside the scope of 
Kirchhoft’s theory. 


$4. SMALL PERFECTLY REFLECTING OBJECTS 
If the diffracting objects are small compared with the wavelength, the appli- 
cation of formulae (8), (8a) is straightforward. Let us assume perfect 
reflectivity by postulating the boundary conditions 


= = Ee = Omeat.the surtace, 99 0 ie sos (9) 

n is the direction perpendicular to the surface, and the suffix || indicates the part 
of E which is parallel to the surface. Now here the example of Figure 1 is not a 
good model, because the small object does not succeed in building up a shadow 
or a reflection front, but is completely illuminated by the wave, without distinction 
between front and back, so that these terms are now meaningless. For the 
diffracted wave starting from the surface, we first take over the whole surface 


Pree ASN 5 ea a? (10) 


in order to satisfy the boundary condition (9). In order to obtain zero identically 
when wu, or E, satisfies equation (9), we put in addition 


Us = Ve Tie actnecsitiace. a) = a eee (11) 
ikr 
which gives OU a= = (— SLE GM at CERO en, ee ae ese (12) 
1 eikr 
SE, 4.5 — gcurl [do x B,—. ee (12.4) 


We disregard H in the calculation of the individual steps of approximation— 
finally it can be derived by Maxwell’s second equation. 

For a point P, whose distance is large compared with the dimensions of the 
object, the first order radiation is easily determined. As uw and 7 are regular 
within the object, we can apply Gauss’ theorem and obtain 


Lyey ; exp tkr j 
6u,(P) = ae ie : div, (Seen grad, uf3)) ars 


PJ 
. . ik 
etl { G EADIE TE ay eorad (ae) _grad, uf3)] dr, 
An object Vpy rps 


Rk 
dE,(P)=—- Be curlp curl; — B,(J) dr 


4a object PS 
bs , : R > 
etl | (iopti(3) x grady “PE mime iinterad apn) 
4n object Tpy Yes 
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The integrands and their derivatives are practically constant inside the object» 
consequently the integration yields just a factor V equal to the volume of the 
object, and the result is: 


V exp tkrp exp tkrp 
du,(P) = — a {eee Uy t+ gradp = ee grad uo} ketene (14) 
tk exp zkr 
dE,(P)= S {iy Hy x gradp ea —curlp (E, x gradp Stee), : 
7 


ee: (14) 


By omitting the suffix J we take every argument not marked by Pat the place of 
the object, which is almost just a point. In this approximation the object 
enters only by its volume, and its shape has no effect. ‘This is different in the 
second and higher orders, where the shape of the object is of decisive importance. 
In the case of a spherical object the result of the yth approximation is 


LP) Pe ees {1 z 3(1 z >) 8 8} ee (15) 
E (P)—E,(P)= 7 ee kx {FE < ((1 - (3) )&- 5(1 a >) ke) 
ee: (15a) 


Here k is the vector wave-number of the diffracted wave in P, ky the vector 
wave-number of the primary wave at the object, © the angle they include (for 
the calculation see Franz 1949, §$4 and 8). Equations (15) and (15 a) show good 
convergence towards the well-known exact solution of the problem. The para- 
meters of approximation are pure numbers; this somewhat surprising fact 
might have been expected since our successive corrections have to amend the 
violation of the boundary condition, and must therefore be of the same order of 
magnitude as the starting function, even if numerically smaller. Another puzzling 
point is why the first step does not itself produce the exact result, seeing that the 
correction has just the correct discontinuities in 0u/dn and in E, to cancel the errors 
in the surface values of the primary wave. The solution of this paradox is, of 


course, that d5u,/dn and dE,, are not zero inside the object, so that there is a 
residual error outside it. 


$5: PERFECTLY REELECDTING HALE-PILANE 


An instructive example is the diffraction of a plane wave by a perfectly 
reflecting half-plane. ‘This emphasizes the most characteristic features of the 
new method. Each step of approximation can conveniently be compared with 
the well-known exact solution of Sommerfeld. 

For simplicity we assume cylindrical symmetry, choosing as primary radiation 
a plane wave incident at right angles to the edge of the half-plane. We describe 
the points of space by cylindrical coordinates z, p, d, with the axis in the edge, 
¢=0 in the open half of the diffracting plane (see Figure 2), and |¢| always less 
than or equal to 7. We have to consider only the scalar wave equation, since for 
cylindrical symmetry the electromagnetic field can be reduced to a scalar one 
(with v= E if E is parallel to the edge, and u =H if H is parallel to the edge). 
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Let us first assume as boundary condition 
is) eat, Doth) sides Of the screen, sees (16) 


The primary plane wave may proceed into the direction ¢p, and its wave function 
may be written 


Up, p)=exp{tkpcos(d—do)}; we (17) 


time dependence, as mentioned above, is as exp (—iw?). 


Figure 2. Diffracting half-plane. 


Screen Cut 


Figure 3. Zero approximation for half-plane. 


In the first intermediate step we add a reflected wave at the near side of the 
screen, and use our scissors to cut away everything beyond the plane of the screen 
(Figure 3). This gives 

cA) gradu, =0 for ¢>0 


Up (p, $) = Ud(p, $) —U)(p, — ¢) l Bee (18) 
grad uy’ (p, 6) = grad ue(p, 4) + grad uo(p,—4) | 


(The condition ¢) >0 is essential only for Figure 3 and formula (18), but can be 
dropped for the rest of the calculation.) 

To get the first order correction du, we have to apply (8), integrated along the 
surface of discontinuity ¢=0, and insert the discontinuity of wu’, which is zero 
for up’ and 2gradu)(p, 0) for gradu’. ‘This gives 

exp thr 


Re [ere = 
du, (p, ) = Ini K i dzy ib dps 7 Sin P41 eXp (thpyCOSho)" «swe (19) 
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Integrating dz, brings in the first Hankel function, which we represent by 
Sommerfeld’s integral 


+0 ikr 12) 
J dz; = =tnH ({R(p? + ps —2pp, cos $)"*5 


=e) 


=i | dbyexp {ik(o? + ps? 2pp, cos $)!* cos $3} 
Cc 


=) { dé, exp { —ikp,cos$, +ikp cos(|¢|—¢1)}.  -+--e- (20) 
C+(a+[o)/2 

The complex path of integration C is shown in Figure 4. The last form of (20) 

results from changing the variable of integration; the path is shifted into 

C+(7+|4|)/2, which makes the integral convergent for all values of p and py. 

If we now deform the path once more to cross the real axis at 6; =7, we can carry 

out the p; integration with the result 


in 
dd, exp {tkp cos (|¢|—4,)} Seb Pome ile es (21) 
| 


tf 
AGE = ai cos, — cos do’ 


crossing at 7. After integrating we can shift the path completely into C + Id |, if 
we respect the rule of crossing the real axis. Formula (21) is a very cony enter 
integral representation, which can easily be worked out for large as well as for small'p. 


C+\¢ 


Figure 4. Integration path. Figure 5. Half-plane, boundary condition u=0. 
Exact solution 
First approximation ——-——-—-—-—-—- 
Second approximation —.—.-—.—.-.-— 


Next we have to calculate the second approximation. The situation is that 
Uy =Uy + du, satisfies the field equations everywhere, 1,’ also satisfies the boundary 
condition (16), but dw, does not, which makes a second approximation necessary. 
We have to form a reflection zone according to Figure 1, and to take for u,’ there, 
of course, —8u,. We need not bother with gradu,’, as graddu, is zero for 
¢=+7. But the screen is not only reflecting, it also throws a shadow, and this 
contributes, as du, is penetrating the screen (it does not ‘know’ anything about 


‘ 
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the screen), another —8u, at the surface of the cut at p=n7. So we get first 
for d>0: 


(ke e kr et kr 


eyes it 
y= 


Sox (22) 


This formula accounts for the boundary condition at the screen, since the dis- 
continuity is now found to be —25,, which gives by reasons of symmetry just 
half the amount, i.e. —8u,, at 7+0, and +8u, at 7—0. (If we calculate for the 
region ¢ = — 7, we have the opposite sign of n, and so there also the correct surface 
value.) We carry out the dz;-integration as in (20): 

r+ © ikr F : 
[bss =i (Rp py?-+ 2ppy cos} 


= | dy .exp {ikp; cos by + ikp cos(|$|—$y)}. se eee . (23) 
C+|¢]/2 


The path C has to be shifted by |4| /2 to ensure convergence for every p and py. 
Again we have a pure exponential of p,; and can integrate if we cross the real axis 
at ¢.=0. The direction of do, is towards diminishing 4, and therefore div 
equals — (sin $(0/dp) + (cos ¢/p)(0/ed)) of the scalar integral. So we get ford >0: 


Co er te fi Re) ee 
dus(p, 6) = Dene? cos , + cos $y Jap does tkp cos (| 4| DD cre aeal coage 


Pee 


d¢, crossing at zero. Here again we insert the modulus of ¢ in order also to include 
negative values of ¢. We can move the path for d$, completely into C+|4| 
without touching a singularity (¢,= +4,+7), as $,+7 at C is certainly out of 
reach. ‘The dd, integration is elementary: 


sin Jp SiN do 
oer 


+ cos ,)(cos f, + Cos dg) 


—%0 


sin bs e'Ps + iho sin Jo Ss ett + elds ] 
=> tea eS ee ee n IP Gey Se Ie! Et RS Ee ee 3S bas mae 
cosdg—cosdy e+e cosd,—cosdy e+e |... 


= do sin $y — $2 Sin bo ; 


cos dy — COS do 


For the evaluation of the last logarithm it is essential that 4, crosses the axis at 
zero. We get finally 

1 : do Sin dy — by Sin hg 2% 

Sualer$) = za; |, ba exp {ike cos 6] — dy) ore EY se (26) 
Now no longer does it matter where we cross the axis as the integrand has no: 
singularity except at infinity; the numerator vanishes at ¢,= +49 as well as the 

denominator. 
To illustrate the result let us calculate uw, and u, for large p and compare with 
the exact solution. In this case only the very neighbourhood of the saddles 
1=|¢| and 4, =|¢]| contributes to the integrals in (21) and (26) (in (21) we pass 
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a pole of the integrand when |¢9|>|4| but this contributes only the geometrical 
wave in which we are not interested). This gives 


exp {i(kp + 77/4)} sin do 
du,( 0, ) = — en Ricoh, Seiecer (27) 
; xp filke +7/4)} 1 dosin| 4l—[d]sin dy 
ee = a cet erry (28) 


The exact Sommerfeld solution is 
ei(ket+7/4) 2 cos 4¢ sin $9 
geo (2rkp)'? cosP—cosdo — 


Figure 5 shows the angular factors of u,, Us, us, Figure 6 the errors of the first and 
second approximation. Near the screen the improvement gained by the second 
approximation is very appreciable. Near the opening the improvement is small; 
on the contrary we pay for the improvement at the screen by a discontinuity of the 
normal derivative in the opening, which is the reason that we need a third 
approximation. 


' 
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Figure 6. Errors of the first (- --—--— ) and Figure 7. First approximation for half-plane 
second (— -—-—-~—-~—) approximation. at the edge, boundary condition u=0- - - - 


a u 
boundary condition — =0 —— 


on 


Now if we proceed to the third approximation we meet with a characteristic 
difficulty: at the edge the function du, is finite, but not continuous—its value 
depends on the angle. So its tangential gradient is infinite there, and it is just 
this gradient which we need for 5u, according to (8), and its singularity makes the 
integral diverge logarithmically for every pp, dp. As this difficulty already appears 
in the second approximation, if we take the boundary condition éu/dn =0, we will 
now consider this case. 


§6. THE DILEMMA OF THE EDGES 
For the boundary condition 
Ou 


> =3().. at. the'screen, — 9) =o. | eee (30) 


we modify the incident wave (17) again according to Figure 3, only we take for 


the reflected wave the opposite sign of (18) in order.to satisfy the altered boundary 
condition. That gives 


a= OE grad uy =0 for 6 >0 
Uo (p, “) = Uo(p, ) ar Uo(p; ii f) 


fo 0, are eames: 
grad uq (p,$) =grad up(p, $) — grad u(p, — 4) rgo< (31) 


‘ 
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Here gradu,’ has no discontinuity, but the discontinuity of uy’ equals 2up. 
Therefore we get according to (8) 

1 


ikr 


5 el) | ee 
uy(p, $) ee da. grad, ( A ) wale 0) 
be alae | pe exDAZKT ar tkpy COS Qo) 
= r div | do, ; ey ees ts eae (32) 
The dz, integration is carried out as in (20), div this time equals 
sin $(0/0p) + (cos f/p) (0/0) of the scalar integral, so we get 
jira er ; sing 
Suled)=5- 7 | a k — $,)} —— "1 __ 
iPr?) == I¢] Jos+isi Prexp Uke cos(lol—d1)} cos $,— Cos bo’ 
sits oe (33) 


crossing at 7. 
Now let us consider the behaviour of (33) as well as (21) at the edge. In (21) 
for p=0 we can put | for the exponential and the integral becomes 


aay d sin do ees | sin do 
277i} cin |) Ccos¢,—cosdy 2nt Jy '*cosd,+cos dy 

1 erhr + eibe ] —10 

Shee ng | = bo ie eve (34) 

In (33) we cannot simply put p=0 under the integral as the integrand would not 
then vanish at infinity. We have first to transform the integral slightly. We 
shift the crossing point for z to |¢|, which makes an extra contribution of 
(¢/|4|) exp {tkp cos (|¢|— | 4o|)} if |6|<|¢ |. Then we transform the variable 4, 
by subtracting |], after which the path of integration is C. As this path is 
symmetrical, we can replace ¢, by —4, in the integrand and take the average 
between this integral and the original one. Now we can put p=( and get 


pees. in oat fem ecete | a {OF for |¢ol>|¢| 
4ni | | cos (¢,+|4|)—cos do 0 for |¢o|<|4| 


? ? (35) 


aT ead ee ee ran 


~ [él 


Figure 7 shows the functions (34) and (35) of the angle ¢. Inspecting it we must 
admit that our functions for the boundary condition du/dn =0 cannot but behave 
in this way, because it is just what we demanded: an unambiguous function with 
a jump of 2u,) at 6=0 and regular everywhere else. To compensate the discon- 
tinuous increase and to join regularly at 6= +7, du, must decrease continually, 
so this troublesome singularity of thé tangential derivative of our function is 
essential and necessary. For the other boundary condition «=0 such a thing 
does not happen as there is no discontinuity of uy’ at 6 =0, only of the gradient, 
which gives no difficulties in the next approximation. 

Now the infinity of the gradient at p=0 not only frustrates the next 
approximation, but is a very serious fact in itself, as it makes the energy of the field 
infinite—there is in the expression of the energy density a term proportional to 
(gradu)2, and as this behaves along the edge like 1/p’, the integral diverges 
logarithmically. This must be, of course, excluded,* therefore we must subtract 


+400 


* The necessity of such a condition was first pointed out by Meixner (1948), who also proved that 
this edge condition is sufficient to determine the diffraction by an edge unambiguously. 
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something from (33) which removes the singularity without adding a singularity 
anywhere else. In our example of the half-plane, this may be easily done in 
different ways; we have only to replace the fraction of the integrand (33) by another, 
which has the same residua for ¢, = + ¢p, and vanishes in an integrable manner at 
the boundaries of C. Such a substitution is for instance 


sing, sin po 1 (36) 


> Py ol eager ees 
cos d, — Cos dy dy Ccosd,—COSs do 


Another possibility is the substitution 


Smid 2 SINGS = SCs ee (37) 
cos J, — COS do sinddy cosd,—cos dy 


This yields the exact solution of Sommerfeld. 

For a general kind of screen a regularized solution cannot be found in such a 
simple way. But in any case we must insert a third kind of step into our approxti- 
mation method, if we are dealing with diffraction at edges, as in this case we have 
three conditions to fulfil—wave equation, boundary condition, edge condition. 
If we ignore the edge condition, there is another method of carrying out the second 
approximation, by a limiting process, which will be shown in the following section. 
This process can be transferred to arbitrary edges. 


§7. IRREGULAR SECOND APPROXIMATION 

Though 6u, in equation (33) displays at p=0 the singular behaviour (35), 
there is, surprisingly, a way of making the integrals of the next approximation 
converge without really altering 5u,. We note that 64% (Rp) 1s a solution of (1), 
which is regular everywhere except at the edge (and at infinity of course) and may be 
suitable to compensate the singularity of 5w, without violating the conditions of 
the problem. But for very small p the Hankel function tends to infinity, thus 
any finite multiple of 6% )(p) is far too large to compensate the singularity (35). 
But we can carry out the second approximation by the following limiting process: 
we choose a length « which is extremely small compared with the p of the observer 
and the wavelength, and replace du, by 


24 H kp) 


du, + nt H 6 NRe)* ataiouelkele 


The addition is twice the amount required to cancel the singularity of (35) for 
p=e, it satisfies the wave equation everywhere and produces no additional 
discontinuity at 6=0; thus our new function is suitable in every respect, and 
alters du, at the place of the observer only by a negligible amount. Inserting (38) 
into (8) we get the second approximation for the boundary condition (30): 

ikr oO 


il 0 + e 2 Wl (1) D 
oun d)= 35), aes | dex ae (Bulon dd 282 Sere) 


gy =x 


At p=e we cut off the surface of integration; in consequence Su, added to (38) 
will satisfy the boundary condition (30) everywhere, if p is greater than the chosene. 


, 
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If we calculate the normal gradient (= — (1/p,)(0/0¢,)) for (33) at 6, =7, we obtain 
an integral which we oe in the following way: 

; sind, 798 
J G44, loos ie COS bo Ps exp {tkpj cos () xy Pr) | 


ar f sing, | = 
J. ,de AH ee Beagcoe 1) 


aa sin $y 
pe | =o dd, exp (—ikp; cos ¢, Dar ree ery 


Equation (39) now = using (23) for 6>0: 


dp F He : 
dus(p, $) = I. aa = | svgyn U2ESP {ik0s 608 $+ thp cos ($—dy)} 


ely a 7] sind, 2 H \™(Rps) 
‘ ta | Coes BEER D roPps 08 G3) 0d | cosh, — COs by re H \(Re) } d 


Aer (41) 


We may neglect some terms which can be made arbitrarily small if we choose « 
small enough compared with the wavelength. We have first * 


1 
lim exp {7kp,; (cos dy — cos : =—C=-—lIny, 
erg [2 py “e | Pp Uy ps ( fe 1)} i] —tkp; (cos dy — cos 1) y, 


where C = the oe, constant 0-577 and y=exp es Note also that 


k 
Jim Fp ee “P. exp (ihps 008 bs) | (eps) — + aa oer Oaeatgrs fis (43) 


as the integral ee for «0 a finite value Rae and Oberhettinger 1943, 
chap.3, §1). This transforms (41) into 


Li 
Suis» 6) = 57 |, yn WPeerP (tke cos (b—$s)} 
(cross at 0) 
sa | oO sing: le dps See -Iny) 
\27tJo+n ''0¢,cosd,—cos bo \Jo ps(1 —tkp3(cos bcos $,)) 
47” dps \n(ykpy/22 ; 
=| = ee a exp (tkp; cos $,)} Si aaget ce = od ene Pree e (44) 


Integrating the last term by parts and neglecting again terms small for small ke, 
we get 


41 (dos, (ykps 
PH Re) | ger na (ae 9; ) &XP (tkpy Cos hs) 


47 1 ykpy : is 
= 3 f.O(he\ Whe) \2 15 | ( ey exp (tkp, cos 5) 


— tk cos $y i (In n oe) exp (tkp, COS fo) dps} 


bao tin ue Ta a hae a eee (45) 
ie 


* e.g. Magnus and Oberhettinger (1943, end of Chapter Ds 
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The remaining dp; integration in (44) is elementary: 


I remem ("+ Reged) 
« py{1—tkp;(cos $y — cos $;)} ke(cos dg — COS ;) 
—> — in +Iny—In2(cos¢,—cos¢;).  — - - +» (46) 


In the last line again we have omitted terms small for small ke and split the result 
in a useful way. Inserting (45) and (46) into (44) yields: 


Buster) =— 55 pa | aay (sa) 


cos ¢, + cos Jo 


< gh, [_ dbyln2(c0s 4, +008 4) exp fikpcos (| $]—4a)} 
TS E419) (47) 
| 
ral 
f 4 j! 
f Gd 
/ hy 
: a f Re 
T =A 
TT A ASG f 
“7 y re 
he if 
fe 


Figure 8. Half-plane, boundary condition be =0. Exact solution 


First approximation — —— -— ——- — . Second approximation —.—-—.-.---. 


Figure 9. Errors of the first (~-----— 
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By inserting the modulus of ¢ at the proper places we make the formulae valid 
also for negative ¢. The ¢,-integration can be carried out with the result 
Lee ahs dy Sin dy — dg sin d 
5 ona Papas 2oS81 Po 
ua(p, d) yy | | ; dd 


Brat 
ae > cos by — C08 $y » exp {tkp cos(|¢|—¢s)}. 


This is now independent of ¢, thus our limiting process works. As only 
the immediate neighbourhood of the edge is concerned with this process, it can 
be applied also to arbitrary screens, at least if the curvature of the edge is continuous. 

To illustrate the result (48) we calculate the value of the wave function for a 
distant observer, from (33) and (48): 


{u(kp + 7/4)} sind 

: __ exp filkp+7/4)} 

1y( oO, p) (2irkp)¥2 cos db a anS dbo DO, | oso oro (49) 
exp {u(kp+7/4)} 6 1 dosindy—¢sind 
dU, (00, 6) = — ce I ae 
ee Qnkpy® Ts] = —cosd—cos¢, 
ae (50) 
The exact Sommerfeld solution is 
Ces {i(kp+7/4)} 2sin4¢ cosddy 

Us( 0, p) = TAOS AUS Lane ne ere (51) 


cosh—cosdy — 


Figures 8 and 9 show the values and defects of the first and second approximation. 
By the second approximation we gain a reasonable improvement near the screen 
(mainly with respect to the boundary condition (30), which we tried to improve 
by this step), but we have to pay for it by a discontinuity of the function in the 
opening, which makes the third approximation necessary. 


Note added in proof (see p. 928). ‘This statement is not quite correct. As 
I have learned in the meantime, L. de Broglie in Problemes de propagations guidées 
des ondes électromagneétiques (Paris, 1941) arrived at an equation equivalent to (8a) 
as a preliminary form of Kottler’s formulae. Obviously he did not realize the 
importance of this formula, and it remained generally unknown until the 
present author stressed its usefulness (Franz 1948). 
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ABSTRACT. ‘The intensity distributions for three oxygen band-systems, the Schumann— 
Runge and Herzberg systems of O, and the second negative system of O2*, are calculated, 
and compared where possible with experimental results. Some suggested identifications of 
bands observed in emission from the night sky are checked by comparison with the 
calculations. 


§1. INTRODUCTION 


N an earlier paper (Pillow 1949) an attempt was made to estimate, by comparison 
] of calculated results with experimental observations, the accuracy with which 

the distribution of intensity among the bands of a system for a diatomic 
molecule could be predicted by an approximate theoretical method. The 
calculations were there carried out. on the y-system of MgH*, and were taken 
only as far as the v=8 vibrational levels. Within these limits, observation and 
calculation were found to be in close agreement. 

It has been suggested that in view of the probable existence of oxygen band- 
systems, both in emission in the night-sky spectrum and in absorption in the sun’s 
atmosphere, a theoretical estimate of the intensity distributions to be expected 
in these spectra might be of some use as a guide to the identification of suspected 
bands. Accordingly, the intensity distribution schemes have been worked out 
for the Schumann—Runge and Herzberg systems of the O, molecule, and for 
the second negative system of O,*. For the potential energy curves the Morse 
function has been used, and the wave functions for the vibrational levels are 
obtained as before by graphical distortion from those for a harmonic oscillator. 
For this purpose, the wave functions corresponding to a quadratic potential 
function have been calculated for vibrational levels up to v=16. ‘Table 1 is 
therefore complementary to the table of wave functions previously published. 

The intensity in emission of the (v’, v”) band is taken in the first place to be 


proportional to 
/ uw 2 
yt ( | bie bn dr) 


and the intensity in absorption to 


»( [Wewear), 


v being in each case the wave number of the origin of the band considered. This 
neglects the effect of vibration on the electric moment of the molecule, regarding 
this as constant. The general agreement of the results with those of experiment 
suggests that this is permissible. Also, the formula as it stands assumes that the 
vibrational levels in the initial state are all equally populated—that is, it assumes 
infinitely high temperature and the absence of any other factors causing preferential 


: 
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distribution. Where it seems likely that thermal motion is the main factor affecting 
the distribution this has been introduced in the Boltzmann factor 


exp[—hew,/(v! +4)/RO] 
for emission, or exp[—hew,"(v" + 4)/k6] for absorption. 
Itis realized that at the higher levels the effect on the calculated intensity values 
of the departure from normalization caused by distorting the wave functions may be 


appreciable, and the effect of this has been considered when dealing with the 
Schumann—Runge system. 


Table 1. Wave Functions for a Simple Harmonic Oscillator 
le. (2%! 4/77)" exp { — 457} H,(s) 


3 9 10 11 12 ils) 14 tS 16 

3 

0 0-393 0) —0-373 0 O23 517 0 — 0-344 0 0-333 
+ 0-128 0:339 —0-154 —0-332 Oeid2 0-341 —0-:077 —0-368 0-042 
4 —0-184 0-315 0-246 —0:248 —0-285 0-184 0-294 —0-122 —0-307 
3 —0:395 —0-044 0-357 0-157 —0:297 —0:241 0-218 0-290 —0-135 
il —0Q-234 —0-358 0-061 0-368 0-091 —0-321 —0-208 0-231 0-283 
1} 0-144 —0-307 —0-315 0-129 0-361 0-054 —0-331 —0-202 0-230 
14 02397; 0-051 —0-342 —0-267 0-164 Of3517 0-042 —0-318 —0-209 
13 0-310 0-353 —0:018 —0-356 —0-239 0-203 0-368 0-078 0-331 
2, —0-028 0-344 0-334 —0:043 —0-355 —0-254 0-155 0-356 0-102 
24 —0-349 0-026 0-359 0:320 —0:093 —0:394 —0-241 0-177 0-375 
24 —0:428 —0-318 0-051 0-357 0-316 —0-034 —0-327 —0-267 0-082 
23 —0:235 —0-427 —0-303 0-053 0-349 0-330 0-002 —0-312 —0-301 
3 0:094 —0-249 —0-424 —0-304 0-027 0-328 0-341 0-061 —0-266 
3h 0-366 0-081 —0-246 —0-418 —0-319 —0-005 0-301 0-352 0-113 
34 0-526 OESs 7/5) 0-087 —0:226 —0:399 —0-347 —0-104 0-198 0-346 
32 0-520 0-522 0-382 0-114 —0-192 —0-397 —0-373 —0-132 0-187 
4 0-418 0-513 0-522 0-400 0-154 —0:144 —0-364 —0-394 —0-204 
44 0-400 0-493 0-512 0-416 0-204 —0-:077 —0-315 —0-398 
43 0-174 0-277 0-392 0-487 0521 0-456 Mess: (0) == (517 
43 0-160 0-254 0-369 0-457 (ors i 0-466 0-320 0-087 
5) 0-045 0-087 0-151 0-238 0-342 0-447 0-506 0-497 0-389 
=1 

52 

54 0-008 0-018 0-036 0-063 0-110 0-178 0-207 0-353 0-434 
53 

s 0-001 0-003 0-006 0-012 0-024 0-046 0-064 0-131 0-201 
6L 

64 

6} 

i 0 0) 0 0 0 0-001 0-001 0-003 0-007 


Of the systems considered, the Schumann—Runge has been observed in the 
laboratory both in emission and in absorption, the Herzberg system, so far, in 
absorption only, and the O,* system, naturally, in emission only. ‘The experimental 
results in each case have been compared with those obtained from the calculations 
made on the corrected wave functions, and the agreement is, on the whole, good. 
This seems to indicate that the method of calculation is, with certain limitations 
to be discussed later, still sound even at fairly high energy levels, and that the 
estimated distribution of intensity under appropriately chosen conditions can be 
taken as a fairly reliable guide to the identification of bands suspected of belonging 
to one of these systems. _ 
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§2. THE SCHUMANN-RUNGE SYSTEM FOR Os 


This system of bands has been photographed in emission by several workers. 
The results were analysed, and further measurements made, by Lochte-Holtgreven 
and Dieke (1929). The system has been further studied in emission by Feast 
(1949), who has analysed a number of new bands. (Millon and Herman (1944) 
and Lal (1948a,b) suggest further bands, but in neither of these cases 1s an 
analysis published, These are discussed by Feast.) 

The most reliable observations on the same system in absorption in cold gas 
are those of Knauss and Ballard (1935). Fiichtbauer and Holm (1925) photo- 
graphed and measured bands appearing in absorption in heated gas up to temper- 
atures of about 1,050°c. 

The constants for the two electronic states involved are taken from the 
measurements of Lochte-Holtgreven and Dieke and of Knauss and Ballard. 
These give, with the usual notation: (i) for the excited state, w, = 700-36, 
w,/x’ =8-002, w,’y’ = 0-375 (neglected in this calculation), B,’=0-801; (1) for 
the ground state, w,” =1,584-9, w,"”" =11-645,-B,"” =1-46. 

The absorption bands in the v” =0 progression converge to a limit at about 
56,877 cm-, and the (0, 0) band of the system is calculated to have its origin at 
49,363cm-!. For the excited state, therefore, the value of D’ in the Morse 
function is 7,862 cm. 

The value of D” for the ground state can be calculated from the limit of the 
v’ =( progression of the Herzberg system, which has the same ground-state as 
the Schumann—Runge system, and whose two potential curves approach the same 
limit. It may also be found by considering the two atomic states into which the 
molecule can dissociate. This is discussed fully by Gaydon (1947). ‘The value 
adopted here for D” is 42,070 cm"!. 

‘The intensities in emission for the bands of the Schumann—Runge system, 
calculated as described above for infinite temperature and for 3,000°K., are set 
out in Tables 2 and 3. The bands recorded by Lochte-Holtgreven and Dieke 
and by Feast are shown in the tables. In these cases the emission occurred in 
some form of gas discharge, so that any assumption of thermal equilibrium is 
invalid, and for comparison with the experimental results it is probably useful to 
consider mainly the table corresponding to infinite temperature. 

As has been mentioned, it seemed advisable to try to estimate, at least approxi- _ 
mately, how far the distorted wave functions used here depart from normalization. 
This was done by making a numerical computation of the quantity fy,?dr for 
each of a number of wave functions, and comparing the values with those for the 
undistorted functions. The results for the wave functions of the ground state 
showed a gradual, almost linear, increase of the ratio from 1 to 1-5 as the 
vibrational level rose from 0 to 16. This means that the intensities in any v’ 
progression should be divided by the factor, between 1 and 1-5, appropriate to 
that value of v’. ‘This has been done for the first few progressions in emission, 
and the results are shown in Table 4. 

For the excited state, the correction is much more difficult to estimate, since 
the wave function curves are very greatly distorted for high r-values, and the 
quantity Jy,"dr tends to very large values when calculated from the graphs. 
But this great distortion occurs only in those parts of the wave function curves 
that are not in fact used, since they do not overlap the curves for the ground 
state. It seems wise, therefore, to attempt no correction for the intensity values 


; 
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in the v” progressions, but to bear in mind that the intensities calculated for high 
values of v are probably slightly too great. This is of importance only when 
considering absorption at very high temperatures. 


It is seen that of the bands analysed the great majority are those to which 
calculation attributes a high intensity. There are a few exceptions, the most 


Table 2. Schumann—Runge System 
Calculated intensities in emission at infinite temperature 
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Table 3. Schumann—Runge System 
Calculated intensities in emission at 3,000°K 
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Table 4. Schumann—Runge System 
Calculated intensities in emission at infinite temperature, corrected for 
failure of normalization 
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L=bands recorded by Lochte-Holtgreven and Dieke. 
F=bands recorded by Feast. 
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obvious being the (1,12) band claimed by Feast and the (2,15) of Lochte- 
Holtgreven and Dieke. Feast considers the (1, 12) band to be definitely present, 
though it is very weak, and exists in a region of many blended lines so that he has 
been able to measure fewer lines of this band than of others. He does not find 
the (2,15) band at all with the exposures he has used, and it is noteworthy that 
Lochte-Holtgreven and Dieke have recorded only lines corresponding to high 
rotational levels in this band. This suggests a possible explanation of the 
appearance of some bands which lie near the main Franck-Condon parabola in 
the vibration table, but whose calculated intensity is low. The method of 
calculation used here takes no account of rotation, and the intensities given are 
therefore those corresponding to the lines near the band-origin. — If rotational 
energy is taken into account, the potential curve is slightly displaced, and this 
could produce appreciable changes in band intensities. It seems safe, then, to 
say that if a band has a high calculated intensity some part of it should certainly 
be observed, but that a band of low calculated intensity, particularly if near the 
main parabola, may yet have some lines visible. If this is the true explanation, 
no such discrepancy should exist in a band system where measurements are made 
mainly on or near band heads, as in the Herzberg and O,* systems. 

Flory’s suggestion, discussed by Feast, that among other bands the (3, 15) 
and (3, 16) bands should appear, is not supported by the calculated intensities. 
The (1, 15) band, on the other hand, should according to calculation be readily 
observed, but apparently it has not been recorded by any workers. 

In considering absorption, the effect of temperature is very much greater than 
in emission for two reasons; the temperatures involved are, as a rule, low 
compared with those in emission, and the constant w,” is generally greater than 


t 


w,’. On the other hand, provided that the temperature is known with reasonable 
precision, the effect can be predicted fairly certainly, since there is little likelihood 
of the distribution of molecules among levels in the ground state being affected 
by anything other than thermal disturbance under laboratory conditions. (In the 
sun’s atmosphere this assumption is not justified.) The intensities in absorption 
for the Schumann—Runge system have been calculated for 300°xk., showing, as 
expected, only the v" =0 progression, for 1,300°K., this being about the highest 
temperature used by Fiichtbauer and Holm, for 2,300°x., which can reasonably 
be expected in a laboratory furnace, and for infinite temperature, which should give 
some guidance in examining the solar spectrum. These intensities are shown 
in Table 5-8. 

The slightly irregular variation of intensity in the v” =0 progression is of no 
significance, since the fundamental quantity 


[| verbo” ar)” . 


is very small for all the terms in this progression, and the errors involved in reading 
from the graphs are enough to account for the irregularity. 


§3. THE HERZBERG SYSTEM FOR O, 

These bands have so far been observed only in absorption in the laboratory. 
Herzberg (1932) observed seven absorption bands in oxygen at atmospheric 
pressure, and for six of these he made rotational analyses. He considered that 
besides the seventh band there was yet another, too faint for any kind of measure- 
ment. Chalonge and Vassy (1934) obtained the same series of bands in absorption 


Table 5. Schumann—Runge System Table 6. Schumann—Runge System 


Calculated intensities in absorption Calculated intensities in absorption 
at 300° k. atl. 35 0K: 
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5 ies 5 1 3H 6 6H Bi if 
6 5k 6 2 8H 14 128 8 3B 
7 20K 7 8 168 24 19H 8H Sj 
8 29% 8 12 Dip BS DoE 114 she 
9 51K 9 Mil 34 43 24 10 2 
10 50K 10 20 = 43 45 24 8 1 
11 78% 11 31 70 62 31 8 1 
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13 90K 13 36 = 100 a9 23 5 
14 Des 14 38 a7 55 19 3 
K=bands observed in absorption in H=bands observed in absorption in heated gas 
cold gas by Knauss and Ballard. by Fiichtbauer and Holm. 


Table 7. Schumann—Runge System 
Calculated intensities in absorption at 2,300°K. 
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Table 8. Schumann—Runge System 
Calculated intensities in absorption at infinite temperature 
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in the atmosphere, and their microphotometer trace shows clearly the eighth 
band which Herzberg presumed to exist, and gives a definite indication of one 
more. Their measurements did not extend to rotational structure. Herman 
(1939) made a large number of observations on absorption in oxygen at various 
pressures, and found that as the pressure increased the Herzberg system gave 
place to a triplet system of longer wavelength, the range covered by one system 
being apparently contiguous to that covered by the other. Herman observed 
and measured fully all the bands recorded by Chalonge and Vassy, and considered 
that there might be two further bands, of still longer wavelength, the rotational 
constants for which he obtained by extrapolation. Morguleff and Vassy (1939) 
also made measurements on the system in absorption in the atmosphere. Their 
record is very similar to that of Chalonge and Vassy. ‘They have made rotational 
measurements, which, however, give more erratic values for the rotational 
constants B’ of the upper electronic state than do those of Herman, or even the 
original values of Herzberg. Stopes-Roe (1948) photographed the system in 


Numbers in brackets are extrapolated values. 

? indicates a band considered to exist, but not measured. 
+ lines present, but too faint for useful measurement. 

* measurements doubtful near band origin. 


oxygen in the laboratory, and he has kindly allowed the writer the use of his plates 
to make yet again a rotational analysis. ‘The comparison spectrum used is that 
of a copper arc. All the values for band-head wavelengths and for B’ are shown 
in Table 9. : 

It can be seen that there is room for considerable doubt as to the wavelength 
of the (0, 0) band, and the constant B,’ to be used in drawing the potential curve 
for the upper state of the system. Herman’s B’ values are very much the most 
regular in their variation with v’. The values deduced by the writer from 
Stopes-Roe’s plates agree better with these than with those of either Herzberg 
or Morguleff and Vassy, and Herman’s values have been adopted here as far as 
the 2685 a. band—that is, as far as his actual measurements go. For the two 
doubtful bands of longer wavelength, Herman gives extrapolated values of B’. 
In view of the almost exactly linear variation with v’ of his measured values it is 
difficult to see why he has adopted these two values, which are well off the linear 


Table 9 
Values now obtaia 
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relation, and changes in the extrapolated values can so alter the calculated values 
of B,’ that, if the two long-wave bands are accepted, the equilibrium distance apart 
of the nuclei, 7,’, can be anything from 1-30 to 1-40. 

The ground-state of the system, since the absorption takes place in cold gas, 
is the same as for the Schumann—Runge system. Since the bands consist of a 
Q-series of lines only, with alternate lines missing, the value of B,” for the ground- 
state has to be used in finding B’ each time. 

Herzberg considered that there should be an absorption continuum at the short- 
wave end of the system, though his plates did not show it. Chalonge and Vassy 
(1934) and Gotz (1934) state that it certainly is there, beginning at about 41,300 cm-! 
(A=2420a.), though Gotz concedes that it may not be wholly due to oxygen. 
Stopes-Roe’s plates leave the existence of this continuum uncertain, as his source 
was weak in this region. 

Calculations have here been made on the intensities of the absorption bands 
in cold gas, using various assumptions as to the origin of the system, and the 


Table 10. Herzberg System 


Calculated intensities in absorption in cold gas, for various systems of 
numbering of bands 


v according Chalonge Herman 
to: oe and Vassy c “~ 
Ao (A.) To.=1°58 7e=1-44 re=1-30 71638 ro= 1-40 
v E. Uv I v I v I v iE 
2795:2 0 64 0 12 0 5 
2737°6 1 100 1 39 1 17 
2685-0 0 1 2 83 2 69 2 35 
2631-7 0) 1 1 5 3 42 3 48 3 OY, 
2594-0 1 3 2 12 4 17 4 92 4 71 
2554-6 2 i 33 26 5 18 5 91 5 82 
Doo eAl 3 15 + 58 6 9 6 79 6 79 
2489-1 + 28 5 74 7 3 7 TT 7 73 
2463-4 5 50 6 95 8 _ 8 100 8 80 
2443 6 83 7 99 9 — 9 41 9 63 
2429 a 100 8 100 10 _- 10 40 10 100 
2425 9 89 


results compared with visual estimates of the intensities in the photographs 
available (plates of Stopes-Roe and print of Morguleff and Vassy), in order to 
give some guidance in deciding between the various suggested positions for the 
origin. The photographs show a rapid increase of intensity from the long-wave 
end of the system to the 2489 a. band, and then nearly equal intensities, perhaps 
falling off slightly at the 2429. band. ‘Table 10 shows how very sensitive are 
the calculated intensities to change of r,’ and also to change of numbering of the 
bands. It appears from this table that the numbering of Chalonge and Vassy, 
which includes all the bands actually measured by any observers, taken with 
Herman’s rotational figures, gives calculated intensities whose order is in very 
close agreement with that observed, while the assumption of the existence of two 
further bands makes the long-wave bands more intense in theory than they 
actually appear to be, and it is difficult to see why, if the two bands exist, they do 
not appear clearly. Hence the numbering has been adopted here which makes 
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the (0,0) band of the system that with Aj=2685a. and B= 1:01, so that 
r, =1-44«10-8cm. The formula for the band-origins is taken as 
v = 36,791 + 787-2(v' + 4) —31-9(v' + 3), 


and with the band-limit at about 41,280 cm this gives D.’ = 4,433 for use in the 
Morse potential function. With these assumptions, the intensities in emission 
at infinitely high temperature have been calculated, and are shown in Table 11. 


Table 11. Herzberg System 
Calculated intensities in emission at infinite temperature 
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(For purposes of comparison with measurements, the wavelengths of the band heads 
are given.) 


Now a number of bands appearing in the night-sky spectrum have been 
tentatively attributed to the Herzberg system. In particular, Barbier (1947) 
made a number of exposures in the region 3100—4600a., and estimated the 
variation of intensity with wavelength, showing this graphically. Many of his 
peaks correspond to wavelengths appropriate to Herzberg bands, and in most of 
these cases he considers also that the profile of the band is of the right form for 
this system. Table 12 shows those of Barbier’s bands which can be identified 
thus, and beside each is given the intensity as calculated here, if this is sufficient 
to allow the band to be visible. It can be seen that there is good correspondence. 
Moreover, the bands which, according to calculated intensities, might be 
expected to appear, can, at least within the range 3100-4100a., all be identified 
with one of Barbier’s bands. (‘The calculations have not been carried to longer 
wavelengths.) ‘This is shown in the Figure, where the heights of the vertical 
lines represent the intensities of the corresponding bands. Above each band, or 
close pair of bands, is written one of Barbier’s wavelengths. It can be seen that 


there is an observed band to correspond with each of the calculated bands of 
appreciable intensity. 


, 
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In identifying night-sky bands, recourse has been had in some cases (Dufay 
and Dejardin 1946, Déjardin 1948) to the two further progressions, which are here 
rejected, and even to others beyond them. ‘The identifications are, however, 
very doubtful, for the photographs are all taken on very small dispersion, and in 


Table 12. 


Bands in Barbier’s Plates of Night-sky Spectrum which have Heads of 
Wavelength appropriate to Herzberg Bands 


Suggested Log I 
Ao (v’, v’) Semel (Barbier) 
3084 (2, 4) 23) 0:5 
3144* (1, 4) 40 1-0 
3164 (O52) 6 0-6 
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3228 @55) 13 0-9 
$295* (es) 22 ie? 
ae (0, 5) Tes 1-4 
3386 (2, 6) 9 het 
3453 (1, 6)? ‘lex 
3483 (Gaz) d-3 


Bands marked * are those considered by Barbier to have the correct profile for Herzberg 


bands. 


Values in columns 4 and 8 are read from graph published by Barbier, and give an 


indication of relative intensities. 
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the v’ =0, v’ =1 progressions as used here there are band origins of wavelengths 
sufficiently close to those in the two suggested progressions to be distinguished 
It seems, then, satisfactory to adopt this system 
of numbering, at least until further evidence is available, and if this is done the 
intensities as calculated here give some indication of the bands to be expected 


from them only very uncertainly. 


in emission. 
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§4. THE SECOND NEGATIVE SYSTEM OF O,* 

The second negative system of O,* consists of a large number of double- 
headed bands, the uniform interval between the heads being 200cm™. These 
have been observed by Johnson (1924), Ellsworth and Hopfield 1927), Stevens 
(1931), Mulliken and Stevens (1933), Lal (1948a,b) and recently by Feast. 
Mulliken and Stevens measured the lines of some of the bands and obtained the 
molecular constants for the two electronic states a?I],—x?H,, and discussed 
the dissociation energies, as Gaydon has also done (1947). In plotting the 
potential curves and wave functions for these calculations, the following values 
have been adopted for the various constants : 


(1) upper state: w,’=898-9, x,’w,/=13-7, D’=14,622cm™4, a=2-78 x 10°, 
reall On® creat 

(2) lower state: w,” =1,876-4, x,"w," = 16°53, D” =53,218 cm, a=2-81 x 10°, 
Pool TASS xe LO cine? 


Working with these in the same way as before, the intensities to be expected in 
emission have been calculated on the assumption of equal population in all the 


Table 13. Second Negative System of O,* 
Calculated and observed intensities 
Me 0 I 2 3 4 5 6 7 8 9 10; 5. Bip ae ena 
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9 9 11 1 1 3 1 1 1 
a¥ aF 4F 
10 Sal e7) 9 5 7 1 2 1 1 1 
aF aF 4F |F 
‘leds 6 Numbers in heavy type show calculated intensities. 
aF aF J=intensities observed by Johnson. 
12 | 14 4 F=intensities observed by Feast. 
aF 


a=bands observed by Feast, but no estimate of intensity made. 


[All the bands reported by Lal are included among those measured by Feast. 
Lal does not quote intensities. ] 


vibrational levels of the excited state. For the wave numbers of the band origins 
the component of longer wavelength has in each case been used. ‘The small 
difference between the components will make no appreciable difference to the 
calculated intensity of a band. 

Now Johnson estimated the intensities of the bands he observed, on a scale 
which made J=9 for his strongest band. Feast has made visual estimates of the 
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intensities of the same bands and of anumber of others, expressing them as nearly 
as possible on the same scale. The calculated intensities have therefore in this 
case been reduced to a scale which gives the observed values for a few of the 
strongest bands, so that comparison can readily be made. ‘T'able 13 shows the 
calculated intensities, those observed by Johnson and those observed by Feast. 

It can be seen that there is very good agreement between theory and obser- 
vation over the progressions involving low vibrational levels—that is, over most of 
the main Franck—Condon parabola. Even here, however, to a slight extent, and 
more noticeably over the subsidiary parabolae, the observed bands correspond 
to slightly higher vibrational levels in the excited state than do the calculated 
maxima, indicating a denser population in the higher vibrational levels than even 
infinite temperature would give under a thermal distribution. Feast has pointed 
out, however, that this is a likely result of the production of excited O,+* ions 
from molecules in the ground state of Oy. 

A number of bands in the lower left-hand corner of the vibration scheme, 
with high calculated intensities, have been observed by Mulliken and Stevens, 
and by Ellsworth and Hopfield. These are in the vacuum region of wavelengths, 
and the fact that they are reported as being faint does not constitute a serious 
disagreement with calculation. 

Lal’s suggestion that some of the observed bands would fit better a narrow 
Franck—Condon parabola receives no support from the calculated intensity scheme. 


$5. CONCLUSIONS 


This method of calculating intensities appears to involve many approximations, 
so that it might be thought that the results obtained by it would be unreliable. 
It has now been tested, however, by comparison with observation on several 
band systems and over a wide range of vibrational levels, and it seems to be 
justified by the very good agreement found. The results may therefore fairly be 
used as a guide in looking for new bands, or in identifying bands suspected of 
belonging to one of the systems discussed, and this should prove of some use in 


the investigation of spectra such as that of the night sky, the aurora, and the sun’s 
atmosphere. 
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Communicated by R. W. B. Pearse; MS. received 24th March 1950 


ABSTRACT. Bands of K, have been photographed in absorption in the region A 3924— 
X 3420 a. on a Hilger E1 quartz spectrograph. The bands lie in two groups, one between 


A 3924 and A 3686 a. and the other between A 3640 and A 3420 a., and appear to belong to 
two different systems. 


The bands in the first group have been analysed. They can be represented by the 
equation : 


v= 26493 -0+ 60-6u’ —0-15u’2—92:64u” +0°354u”?, 
where w=v+#. 


The upper state of this system is considered to dissociate into 47S-+5?P atoms of 
potassium. 


No analysis for the bands lying between A 3640 and A 3420 a. has been possible. 


§1. INTRODUCTION 
OTASSIUM has several systems of bands in the infra-red, the visible and the 
Pitesiote regions. Although the study of the infra-red and the visible 
bands is fairly complete as far as their vibrational structure is concerned, 
the position with regard to the ultra-violet bands is still far from satisfactory 
(Sinha 1947, 1948). 

The ultra-violet bands of K, have been reported in previous work (Chakraborty 
1936, Yoshinaga 1937, and Sinha 1945 a,b) to extend to wavelengths as low as 
A2900a. ‘The analysis of these bands provided by the various authors does 
not appear satisfactory and more work was therefore considered necessary. 
Consequently the K, bands between 44510 and A3940a. were photographed 
in the first order of a 21-ft. concave grating with a dispersion of 1:28 a/mm. 
These bands were found to belong to two systems, one lying between A4510 
and A4220a. (the blue system) and the other between 4160 and A3940a, 


(Sinha 1948). In the present investigation the lower wavelength limit has 
been extended to about A3680a. 
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§2. EXPERIMENTAL 

The bands have been studied in absorption. The details of the absorption 
tube used and the method employed for purifying the metal and heating it are 
essentially the same as previously described (Sinha 1948). An ordinary 6-volt 
motor headlamp was found to give a satisfactory continuum up to about A3400 a. 
which was fairly strong when the current through the lamp was slightly higher 
than the rated value. ‘To prevent rapid distillation of the metal to the cooler 
parts of the tube, nitrogen from a gas cylinder was introduced inside the 
absorption tube and maintained at suitable pressures. A Hilger El quartz 
spectrograph was used to photograph the bands. 

With a column of the hot potassium vapour about 20 cm. long, the bands in 
the region A 3900 to 43700 a. appeared most satisfactorily at temperatures between 
700 and 800° c., when the pressure of nitrogen inside the absorption chamber was 
about 15 cm.Hg. 


§3. APPEARANCE OF THE BANDS AND THEIR MEASUREMENTS 


The bands appear in two distinct regions : A3924—A 3686 a.and A 3640—A 3420 a. 
In general they are weakly developed although in the former region they appear 
somewhat stronger than in the latter. ‘The bands are all degraded to the red. 

Measurements of the band heads in the two regions are givenin Tables 1 and 2. 
Intensities given are visual estimates on a scale of 1-10. ‘The asterisks in Table 2 
indicate that these bands were in the region covered by the principal series line 
at \ 3447 a. which had broadened to cover a region of about 60.4. between A 3430 
and \ 3490 a. under the conditions of the experiment. 


Table 1. K, Bands: A3924—A3686a. 


Aaix (A-) Int. viv’ See Nair (A-) Int. Uv mae 
3923-6 1 — — 3784-4 7 D2 2:6 
3906-5 3 0, 10 1-2 3776-0 7 32 il) 
3894-0 3 OFS —2-+1 6775 8 2a 1-7 
3880-8 5 0,8 —1:1 3767-0 6 = = 
3867-6 5 On —O-1 3762°8 7 35 33 
3854-9 5 0,6 —2:5 3759-1 5 — = 
3841-3 3 0,5 0-9 3754-7 6 4,1 1-2 
3832-6 © 5 15 —0+3 3748°8 + — = 
3828-1 5 0,4 1-5 3746-6 7 Sell —0-3 
3819-4 5 1,4 —0-7 3738°1 7 Gnal 1:6 
3813-9 3 -—— — 3733-8 7 DH) —0:7 
3806:2 7 15.3 1:8 372558 6 6, 0 1:6 
3801-7 5 — _— Ses 5 — = 
3797-6 10 aes, 1:3 SND 6 750 1:5 
3793-7 10 il, 2 —2:2 3708-9 4+ 8, 0 3°5 
3789°2 7 by. os —0-1 3701°1 ? — — 
3686°1 ? — — 


§4. VIBRATIONAL ANALYSIS AND MOLECULAR CONSTANTS 

The bands between A 3924 and A 3686. have been found to belong to a single 
system. The assignment of their vibrational quantum numbers is given in 
Table 1. It has not, however, been possible to propose any satisfactory analysis 
for the bands in the region A 3640-A 3420. 
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Table 3 gives the (v’-v") scheme for the bands in the region A3924-A 3686 a. 
In conformity with the name proposed for the K, bands between A4160 and 
A3940.a. as the first ultra-violet system of K, bands (Sinha 1948), the bands 
between 3924 and A3686a. will be named the second ultra-violet system of 
K, bands. 


Table 2.. K, Bands: A3640-A 3430 a. 


Aair (A.) Int. Nair (A:) Int. air (A.) Int. 
3640°3 ? 2 3568-4 4 3503°3 = 
3632:2 ? 1 3560°9 3 3501:-1 e 
3612-7 ? 2 3557-0 1 3491-0 s 
3603°5 ” 3552°5 1 3483-6 e 
3600-0 2 3545-4 2 3474-6 = 
3594°8 3 303755 1 3461°8 A3 
3590-7 5 3529-9 ? 1 3451-6 e 
3586-4 5 3520-0 1 3441-5 = 
3583°5 = Soule 7a 1 3434-3 % 
3575-8 4 3505-0 ? 1 3429-9 a 
* See text. 


‘The bands measured in the present investigation between A 3924 and A 3686. 
can be represented by the equation 


v =26493-0 + 60-6u’ — 0-15u’2— 92-64u"+0-354u"2 (1) 


whereu=v+4. ‘Thetermsinw” are those from the work of Loomis and Nusbaum 
(1932). ‘The values of v.,0— rea. given in Table 1 indicate that the agreement 
between the observed values of v and those calculated from equation (1) is fairly 
satisfactory. 

The constants of the six states of K, now known have been summarized in 
Table 4. The constants v,, w,, WX, WeVey We% are those in equation (2), and 
are given incm'!: 

V=V,+@U— WX UO Val — WS Ut. a eee (2) 


Table 4. Vibrational Constants for the Ground and the Excited States of K, 


Dissociation 
States Ve We WeXe WeVe Wese D products 
1p + ~ 92-64 0-354 — — 4135 4°7S+4?S 
1p,+ 11683°6 69:09 0-153 — — S108 = 4457-4 2k 
aI. 15378:0 75-00 0-388 0:00437 0-00018 1780 4°S-+42P 


Dew 22970:0 60:60 0:20 = = 4560 »4 25-25 2P 
_— 2642777 ~=661-60 }§©3=s-« 0-90 0-0010 0-00030 OS OA aie 
oy 26493-:0 60:6 Oils = — 2309.8 Bao a oe be 


§5. DISSOCIATION PRODUCTS 


The energy of the dissociated atom is calculated from Vatom=V%o9+ D’—D", 
where the terms used have the usual significance. In order, therefore, to calculate 
Vatom one should be able to obtain a value of D’ from equation (1), the value of IBY 
being already known from previous work to be equal to 4135 cm™ (Loomis and 
Nusbaum 1932). D’ may be determined by applying Birge and Sponer’s 
extrapolation rule. The value thus obtained is approximately 6000 cm". 
Such a determination may be altogether out of order, first because only a few 
levels have been observed and consequently the range of extrapolation is too 
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large, and secondly because the values of mean differences given for the excited 
state in the last column of Table 3 are so irregular that it is not possible to judge 
from them whether the separation between the excited state vibrational levels 
actually decreases linearly, even over the short range observed, a condition which 
is essential for the application of Birge and Sponer’s extrapolation rule. 

With such a doubtful evaluation of D’ an attempt to determine the states of 
the atoms into which the molecule in the excited state dissociates may not seem 
justifiable, especially in view of the fact that the atomic states in this region are 
rather crowded. However, one may try to obtain the probable dissociation 
products from other considerations. 

It has been shown earlier (Sinha 1948) that, of the four systems of bands, the 
infra-red and the red systems are due to transitions from the ground state to the 
excited states disso¢iating into 42S and 42P atoms. The upper state of the blue 
system dissociates into 42S and 5?P atoms and that of the first ultra-violet system 
(A4165-A 3940.) into 47S and 3?D(?) atoms. Theoretical considerations 
indicate that there should be one more system of bands whose upper state would 
dissociate into 42S and 52P atoms, another whose upper state would dissociate 
into 42S and 32D atoms and yet another whose upper state would dissociate 
into 42S and 52S atoms, in addition to several other systems whose upper states 
would involve potassium atoms in still higher excited states. We have now to 
decide which of the above sets of dissociation products would be probable for 
the system of bands investigated in the present work. Now from the line spectra 
data for potassium we have 47S—5?5=21,025 cm“, 425— 32D =21,540 cm 
and 425—52P=24,710 cm. In the equation v,.,=%,9+ D’ —D”, assuming 
the value of v9 to be 26,480 cm™! approximately (see equation (1)) and 
D” =4135 cm! (Loomis and Nusbaum 1932), we can calculate the values of D’ 
necessary in order that the dissociation products may be one of the sets of atoms 
for which the line spectra data are as quoted above. ‘These values are found 
to be — 1320, — 805, or 2365 cm“! for dissociation products 478 +525, 42S +4 32D 
and 475+5?P atoms respectively. ‘The first two sets are thus automatically 
ruled out as D’ cannot be negative, and hence the only probable dissociation 
products are the atoms of potassium in the 47S and5?P state. The corresponding 
value of D’, viz. 2365 cm“, is low compared with the value obtained from Birge 
and Sponer’s extrapolation rule. However, this is not very unusual, as the extra- 
polation rule often gives a higher value (Gaydon 1946). A comparison with the 
values of the heat of dissociation for the 1X%,,+ and 1II, states of the molecule 
resulting from the 47S and 4?P atoms of potassium, given in Table 4, also indicates 
that the value 2365 cm~! for D’ is not unlikely. 
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A New Ultra-Violet Band-System of GeBr 


By E. B. ANDREWS anp R. F. BARROW 
Physical Chemistry Laboratory, Oxford University 


MS. received 12th April 1950 


ABSTRACT. ‘The vibrational analysis of a new system of GeBr has been made. The 


constants derived in cm are: 


Ve We XeWe 

: : 27252 197+, ee 

gi aa 27156 190-, ae 
2115/5 1150 

aT 295+, 0:75 


The upper state, probably ?A or “II, appears to converge rather rapidly to a limit at about 
S-oreveabove © —0) of 7175: 


(1937). They photographed a system, probably 72 —?II (ground state), 

lyingin the region 2950-3260 a., and gavea vibrational analysis. By analogy 
with GeCl, for which two ultra-violet systems are known (Jevons, Bashford and 
Briscoe 1937, Barrow and Lagerqvist 1949) it was expected that a second system 
of GeBr should exist at somewhat longer wavelengths than the 2X —II system. 
Spectrograms of this system have now been obtained, and the vibrational analysis 
is given below. 

The spectrum was excited in a stream of GeBr, vapour by an oscillator 
dissipating about 400 watts at a wavelength of about 30 metres. ‘This gave a very 
pure spectrum, consisting of GeBr bands, weak emission at shorter wavelengths 
from GeO, and Ge and Br lines. It was photographed on a Hilger quartz Littrow 
spectrograph (E.478) on Ilford Process plates. Exposure times were about 
30 minutes. The band-heads were measured against Fe arc lines whose wave- 
lengths were taken from the M.I.T. tables (1939). 

A reproduction of a spectrogram is given in the Plate. ‘The system splits 
obviously into two sub-systems, of which the longer wavelength component 
appears to be somewhat the more intense. ‘The vibrational analysis is straight- 
forward, and the results are given in ‘Table 1. 

Little can be said of the structure of the bands. ‘The 0,0(1) band certainly 
possesses two heads, at 27107-9 and 27102:1cm™, but the 0,0(11) head appears 
single under the present resolving power. Away from the origins, the vibrational 
isotope effect becomes marked, and on the long-wavelength side, four heads are 
often measurable. Comparison with the calculated shifts showed that the heads 
observed here were due to “Ge®!Br, “Ge7®Br, a blend of ”Ge”Br and Ge*!Br, 
and 7GeBr, in order of increasing wavelength. The measurements given in 
Table 1 are considered to relate to “Ge7Br (u =63-38 x 10-*4gm.), which is close to 
the average GeBr molecule for which 4. = 63-16 x 10™. 

PROC. PHYS. SOC. LXIII, 9—A 65 


Akt band-spectrum of GeBr was first studied by Jevons, Bashford and Briscoe 
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The derived constants are given in Table 2. Those given for the ground state 
have been obtained by treating mean AG” values from Jevons, Bashford and 
Briscoe and from the present work. Figures for the 2X state are included. 
The *I13/. 1). separation, from the analysis of the 2S — IT system, is 1150cm~ : 
using this value, the multiplet separation in the new upper level is 96cm~ as 
compared with 62 cm“! for what appears to be the corresponding state in GeCl. 
The figures are compatible with these states being either 2A or 211. 


Table 2. Summary of Constants for GeBr 


State Ve We XeWe 
“OE 33413 383°, 0:7 
22A 2 25252 Oven ayn 
ees 25156 190°, 795 
"TI 5/0 1150 ‘ - 
211 5 0 Don 0:7, 


The vibrational levels of the upper state appear to converge rather rapidly : a 
linear Birge—Sponer extrapolation gives D,(i) =1250, D, (ii) = 1360 cm-1, but it is 
probable that both estimates are a little high, for the fairly well identified position 
of the 3, 1 (11) band suggests that y.w/is not negligible. If we assume D,(i) =1100, 
D,.(11) = 1200, we get dissociation limits at 28,110 and 28,300 cm“! above v” =0 of 
II. The difference between these values is too small for it to be regarded as 
certain, but we can be fairly confident that there is at least one dissociation limit at 
about 28,200 cm! or 3-5ev. ‘The energies and products of dissociation in this 
group of molecules will be considered in a later paper. 
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A Large @-Spectrometer with Two-Directional Focusing for 
Precise Measurements of Nuclear Radiation 


By A. HEDGRAN, K. SIEGBAHN anp N. SVARTHOLM 
Nobel Institute for Physics, Stockholm, Sweden 


Communicated by F. D. Cockcroft; MS. received 14th December 1949 


ABSTRACT. In-the present paper a large double focusing spectrometer is described, 
which is convenient for precise measurements of B- and y-radiation. ‘The magnetic field 
form is discussed and also the arrangement to measure the value of the field continuously 
with an accuracy appropriate to investigations of this kind. The difficulties due to electron 
straggling in the B-source or converter are emphasized, and this effect is for the moment the 
factor which prevents the full utilization of the high resolving power of the apparatus. 
A number of typical spectra are given to illustrate the applicability of the instrument to 
different cases. 


IN ERO DG Calon 


Tt has long been recognized that one great disadvantage of the ordinary 

semicircular focusing principle, formerly extensively used in f-ray spectro- 

scopy, was the inability of focusing electrons in more than one plane, 
i.e. that perpendicular to the lines of force. For this reason the semicircular 
method has been almost abandoned in favour of the magnetic lens method, 
which has the important property of space focusing. 

For very precise measurements, however, the ordinary lens method may no 
longer compete so favourably with the semicircular one, since at higher resolving 
powers the use of ‘ central rays ’ in the electron-optical image formation becomes 
of increasing importance. ‘These rays have, as is well known, to be screened off 
by a central shutter in the lens case, but will always constitute the most valuable 
part of the beam in the semicircular case. It has been shown in earlier papers 
by two of us (Svartholm and Siegbahn 1946) that it is possible to combine the 
advantages of space focusing (high transmission) with that of using central rays 
(high resolving power). The basic idea was to use an inhomogeneous field of 
axial symmetry and of such shape that there is axial focusing as strong as the 
radial focusing. In the present paper we shall give, in § 1, a short review of the 
characteristic properties of this focusing method including some alternative 
subsidiary cases with reference to the special demands of the f-spectroscopic 
technique. ‘The aberration effect is treated in some detail in §2. In §3 the 
electron-optical arguments of the first section are faced by the specific problems 
of magnet construction, §4 gives a general description of the apparatus, while §5 
deals with the equipment for the magnetic field measurement to which some 
attention has been paid. A record of some of the adjustments and calibrations 
made is given in §6. Finally some results some of which have been published 


in other contexts, are reproduced in $7 in order to show some of the facilities of 
the instrument. 
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§1. THE ELECTRON-OPTICAL PROPERTIES OF THE 
TWO-DIRECTIONAL FOCUSING MAGNETIC FIELD 
According to the theory of the focusing properties of an inhomogeneous 
magnetic field of cylindrical symmetry, there is, in every field, an axial focusing 
action which decreases outwards from the axis of symmetry, and a fadial 
focusing action provided that it does not decrease more rapidly than the 
reciprocal of the distance from the axis. Consequently there is a two-directional 
focusing in every field or region of a field, which fulfils both these conditions. 
Or, more precisely, with a view to the application of this principle in a 
B-spectrometer, if the magnetic field has a symmetry plane perpendicular to its 
axis and if the axial component of the field in this plane at the distance r from the 
axis is denoted by H(r) and its derivative with respect to r by H’(r), then an 
electron leaving a point on the circle y=7, in the above-mentioned plane with 
a velocity vp corresponding to a momentum mv, =eryH(ry) and at a small angle 
to the circle, will first strike the cylindrical surface r =r, after an angle ¢,, the radial 
focusing angle, and the symmetry plane after an angle ¢,, the axial focusing angle. 
These focusing angles are determined by the equations : 


elas E fe ee | on (1) 


Pe ee Oe pee 
Hrs) | 


with the obvious relationship 


See 


Generally this focusing system is astigmatic, but putting ¢,=¢,, which gives 
; 1 
HI’ (ro) = — rm Rabe (3) 


makes the radial and axial focusing equal in strength and gives the focusing 
system all the characteristic properties of a Gaussian imaging system, the 
only difference being that the optical axis is a circle with a definite radius. 
Therefore, an extended object gives an inverted but true image in natural size 
after the angular distance 74/2=255° from the object. 

In so far as only central rays are considered the focusing conditions (1) 
and (2) are identical with the stability conditions in the theory of the betatron 
by Kerst and Serber (1941), cf. also Kurie, Osoba and Slack (1948). The 
condition (3) is then equivalent to the requirement of a field varying approxi- 
mately as 1/,\/r. As in our case, however, non-central rays must also be taken 
into account, we have to define the field shape more closely, for example, by 
means of a series expansion of the axial field component in the vicinity of the 


optical axis 
oe. ee 
H(0)=Ho| +a (“ =") + bh (‘ =) + ol Rr (4) 
1 


where the coefficient «, according to equation (3) has to be chosen equal to — 3 
while the other coefficients, such as «», are at our disposal. It was shown by one 
of the present authors (Svartholm 1946) and independently by McMillan (1946, 


962 A. Hedgran, K. Siegbahn and N. Svartholm 


private communication) and Shull and Dennison (1947) that the coefficient a, 
is the factor determining the aberration properties of the system and, conse- 
quently, the resolving power of the spectrometer. ‘Though the requirements 
of B-spectroscopy must not be identified a priori with those of faithful image 
formation, we may state that the spherical aberration is the most important 
factor, while other factors, which must be taken into consideration, are the sizes 
of the radioactive sample and of the slit of the G-M tube, astigmatism of the 
field, disturbances, and so on. 

Furthermore, the accessible resolving power of a /-spectrometer must 
always be considered in relation to its transmission power, a fact which must be 
remembered when different types of spectrometers are compared. It is also 
advisable for full utilization of any apparatus to balance some of the factors 
contributing to the resolving power, such as the spherical aberration, the size 
of the sample and the slit width of the G-M tube, in such a way that the trans- 
mission power, or more precisely, the counting rate in a definite momentum 
(or energy) range, will be as high as possible. Now the transmission power is 
usually defined as the solid angle in which the radiation from the sample is 
utilized in the apparatus, a definition which is adequate only in so far as the 
influence of the size of the sample upon the resolving power can be neglected. 
As this condition is not fulfilled in most applications, especially in investigations 
of primary B-spectra where the specimen must be made very thin and, therefore, 
large, it seems preferable to introduce a quantity defined as the solid angle 
multiplied by the area of the sample, and in accordance with a recent paper 
by Geoffrion (1949) on optimum conditions in the semicircular B-spectrometer 
we shall call this quantity the luminosity of the apparatus. 

The resolving power of a f-spectrometer is. defined empirically as the 
relative half-width of a f-line. ‘The possible contributions to the line width 
can be referred either to electron-optical imperfections, such as aberrations, 
or to effects inherent in the sample, such as the straggling effect. In a broad 
sense the latter effects constitute the natural line width of the sample. Since the 
latter may be highly different in different investigations it is desirable that the 
electron-optical line width (or the proper line width of the apparatus) be known 
for different arrangements of the shutters, slits and so on. This can be done 
either experimentally (by using samples with.a negligible natural line width) 
or theoretically. There are at least two reasons for making a careful study of this 
question. On the one hand, knowing the electron-optical line shape promotes 
greater accuracy in the energy determination of a f-line. On the other hand 
it is possible on the same basis to determine the natural line shape in order to 
investigate the influence of straggling and similar effects upon the primary 
radiation. 

As the electron-optical line shape may be very different in different cases it 
seems to be more practicable in a theoretical discussion to use the total relative 
line width instead of the former. The total relative line width is determined by 
the total image width 6 (in which the width of the G-M slit may be included in 
an incorrect sense) and the dispersion. ‘The length of the G-M slit is assumed 


to be so large that no radiation is lost. The dispersion d is expressed (Svartholm 
and Siegbahn 1946) by ; 


d =4r, Av/v, 
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which is twice as large as the corresponding quantity for a homogeneous field. 
This means that the total line width expressed for example by AH/H 1s 
expressed as 

AH/H =8/4ro. 


The proper image width 6 is composed of two terms at least, the object or sample 
width 59 and the aberration width 5, When the double focusing condition 
is not exactly fulfilled on the optica! axis a third term must be included, namely 
the astigmatic line width 8,,.. The aberration width 8,,, comes only from the 
spherical aberration. It is a characteristic property of the apparatus as a 
dispersive instrument with geometrical focusing that the object width enters the 
total line width to the first power while the aperture angles appear as the 
second power. 
The aberration width may be expressed (cf. Svartholm 1946) by 


Pcl V 0 cine 0 | a (5) 
where y, and y, are the radial and axial aperture angle, respectively. The 
numerical coefficients p and q are both related to the second derivative of the 
magnetic field according to the formulae 


In the region }<a,<% both p and g are positive. In the case of a square 
aperture y;=y.=y; we then obtain from equation (5) 


Sap =(P+ QF = 3Y'0» 
which is independent of the value of «. For this reason it is impossible to 
diminish the aberration line width by a suitable choice of «, when the radial 
and axial aperture angles are of the same order of magnitude. A considerable 
reduction of the aberration can be achieved, however, by a radially or axially 
extended aperture choosing «, in these cases equal to $ or 3 respectively. This 
question will be discussed in more detail later. 

The astigmatic line width 5,, derives from an incident deviation of the 
gradient of the field from its prescribed value. According to equation (3) the 
quantity «,, defined by 

dy =H (r9)/H(79); 


should have the value —4. If the relative deviation of «, from this value is e, 
so that «, = —4(1+<), then the path of a paraxial] ray in the symmeiry plane of 
the field is determined (Svartholm 1949) by 


1l+e)\-!? . 1+e\12 
r-n=ron (5) sin ( a ) ’ 


where ¢ is the angular distance from the object. In the image plane $=7/2 
we obtain in the first approximation 


r—1y=—Noyye7|V/2, 
which means that the astigmatic line width is 
Sy =ronnen V2 
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and that the intensity distribution within the line is the same as that within a 
line from an object of the same radial width. As a rule the astigmatism should 
be negligible, but in all cases where iron is used it is unavoidable that «, depends 
upon H in the region of weak fields. Moreover, if the sample and the G-M slit 
are situated on a circle of a radius 7,, which differs by Ar from the radius 7, of 
the optical axis, then there is an astigmatism of the first order expressed by 


 =(3—8ay)Ar/2r, , 


which vanishes if z,=. The meaning of this result is simply that there is only 
one field form, which gives double focusing for every value of 7, namely that 
varying as 1/,/r. In all other cases the double focusing condition is satisfied 
only for a special value of r (in the general case, however, for two different values 
of 7). 

While the question of the greatest utilizable aperture angle is intimately 
connected with the problem of magnet construction, which will be treated in $3, 


fs) 


1.33 


Q 


1 1.27 


Figure 1. The relation between solid angle () and total line width (8) for : (a) rectangular apertures, 
(6) circular apertures. 


it is useful to study to what extent the resolving power can be increased by means 
of shutters at the least possible cost of the aperture angle. In the cases 
% = and § where the aberration line width depends upon only one of the 
angles "1 and 5, it is evidently possible and advisable to diminish only the angle 
in question. In this way the solid aperture angle is varied linearly, while the 
line width is varied according to a square law. In our case where «,=% and 
the maximum values of y, and y. are equal (=0-2 radians), the relation between 
the relative line width 5 and the solid angle Q is, for radial shutters, 


§ =7?202/3y92 


where { is expressed in fractions of 47 (Figure 1). According to this formulae 
the line width is § per cent for the solid angle 4/7 per cent, but only 4 per cent 
for the angle 2/7 per cent. A similar result is obtained in the case %) =¢ if axial 
shutters are used. For other values of «, one needs both radial and axial shutters 
to obtain a solely linear dependence between 8 and Q. 


, 
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§2. THE THEORETICAL LINE SHAPE 


The electron-optical line shape, which will be discussed now, 1s defined as 
the number of electrons of a definite energy entering the slit of the G-M tube 
as a function of the magnetic field strength. Among the different contributions 
to the line shape only one has to be discussed in some detail, namely that due 
to aberration. 

When the line width is caused by more than one factor a composition of 
these is accomplished mathematically by the well-known ‘ Faltung’ process. 
If for example a finite object width gives in itself the intensity distribution J,(H) 


Figure 2. The aberration pattern of a point object and circular aperture in the cases: (@) %.=3/8, 


(b) 5/16, (c) 1/7, (d) 1/8. 


and if similarly the spherical aberration gives J,,(H) then the resulting intensity 
distribution is expressed by . 


I(H) = | I,(H —H')L,,(H')dH’. 


If the number of electrons emitted per unit area of the sample is constant so that 


Ti(iy=t5 for H,<H<H,+AH, 
and 1,(d2) =0 elsewhere, 
H,+ AH, 
then (HM) =) I,(H—-H" dH”. 
H, 


The function J,,(H) is calculated from the intensity distribution in the 
image plane of a point object emitting radiation isotropically. Since every 
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primary direction determined by the aperture angles y, and y, is represented 
by a point of the image in the image plane according to the formulae 


yt) = (pyr +o oe ae eee (8) 
B= = LOY V7. ene ee (9) 


it is possible to calculate the shape of the aberration pattern as a function of the 
aperture shape. ‘The representation is visualized schematically in Figure ae fOr 
four different values of the coefficient a. The marginal rays of a circular 
aperture y are indicated by two semicircles in the upper part of the Figure, one 
corresponding to rays inside, the other to rays outside the cylinder T=. In the 
general case both semicircles are represented by an ellipse with its centre in 
the point r=ry—3y79, x=0 and the semi-axes a=}|p—gq|y*7o and be |q| y?7o- 
The ellipse degenerates (cf. equations (6) and (7)) into a straight line in the two 
cases u, = % and «,=} and into a circle in the case «,=5/16. For all values of «, 
in the region }<a <2 the aberration pattern corresponding to the full circular 
aperture is limited by an arc of an ellipse and two tangents through the centre 
¥=1y 2=0. For a= the aberration pattern is a straight line of length $y,ro 


1/2 
v4 Z es 


a 


Figure 3. The aberration pattern of a point object for «.=%. The ellipse corresponds to a circular 
aperture (y), the other curves to rectangular apertures of different sizes (1y, $y, #y, y). Each 
surface element of this pattern is hit by the same number of electrons. 


for any aperture shape with a maximum radial aperture angle y,. For «,=4 the 
aberration pattern for a circular aperture is the ellipse itself while a rectangular 


aperture gives a pattern limited by an arc of the parabola r—7) = —32?/16y,?7_ © 


and the straight line r—r) = —4y,?79/3. ‘The pattern is plotted for four different 
values of y,; and y, in Figure 3. For a square aperture y,=y.=y the largest 
extent of the image in the axial direction is 16y?r,/3, which is twice that for a 
circular aperture with the angular radius y and four times the radial extent of the 
image. For reasonable values of ry and for example 7,=50cm., y=0-2, this 
figure is as large as llcm. ; this is the stronger argument against the case 
%2=% since the axial extent is a function of both angles y, and yp. 

Now we calculate from equations (8) and (9) a relation between the small 
solid angle dy,dy, and a surface element drdz of the image 


1 . 
dy,dy. = ae [g(r —19)*— pqz?}-¥? drdz. premio (LU) 
0 


As the number of electrons emitted within the angle dy,dy, is proportional to 
this angle the right-hand side of equation (10) gives the intensity distribution of 
the image. By integrating the expression (10) with respect to x between the 


4 
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margins of the pattern we obtain the number of electrons hitting a narrow slit 
of radial extent dr. This number, being a function of 7, represents the intensity 
distribution of a B-line measured by means of a G-M tube the slit of which is 
infinitely narrow in the radial direction but sufficiently long in the axial direction. 
The result of the integration will be reproduced only in some especially 
important cases. Thus in the case «= % the intensity distribution, when a 
rectangular aperture is used, is proportional to (ry7—r)-"? within the range 
¥y9>1>1)—4y,"ro/3 and zero elsewhere, while in the case ~,=} the distribution 
has the same shape within the range ry >r>r)—4y2270/3. If a circular aperture 
of the radius y is used, the distribution is 


1 ele 
ES ih BS Ea 


within the range ry >>r>1r)—4y?r9/3 and zero elsewhere. It is noteworthy that 
the distribution function for rectangular shutters is the same in the cases «, = 3 
and § as for the semicircular focusing system. Consequently the well-known 
theory of the composition of object width and aberration width can be transferred 
in full detail. 


(3) PRINCIPE bse ORSVEAGCNER CONSTRUCTION 

It is no doubt possible to realize a magnetic field possessing the required 
properties of inhomogeneity in a number of ways. Some general points can be 
stated, however, if one tries to reconcile the electron-optical demands defined 
in the preceding sections with the practical possibilities of magnet construction. 
As the dispersion is a linear function of the mean radius 7, it is desirable to choose 
the latter as large as possible. Now, measuring f-ray energies of 5-5 Mev. requires 
a value of Hp equal to 20,000 gausscm. Putting 7)=50cm. this gives 
HH, =400 gauss, which is a reasonable value for a magnet having a large pole gap. 
For a value of Hp 40 times smaller, corresponding to a f-energy of about 20 kev., 
the magnetic field strength is only 10 gauss. This is a just tolerable lower limit 
partly because of the remanence effect of the iron and partly because of the 
influence of the earth’s magnetism. Consequently a universally serviceable 
B-spectrometer should not have a mean radius greater than 50 cm., at least if iron 
is used. After having fixed the radius our next task is concerned with the 
dimensions of the pole gap. In our first small apparatus, which was used solely 
for demonstration of the double focusing principle, the magnetic field was 
realized by means of a ring magnet with a central cylindrical yoke surrounded 
by the coil. Though it is possible by using other designs for the yoke to utilize 
the space taken up by the coil for the effective pole gap, the ring magnet type 
was found to have some important advantages, especially the easy accessibility 
of the pole gap. As this arrangement prevents the use of a considerably large 
width of the pole gap, and consequently a large radial aperture angle y,, we 
preferred the case «, = % instead of a,=4. In that case the largest line width is 
determined only by y,, and equal to 4/3°% for y,; =0-2, which corresponds to an 
effective width of the pole gap 2a=28cm., which is computed from the formula 
a=ryV/2: The choice of the height of the pole gap 26, connected with the 
axial aperture angle y, by means of the relation b=rgy.\/2 is then mainly a 
question of power consumption. In other words, the solid aperture angle and 
the transmission power can be increased, the resolving power being fixed, by 


968 A. Hedgran, K. Siegbahn and N. Svartholm 


an increase of power. For larger pole distances, however, the fringing fields 
become more troublesome and require the use of an impracticably large amount 
of iron. For this reason we took y, =0:2, which corresponds to a mean distance 
between the pole faces equal to 28 cm. so that the pole gap profile is approximately 
a square. It should be admitted that resolution twice as good as ours can be 
obtained with a pole gap of the same size by choosing v=} and using circular 
shutters, but in that case we lose the possibility of improving the resolving power 
according to a square law by diminishing the radial aperture, as was expounded 
in §1. This possibility was considered advantageous for precise energy 
measurements. 


§4. GENERAL DESIGN 


The large spectrometer to be described in this paper has, compared with the 
first model, about ten times greater solid angle, this being about or somewhat 
more than 1% of 47. Furthermore the mean radius of curvature has been 
increased to 50cm. We believe that a further increase in radius and dimensions 
would be impracticable since the amount of iron which has to be put into such 
an apparatus rapidly limits these possibilities. 

Generally the resolution of a spectrometer is given by a sum of two terms, 
one depending on the solid angle used and the other on the sizes of the sample 
and the G-M slit compared with the radius of curvature. ‘Thus by making the 
spectrometer large the important advantage is obtained that more extended 
sources can be used at a given resolution. When working at very good resolution 
(<1%), for which the new spectrometer is intended, this condition must be 
regarded as rather essential since in a great number of cases the maximum 
obtainable resolution is limited by the thickness of the radioactive sample itself 
due to straggling. ‘This is especially obvious when investigations of secondary 
photo electrons expelled from a lead converter are made. The number of photo 
electrons increases when the amount of lead is increased. On the other hand 
the photo lines soon get broader than approximately 1° when the thickness 
of the converter exceeds a few mg/cm? of lead. It is therefore important to spread 
out the lead on a large, thin surface (by evaporation) which in turn requires a big 
spectrometer if the geometrical resolution due to the large surface is not to be 
decreased too much. Furthermore, when the dimensions of the spectrometer 
are increased, the G-M counter can be more effectively shielded from stray 
radiation and y-radiation from the source. 


The Magnet 

If the iron permeability is put equal to infinity the surfaces of the pole shoes 
are equipotentials. In this first approximation the airgap is very nearly a parabola 
Z=C v7, i.e. with the vertex at the centre of the magnet. One has to take into 
account the effects of finite iron permeability and also the outer and inner fringing 
fluxes. The magnitudes of these corrections are largely dependent on the type 
of magnet chosen. As was mentioned in the preceding section we preferred an 
alternative shown in the lower part of Figure 4, where the airgap surrounds the 
coil as well as the iron yoke.* Another alternative shown in the upper part of 
Figure 4 proposed to us by Dr. L. Dreyfus at ASEA has the advantage that 


* We are indebted to Dr. Dreyfus for the design of the magnet. 
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the fringing fluxes are more easy to handle. This alternative requires, however, 
more copper and power, but seems to have definite advantages, for example 
in the case of smaller magnets. 

It Hp = 20,000 is put as an upper limit for any 6-spectrum to be investigated, 
it 1s immediately realized that the magnetic fields employed at p=50cm. are 
rather small (less than 400 gauss, usually of the order of 50-100 gauss or still 
smaller). ‘The problems of remanence and of changing permeability at the lower 
part of the magnetization curve have also to be considered. It is therefore 
impossible to base Hp determinations on any approximate linear relationship 
between current and magnetic field since the departure from this linearity is 
far too great, at least in the lower energy region. It is therefore necessary to be 
able to measure the magnetic field directly and for convenience this can be done 
continuously in the present apparatus. The system adopted for this purpose 
will be described later on in this paper. 


MMU YY 
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Figure 4. Two alternatives for the magnet construction. 
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The magnet weighs about 4 tons and consists of high permeability iron. It is 
constructed in three parts : a central cylindrical nucleus (6 =45 cm., h=20cm.), 
and two poleplates (¢=154cm.). The pole gap at the mean radius (p =50.cm.) 
is 28cm. and it is possible to utilize a field gap which does not differ very much 
from this figure in the radial direction. At maximum transmission the cross 
section of the transmitted electron beam is therefore at least approximately a 
square. Since the type of field chosen in this apparatus has the property of 
focusing electrons in the axial direction to a higher order of approximation than 
in the radial direction, the cross section of the transmitted beam at better 
resolution is made rectangular by keeping the height constant and decreasing 
the radial aperture. Thus, at good resolution a departure from the required 
exact field form in the outer parts of the field is comparatively unimportant. 
The magnet is excited by four coils in series having a total resistance of 16 Q. 
The number of turns is 2,000 and the copper wire has a cross section of 2 x 4mm. 
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As the shape of the pole shoes was computed for a maximum field of 400 gauss 
the field form must be expected to deviate from the computed one at very much 
lower fields, but this deviation can to a great extent be compensated for by a 
suitable premagnetization. A deviation in the field gradient from the correct 
one at lower fields causes only a very small shift of the top of the lines, but the 
intensities can be affected more seriously. If good intensity measurements 
are to be performed, especially at low energies, the behaviour of the field gradient 
must be carefully controlled. 

The measured magnetic field distribution is shown in Figure 5. These 
measurements were performed with a fluxmeter and search coil at fields 400, 
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Figure 5. Magnetic field distribution. 


128 and 36 gauss. The fields corresponding to «, =} and 3 are drawn in the Figure 
for comparison. ‘The two field distributions at 400, 128 and 36 gauss are, as can be 
seen almost identical and come very close to the 3 case, especially for r<60 cm. 
where the deviation is generally less than 0:2°%, considered a reasonable tolerance. 

If the magnet is not demagnetized, there is a tendency at lower fields to 
form a slightly more uniform field. 


Design of the Spectrometer 

The interior of the spectrometer can be seen in Figure 6 (Plate). Two 15mm. 
thick brass cylinders form the outer and inner walls of the vacuum tank. The 
pole shoes of the magnet constitute the roof and bottom of it and are pressed 
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against rubber gaskets (cross section 4 mm. x 6mm.) running in grooves along the 
edges of the brass cylinders. 

The sample holder can be seen to the right in the figure. The sample is 
introduced into the vacuum chamber by means of a vacuum lock device. 

There are three shutters. he first of these is permanently set for a little 
more than the intended maximum transmission of the spectrometer. The 
other two shutters are adjustable from outside by means of Wilson sealing 
arrangements and the position of their two vertical defining slits can be read on 
a millimetre scale. Since the focusing in the axial direction, as mentioned before, 
is attained to a higher approximation than in the radial direction, the upwards 
and downwards defining slits for our field can be set for maximum transmission 
in that direction (defined by the pole distance), and need not be adjusted when 
better resolving power is wanted (at least to a certain limit). 

The G-M counter (of bell type) is mounted on an adjustable brass holder 
and different sized G-M slits can be moved in front of the counter (width = 3, 6, 
12, 24mm.). 

Large lead blocks (not shown in the figure) are placed between the source 
and the G-M counter, and because of the large distance between these units 
even very strong sources can be used without increasing the background (about 
20 impulses/minute) of the counter to any noticeable amount. 

The vacuum chamber (approximately 300 litres) is evacuated by means of a 
big molecular disc pump (200 l/sec.) combined with a Cenco rough vacuum 
pump, and the pressure (~15 x10-°mm.) can be controlled by an ordinary 
ionization vacuum gauge, stabilized for constant emission. ‘The lock device for 
the sample can be separately connected to the rough pump, so that the high 
vacuum is not destroyed when the sample is introduced into the spectrometer. 

The magnetic field measuring equipment, which will be described below, 
contains a rotating coil system on a long rod. This rod enters the spectrometer 
through a brass tube, closed at the inner end and screwed on to the spectro- 
meter wall. (It is rendered vacuum-tight by means of rubber gaskets.) ‘This 
tube can be seen in the front of Figure 6. 

Figure 7 (Plate) shows a photograph of the total apparatus. 


§5. MAGNETIC FIELD-MEASURING EQUIPMENT 

The large dimensions of the spectrometer with the correspondingly great 
dispersion make the instrument well suited for precise energy measurements of 
-lines and photo electron lines. ‘The resolution or relative haif-width of the 
lines is in most cases limited to about 1°% or a little less. This limit is not set 
because of the focusing properties of the instrument itself, but is due to 8-source 
thickness and converter thickness in practical cases. ‘Though it doubtlessly 
may be possible in many cases to decrease the * natural ’ line widths caused by 
straggling in the source and converter respectively, the experience hitherto 
gained clearly demonstrates that a number of difficulties arise in connection with 
a full utilization of still smaller line widths. Some of these will be discussed later. 

If the relative line width (measured at half maximum height of the line) is 1% 
it ought to be possible to locate the position of the line within about 0-1%. ‘This 
means that the magnetic field has to be measured with the same accuracy. 
Furthermore the field has to be measured for each setting and controlled 
continuously. This, of course, necessitates that the measuring system should be 
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easy to handle. Ordinary ballistic methods are therefore unsuitable. Various 
of other methods already in existence which allow continuous readings are either 
definitely not accurate enough or else too uncertain to be relied upon for longer 
and repeated runs. We have finally adopted a balance system (null method) 
making use of a reference field produced in a Helmholtz coil, piaced at a distance 
from the spectrometer. 

Two coils, one in the spectrometer and the other in the reference field, rotate 
on the same axis. The magnetic field in the spectrometer is adjusted until the 
E.M.F. induced in the rotating coil in the spectrometer is equal to that induced 
in the corresponding coil in the Helmholtz field. By carefully adjusting the 
rotating coils relative to the two fields once and for all the phases of the E.M.F. 
induced in the two coils can be made very nearly opposite. When the spectro- 
meter field is changed until minimum alternating current flows in the rotating 
system, the current through the Helmholtz coil is a measure of the magnetic 
field in the spectrometer. The calibration of the arrangement is most 
conveniently made by means of a strong f-line of known energy (e.g. ThB F 
line etc.). 

The actual arrangement is shown schematically in Figure 8. The two 
rotating coils (A and B) are situated on a stainless steel rod 2 metres long which 
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Figure 8. Arrangement for measuring the magnetic field. 


has been accurately straightened. A third coil (C,) on the steel rod rotates around 
its own axis and is connected to an amplifier and a microammeter. The coils 
A and B (6=30mm.) contain each 1,000 turns 0-2mm. copper wire, the coils 
C, (¢=44mm.) 2,000 turns 0-2mm., and C, (6=55 mm.) 4,000 turns 0-1 mm. 
respectively. The rod rotates at a speed of about 30 revolutions per second 
and is driven by means of a motor. The arrangement does not require a 
particularly constant speed of rotation. 

In order to obtain sharp minima for the field determinations the rotating 
system has to be very carefully balanced mechanically so that no extra vibrations 
etc. occur. The rotating coils should be adjusted to rotate around a diameter. 
The sharpness of the minimum is also critically dependent on the adjustment 
of the relative phases of the two induced E.M.r’s. The final adjustments in 
this respect are perfected by means of screws supporting the Helmholtz coils 
on the rigid foundation. 

. Figure 9 shows the amplified current through the pick-up coil system on the 
microammeter, viz. current through the Helmholtz coil, the spectrometer field 
being kept constant (25 gauss). According to this curve the position of the 
minimum can be defined with an accuracy within about + Ost 
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Figure 6. The interior of the spectrometer. 
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Figure 7. Photograph of the spectrometer. 
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Most of the experiments so far made with the spectrometer have been 
performed with the system described above. It should be mentioned here, 
however, that recently this system has been rebuilt. The mechanical stability 
is increased and the definition of the minimum is now around 0-01057 the 
description of this improved system will soon be published, when more 
experience has been gained with it. 

The rotating rod is situated in the north-south direction and the earth’s 
magnetic field then contributes by its vertical component (Ie) 40) thet Mie: 
produced by the field in the Helmholtz coil. The magnitude of this component 
has been measured and must be considered when evaluating the Hp values 


is deflection 


——_ ~ ———————— 
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Figure 9. The minimum of the induced E.M.F. in the rotating system, showing the value of the 
magnetic field. 


for different values of J through the Helmholtz coil. The relation (Hp, J) 
(amps through Helmholtz coil) is then expressed by 


Hp (gauss cm.) =constant x /+pH,,..4 


(~50 x 0-4). The last term is, however, always very small compared with the 
first one. 

It is important that the coil rotating in the spectrometer is actually placed 
at the same distance from the centre of the magnet as the sample and detector 
(optical axis of the spectrometer). Small changes in the magnetic gradient 
which, as stated above, by themselves cause no noticeable shifts of the line peak, 
will of course introduce errors in the relative field determinations if the coil is 
not placed on the mean radius used (p =50 cm.). ° ; 

Since the electric power required to focus even high energy f-particles with 
the spectrometer is quite low (for instance 2 Mev. corresponds to about 
200 watts) one can very conveniently use batteries as current supply. ‘Thus no 
extra current-stabilizing system is needed. 

The current and magnetic field regulating system is shown schematically in 
Figure 10. The principle of the arrangement is the following. A certain current 
through the Helmholtz coil is first chosen (regulated by means of the resistance R,) 
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corresponding to the desired field in the spectrometer. This current branches 
into two, one branch exciting the magnet. The Helmholtz coil is so designed 
that it will produce a field equivalent to half that in the spectrometer at the same 
current. By means of the small resistance R, the ratio between the two currents 
in the bridge can be varied, until the magnetic field in the spectrometer just 
balances the Helmholtz field (as indicated by the rotating system). This 
adjustment will not affect the total current through the Helmholtz coil in the 
first approximation, provided the resistances of the two branches are approxi- 
mately equal, and R, is small compared with R;. ‘These conditions are always 
satisfied in our case for high fields. At low fields, when the deviation from linear 
relationship between current and magnetic field is relatively greater, the bridge 
has to be side-balanced more than the small resistance R, would allow. This is 


Figure 10. The current-regulating system. 


accomplished by the multipole switch system, shown in the Figure. By means of 
this the current through the spectrometer can be changed in bigger jumps 
simultaneously with the introduction of new resistances into the main line, so 
that the total current is kept constant. When this switch is far from its zero 
position the first of the above conditions is violated. However, as mentioned 
before, this only happeris for low fields, where the resistance R, is very high 
compared with the resistance in the bridge. In that case the final balancing 
adjustments by the resistance R, will not affect the main current either. The 
practical experience gained with the system here described has shown that one 
never needs to readjust the main current (giving the desired field) after balancing 
the bridge, which simplifies the field setting very much. It should also be 
remarked that within smaller field intervals, such as that occupied by a single 
line etc., the balancing does not need to be changed. 
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§6. ADJUSTMENTS AND CALIBRATION OF THE SPECTROMETER 
After the magnetic gradient has been measured and a convenient mean 
radius (radius of the optical axis) chosen, the positions of the sample and 
detector must be experimentally adjusted relative to each other. ‘This is most 
easily accomplished by means of the photographic method. ‘The exposure of a 
strong f-line (the source being deposited on an aluminium cross) will immediately 
tell if the image falls on the correct height relative to the G-M slit. Simultaneously 

the axial focusing properties are especially well demonstrated on such plates. 
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i 1. The H,I,JandJ,-linesof Figure 12. The X and X,-lines of ThC’” (£y=2-62 Mev.) 
Ts Width of eaeplees 5mm. Width of sample=3 mm., shutter width, left 10°5 cm., 
Shutter widths=10-5cm. right 65cm. G-Mslit, left 12 mm., right 6 mm. 


G-M slit=6 mm. 

The choice of shutters and width of the G-M slit has to be investigated for 
different widths of the samples and the resulting resolving powers, as obtained 
from the half-widths of lines, registered. These investigations must be 
performed by means of weightless samples such as ThB + products since otherwise 
the line widths may be strongly dependent on the weight of samples used, energy 


of f-lines etc., due to straggling in the source itself. Figures 11, 12 and 13 
66-2 
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show a number of f-lines obtained by the spectrometer with different shutters, 
G-M slit width and sample width. The height of the sample and G-M slit was 
kept constant (20mm.). Decreasing the line width below 0-35°% was not tried 
since the field measuring equipment used in these experiments limited such 
possibilities. 

The spectrometer was calibrated using seven different B-lines in the spectrum 
of ThB + products, namely (according to Ellis’ notation) : F, H, I, J, Ja, X, Xa 
lines. These lines cover an energy range from 150 kev. to 2-6mev. The relative 
Hp values of these lines are probably accurate within about 0-1°% according to 
Ellis’ (1934) measurements. The latest absolute determination (Siegbahn 1944) 
of Hp values in this spectrum concerns only the F and I lines and provides the 
values Hp, =1,383-8 and Hp, =1,749-6. These values (based on x-ray data) differ 


P 
x X 
7 os 5 See ai. 3.5 Lee 
° 
—— = T eS al 
42 43 44 amp. 42 43 44 amp. 


Figure 13. The X-line of ThC’”’. Width of sample=3 mm., G-M slit=3 mm., shutter width 
left 10-5 cm., right 6-5 cm. i 


only by 0-15°% from Ellis’ corresponding values (Hp, =1,385-8 and Hp, =1,751-0) 
and were estimated to have probable errors around 0:03°%%. This B-spectrum is 
certainly well suited for calibration purposes. ‘The constant before I in the 
calibration formula of the spectrometer (Hp =constant x +20) was determined 
for these seven lines and the results are shown in Table 1. 


Table 1, Calibration : Hp =constant x [+20 


Line I (amp. at peak) Constant Line I (amp. at peak) Constant 
F 0-590(0) 2313 ie 0:774(0) 23 
; ws 
H 0-719(0) 2320 xX 4-31(0) 2312 
I 0:745(5) 2315 Xs 4-42(0) DENG 
J 0-769(5) 2319 Mean : 2317 


The maximum deviation in any line measurement from the mean value is 
about 0-2%. Furthermore there is no sign of any systematic change in the value 


’ 
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of the constant as a function of the 6-energy. The maximum deviation of 0:2°% 
must be expected since the ammeter used in these experiments (Western Electric 
Instr. Corp. Model 5, No. 1695) to measure the current through the Helmholtz 
coil has a precision of within +0-1°, (calibrated recently at our Institute) and 
the field balancing system was able to give the same accuracy. With the new 
balancing system and a compensation arrangement to measure the current 
there are good possibilities of increasing the precision still more in the 
future. 

On this point it should be stressed, however, that in ordinary cases, i.e. when 
the weight of the sample is not negligible, there are great difficulties which 
prevent full use being made of much better accuracy than that already obtained 
(within 0-1—0-2%). The electron straggling in the sample which is energy 
dependent will not only cause a line broadening, but will also shift the lines 
towards lower energies. To handle this effect and make accurate, well 
established corrections for it is an intricate problem, which calls for a systematic 
experimental investigation in the future, in order to bring down the systematic 
errors in 6-energy determinations to the region of 0:1% or less. Such investiga- 
tions are in progress in this laboratory and valuable contributions have recently 
been supplied by Jensen et al. (1949) and Lauritsen et al. (1949). The straggling 
effect is perhaps most obviously recognized when using converters to obtain 
external photo lines from the y-rays. In these cases one must of necessity use 
radiators of finite thickness in order to get measurable intensities of the photo 
lines. It is found that line shifts due to straggling may very well be of the order 
of 1% even if the radiator is less than 10u. Furthermore the relative line shifts 
are not only dependent on the y-energy, but also on the actual choice of spectro- 
meter resolution etc. (Lauritsen and Hornyak 1949). The most straightforward 
procedure to investigate this problem is to measure first a number of /-lines 
from ‘ThB + products (weightless sample) and then measure the position of the 
corresponding external photo lines from a converter in a given geometry. 
Another valuable method, which has been suggested by DuMond et al. (1948) 
would be to determine the energy of a number of different y-rays from wavelength 
measurements by x-ray methods and then to calibrate the converter directly by 
observing the position of the corresponding photo lines. Hitherto only a few y-ray 
wavelengths are known with sufficient accuracy to serve this purpose and there 
are some difficulties in measuring hard y-rays (i.e. >1 Mev.) accurately with this 
method. The half-widths hitherto obtained by DuMond are around 1%, but 
the wavelengths are estimated to be considerably more precisely known (for the 
annihilation radiation to approximately 1 part in 10*). 

Before these questions are completely mastered the uncertainty in the energy 
determinations from photo lines or from f-lines of non-weightless samples 
generally have to be given with larger errors than those set by the half-widths of 
the lines and the field measuring system themselves. 

In those cases where the energy of the photo line comes very close to the 
energy of the photo lines of the annihilation radiation the energy dependence 
of the line shift may be disregarded and a direct comparison can be made between 
the unknown y-line and the energy of annihilation radiation. If the same 
converter is used in both cases or, better, if the two samples are simultaneously 
put in the same converter, one can measure a close ‘ doublet’ of photo lines, 
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one of which is very well known (Eynnin, a. =510°8kev.). ‘This procedure was 
tried by Siegbahn and Hedgran (1949) with the spectrometer in the case of *Au, 
where the K photo line of the annihilation radiation and the L photo line of the 
gold y-radiation form a doublet consisting of two lines only 28 kev. apart from 
each other, when using lead converter (the doublet would be still closer if uranium 
instead of lead had been used). In this way the y-radiation was determined to be 
411-0 + 1-0kev., in good agreement with DuMond’s corresponding crystal value 
(411-2kev.). It is of course possible to extend this procedure to other energies 
by comparing the K photo line of gold against an L photo line of some other 
y-ray using the same converter etc. The doublet method can certainly be 
made more accurate in future with the new field-measuring system. 


§7. SOME RESULTS 
The spectrometer has been used for some time to study different radioactive 
disintegration problems. Some of these investigations have already been 
published or will appear elsewhere. In most cases the spectrometer has been 
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Figure 14. 14Cs photo-electron spectrum. Lead converter 9m. 


adjusted to give a resolving power of about 1% or slightly less. To make full 
use of this resolving power it has been found necessary to restrict the thickness 
of the lead converter in y-ray investigations in most cases to be less than 5 
(evaporated on an aluminium foil). The heights of the photo lines compared 
with the inevitable Compton electron distributions are then not particularly 

[Continued on p. 984 
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(a) (6) 
Figure 15. 1%4Cs B-spectrum. 


(a) With 1% resolution in the spectrometer. 
(b) With 4:5% resolution taken with another spectrometer. 
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Figure 16. ?°%Hg photo-electron spectrum. Lead converter 3-7. 
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Figure 17. +*8Au photo-electron spectrum. Lead converter 3°7 p. 
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Figure 18. Zn photo-electron spectrum. Lead converter 9 b 
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‘{a) 
Figure 23. 5®Co photo-electron spectrum. 
(a) Lead converter ~50 4, spectrometer resolution ~1%. 
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Figure 20. Co photo-electron spectrum. 
(b) Lead converter 3-7, spectrometer resolution ~1%. 
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Figure 21. K, Land M photo lines from the annihilation radiation (Ey=510°8 kev.). 
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Figure 23. 1°Ag photo-electron spectrum. Lead converter 3°7 p. 
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Figure 24. Internal conversion lines of ®°Co. 
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Figure 25. Internal conversion lines of ®Br. 


great despite the high spectrometer resolution, but instead the lines are sharp 
and the line shifts due to straggling relatively small. In some cases the resolution 
has made it possible to resolve close lines from each other, which with thicker 
converters and with less spectrometer resolution would have been taken to be a 
single line. 

“*Cs has three y-lines (plus a hard one of very low intensity at 1-35 MeV.) 
at 567, 603 and 796kev. (Elliot and Bell 1947, Siegbahn and Deutsch 1948). 
The two softer lines differ by 45° in Hp. These two lines are completely 
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resolved in the photo-electron spectrum shown in Figure 14 (Pb converter Dw). 
The conversion factors for the y-rays are low, as expected (4-4 x 10-3) for the 
two softer y-rays and 2-4 x 10-° for the one at 796kev. The increased resolution 
makes these f-lines easy to detect in the primary f-spectrum (Figure 15(a)). 
For comparison the same spectrum is shown (Figure 15(b)) obtained with a 
spectrometer adjusted for a resolution of 45°. (In that case the two softer B 
and photo lines could not be resolved.) 

Some other typical photo-electron spectra (with increasing energies) are 
shown in Figures 16, 17, 18 and 19 (y-rays from 2%Hg (Slatis and Siegbahn 1949), 
Au (Siegbahn and Hedgran 1949), Zn and **Fe respectively). The corre- 
sponding y-energies are °3Hg =279 kev., Au =411 kev., ®Zn =1-125 Mev. (plus 
a few presenting annihilation radiation * and *°Fe = 1-097 and 1-295 Mev. Except 
for 1*8Au, where the doublet method has been used, the maximum errors of the 
y-energies should not yet be set at less than 0-5%. 

Figure 20 is a good illustration of the importance of using a converter thickness 
suitable for good spectrometer resolution. Figure 20(a) shows the photo- 
electron spectrum of **Co (Deutsch and Hedgran 1949) with a lead converter 
approximately 50u thick and the spectrometer set for a resolution of about 1°%. 
Figure 20(b) shows the same spectrum with the same spectrometer resolution, 
but with a lead converter 3-7u thick. In the last case a K and an L photo line, 
previously masked by the strong annihilation radiation photo lines, can now 
be observed. ‘The existence of this line is essential for the disintegration scheme 
of ®Co. Figure 21 (with a 3-7y lead converter) shows the K, L and M photo- 
electron lines from the annihilation radiation used to calibrate the spectrometer 
for photo-electron spectra around these energies. ‘The calibration factors obtained 
from the three different lines differ among themselves by less than 0:2%. 

Figures 22 and 23 show part of the primary f-spectrum of °Ag (Siegbahn 
1949) and the corresponding photo and Compton electrons from a 3:7 lead 
converter respectively. As can be seen this spectrum is unusually complex and 
good resolution is necessary if any analysis of the radiation is to be tried. 

According to theory the internal conversion coefficients at low Z values and 
hard y-rays must be expected to be very low. Comparison with theory in this 
region is therefore difficult. In the case where the f-lines fall outside the 
continuous f-spectrum the situation is however much simplified. ‘This happens 
for instance in the case of ®°Co (Deutsch and Siegbahn 1950) and **Br (Siegbahn, 
Hedgran and Deutsch 1949). These two spectra are shown in Figures 24 and 25. 
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ABSTRACT. A post-focusing electron accelerator has been developed to reduce the 
absorption of low-energy electrons by the Geiger counter window of a f-ray lens spectro- 
meter. Its useful range is for B-rays having energies between 0 and 30 kev., and although, 
in general, for a portion of this range only partial transmission will be obtained, correction 
curves for the absorption can be accurately determined. 

In order to test the spectrometer, an investigation of the L-Auger lines of thorium-active 
deposit has been made. ‘An identification of the lines has been attempted. 


$1, INTRODUCTION 


O improve the transmission of a magnetic lens spectrometer* for low energy 

electrons a post-focusing accelerator has been fitted to the spectrometer 

tube. An account of a similar apparatus has already been published by 
Langer and Cook (1948); construction of the present apparatus was, however, 
carried out quite independently of this publication, and the tests to which it has 
been subjected are possibly more severe than those applied by Langer and Cook 
to their apparatus. 

The function of the accelerator is merely to accelerate the electrons, after 
they have been focused, to give them additional energy so that they may penetrate 
the Geiger counter window. A diagram of the apparatus is shown in Figure 1. 

The electrons are focused through a hole A, just over 2mm. in diameter. 
After passing through this hole, they are accelerated towards the Geiger counter 
window B. In the present case the window consists of collodion foil, approxi- 
mately 0-4 mg/cm? in thickness, mounted on a brass grid. As will be described in 
more detail later, the accelerating potential normally used is about 18 kv. 

Since the Geiger counter is, of necessity, at this high positive potential, it has 
been found most convenient to raise the counter scaling circuit and power pack 
to the same potential. The general layout of the necessary apparatus is shown in 


* Details of the constants of this spectrometer have been given in an earlier paper (Butt 1949). 
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Figure 2. The scaler circuit and power pack are housed in a bin lying horizontally 
on four insulators. The a.c. supply to these units is insulated from the a.c. 
mains by means of a 1: 1 ratio transformer, the primary of which is insulated from 
the secondary to a potential of 100kv. The dials of the scaler are observed 
through windows in the side of the bin, and the controls are operated by means of 
ebonite rods protruding from the side of the bin. 
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The accelerating voltage is supplied by a 100kv. x-ray transformer, the 
output of which can be continuously controlled over the whole of its range. 
This voltage is rectified by means of a half-wave valve rectifier and smoothed by 
a 0:02 mr. condenser. 

During operation, the pressure in the accelerating tube is kept below 
10-> mm. Hg. 


§2. PERFORMANCE OF THE ACCELERATOR 

When preliminary tests of the accelerator were made, the following character- 
istics were recorded. Initially, with no source in the spectrometer, the spectro- 
meter was evacuated to a pressure of less that 10-° mm. Hg and the Geiger counter 
was filled. The accelerating voltage was gradually raised from zero. No 
spurious counts were observed until the accelerating potential was approximately 
19kv. Above this voltage, however, the number of such counts increased very 
rapidly until, at about 20kv., the scaling circuit was completely blocked. 

To determine whether these spurious counts were due to radiation actually 
entering the counter, a 20 mg/cm? copper foil was placed in front of the counter 
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window. On the application of the acceleration voltage no spurious counts 
were recorded. This points to the fact that these spurious counts are caused by 
electrons or ultra-violet radiation. 

That the spurious counts are due to electrons was finally verified by deflecting 
them away from the counter by means of a magnetic field. These spurious 
counts are presumably due to the formation of positive ions in the accelerating 
chamber which bombard the cathode and cause it to eject electrons; these in 
turn are accelerated to the Geiger counter window. From this theory it is to be 
expected that spurious counts will begin to occur when the ejected electrons are 
accelerated through the Geiger counter window. The fact that, in the present 
case, spurious counts begin to occur at 18 kv. roughly supports this idea, because 
it is to be expected that a 0-4mg/cm? window will begin to transmit when the 
electrons reach an energy of about 18 kev. 

An attempt to reduce these spurious counts by using an aluminium cathode 
and then coating it with ‘Aquadag’, in order to raise the work function of the 
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cathode surface and lower the stopping power to the positive ions, made very little 
difference. A natural limitation of the useful accelerating voltage (at about 
18kv.) was therefore accepted. Since the collodion windows can transmit 
electrons with energies as low as 18kyv., with the aid of the accelerator, a small 
percentage of electrons having energies of only 1 or 2kev. can be detected. 

An advantage of the method of acceleration is that the transmission curve for 
any window can be quickly and accurately determined. A source is put into the 
apparatus, and the lens is made to focus electrons. of a certain momentum. The 
accelerating voltage is then applied and is increased in steps of 1 kv., counting being 
continued until a saturation counting rate is reached. At this stage all the electrons 
incident on the window are transmitted. From one such curve the transmission 
of the film for electrons of any energy may obviously be determined. Such a 
curve is shown in Figure 3 for one of the windows used in the present work. 
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§3. INVESTIGATIONS WITH THORIUM-ACTIVE DEPOSIT 


In order to test the behaviour of the apparatus, the low-energy line spectrum 
of Th(B+C+C’+C”") was investigated. The active deposit was carried on a 
thick platinum button just under 2mm. in diameter. Back-scattering from the 
thick button, of course, made the continuous spectrum in the low-energy range 
almost entirely ‘instrumental’, but it was thought that the relative intensities of 
the lines in the same region should not have been greatly affected. 

To test the possibility of the electric field leaking through the hole A (Figure 1) 
into the spectrometer chamber, and so vitiating the calculation of line energies, 
the low-energy spectrum was observed with an accelerating potential of 10kv. 
and then with an accelerating potential of 18kv. Calculated line energies 
were the same within acceptable limits. 

There is a further effect which might lead to the distortion of the line energies, 
at very low energies. ‘This arises from the fact that the axis of the spectrometer 
is horizontal and only the vertical component of the earth’s magnetic field is 
compensated. ‘The axial focusing field, therefore, is not strictly proportional 
to the focusing current. ‘To test whether this was affecting the calibration of the 
spectrometer in any way, the line energies were observed with the focusing current 
passing first in one direction and then in the other. The energies agreed to within 
1°% within the energy range 6-12 kev., which is about the limit of accuracy with 
which the energies can be determined. 

The low-energy spectrum finally obtained is shown in Figure 4. As will be 
explained, the lines in this spectrum are L-Auger lines of ThC(83) and ThC’ (81). 


Counting Rate 


250 300 350 400 
Ho (gauss cm) 


Figure 4, 


Flammersfeld (1939) published some low energy lines in the Th(B+C+C’+C") 
spectrum that he presumed to be new. He reported four lines of Hp values 284, 
314, 333 and 368 gauss cm. Flammersfeld hinted that these lines were possibly 
due to Auger electrons, but he made no effort to identify the transitions concerned. 

These lines of Flammersfeld are presumably the same as those investigated 
photographically by Black (1925). Strangely enough, Flammersfeld’s energies 
do not appear to agree very well with those of Black ; for instance, his most intense 
line is that of Hp =284, whilst Black’s strongest line is of Hp =303-4. Flam- 
mersfeld found no line near Hp = 303-4. 
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The genera! shape of the spectrum obtained in Figure 4 closely resembles 
that obtained by Flammersfeld, though there is considerably more detail present 
in Figure 4 than there is in Flammersfeld’s spectrum. Black’s spectrum also 
has more lines than Flammersfeld’s and, as the following table shows, our energies 
agree with Black’s in nearly every case. Flammersfeld’s energies are included 
for comparison. All values are in kev. 


Table 1 
Flammersfeld Black Author Flammersfeld Black Author 
6:46 6°71 a® 8-7 9-82 9-82e 
7-06 7-10 b 9-7 10:45 10°39 £ 
7°53 7:°55+0:05 c 10-45 12-44 12:26 g 
7:00 8-11 8:10-+0:05 d 12-69 12:62 h 


8°47 


* The lettering here given to the lines is the author’s own nomenclature. 


Relative intensities are omitted, since Black’s are merely visual estimates, and 
it is obvious from Figure 4 that ours cannot be determined with any high degree 
of accuracy. In respect of energies, our lines a and g diverge slightly from 
Black’s values, and line c is one that he suspected as being present but was unable 
to confirm. Black also queried the genuineness of the 310-1 Hp line (8-47 kev.), 
and Figure 4 shows no sign of it. 

In view of the general agreement between the last two columns of Table 1, 
Flammersfeld’s energy values are rather difficult to explain. One possible 
reason why he obtained lower energies is that he deposited his source on a thin 
Zapon film. ‘The resultant charge lost by such a source is negative, and Zapon 
film in a vacuum is a very good insulator. Braden and his colleagues (1948) 
have shown that in such circumstances the source may become positively charged, 
so that electrons leaving the source lose energy in doing so. ‘The fact that both 
Black’s source and the author’s source were deposited on a conductive platinum 
backing is perhaps significant in this connection. On the other hand Flammersfeld’s 
spectrum is not uniformly lower in energy than ours, which might be expected 
from a simple explanation along these lines. 
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In attempting an identification of the lines listed in Table 1, it is first necessary 
to determine from which elements the lines are likely to have originated. The 
elements in equilibrium on the source have Z=81, 82, 83 and 84. Clearly, only 
those giving strong groups of secondary internal conversion electrons need be 
considered as sources of tertiary Auger electrons. 

From Flammersfeld’s tables it can be seen that by far the most intense 
internal conversion occurs in the K and L levels of element 83 and in the L levels 
of element 81. For the present purpose the contributions of elements 82 and 84 
can be neglected. We shall therefore proceed to estimate the ultimate intensity 
of L shell ionization in Z=83, 81 and its distribution amongst L,, Ly, and Ly, 
as a basis for our discussion. 


Lie 83: 


It is well known that strong internal conversion takes place in the K and L 
shells of this element due to the Dg, E, F, Fa and H y-rays (nomenclature of Ellis). 
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These y-rays produce the K internal conversion lines D, A’*, F and H, which give 
a total K ionization of 33-12 per 100 disintegrations. 

Again, according to Flammersfeld, the total L shell ionization (for all three 
shells) due to the complex lines E, I and Jc is 5-94 per cent. 

A photographic determination, by Surugue (1937), of the relative intensities 
of internal conversion corresponding to the three components of the I line gave 
L,=1, Ly =0-06, Ly,=0-004. If similar ratios hold in the other two 
cases we have L;=5-58%, L,,=0-33%, Ly,;=0-02% for the total intensities of 
primary L ionization due to the three lines concerned. 

Secondary L ionization will result from the transference of most of the K 
ionization, either by Auger transitions, or x-ray emission, to the L shells. 

From Flammersfeld’s figures the total intensity of the K-Auger lines is 
0-89 per cent. This leaves an intensity of 32-23 per cent to be transferred to the 
L and higher shells by x-ray emission. 

From Massey and Burhop’s (1936) tables 53°% of the 32-23 goes to the Ly; 
shell on the emission of the Ka, x-ray, and another 26% goes to the L,, shell by 
the Ka, emission. ‘The rest of the K ionization goes to the M and N shells and 
does not concern us here (L,— K is a forbidden transition). 

Hence the ionization of the three L shells resulting from the primary K 
ionization that gives rise to x-rays is, per 100 disintegrations, L;=0, L,,=8:38, 
Gay =17-l. 

The small contribution of the K-Auger lines to L ionization can be estimated 
from the work of Ellis (1933). From the relative intensities of the K-Auger 
lines the L ionizations are obtained as follows: L,;=0:59%, Ly=0:38%, 
Li =9°55%.- 

Taking the three calculated contributions together, the total intensities of 
iveton ation - are, therefore,» L,—6-17, 1,=9-09, Li, =17:67 per 100 
disintegrations. 

Ze Ol; 

The primary ionization in element 81 is due to the two y-rays A and Ga, 
which have energies of 39-8 and 286-9 kev. respectively. 

The intensity of the K internal conversion line of Ga is 0:11 per cent and the 
L internal conversion lines are so weak that Flammersfeld did not detect them. 
In the following calculations therefore the effects of this y-ray will be neglected. 

y,, having an energy of only 39-8 kev., produces no K internal conversion line, 
but is strongly converted in the L shells. Flammersfeld was unable to resolve 
the lines due to the three shells, but determined their total intensity as 20-9 per 
100 disintegrations. Included in this figure, it must be remembered, is the 
intensity of 3-5 per cent due to the A’ line mentioned previously. ‘Thus an intensity 
of ionization of 17-4 per cent remains to be divided between the three shells. 

Surugue’s (1937) relative intensities for the three L-lines are L; =1, L,,;=0-1, 
Lyt=0-01. Since he did not resolve the A’ line these values are somewhat in 
doubt. For our purposes, however, they are probably accurate enough. We 
therefore divide the 17-4 per cent according to the above ratios and obtain for 
the primary ionizations: L,y=15-67%, Ly=1:57%, Ly =0°16%. 

The relative intensities of L ionization for elements 83 and 81, calculated above, 
are further considerably modified by two very intense Auger transitions that were 

* The A’ line is almost coincidental with the A line, and its intensity is calculated from a knowledge 


of the intensity of the E line (the L. conversion of the same y-ray) (Kinsey 1948). 
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first identified by Coster and Kronig (1935). Coster and Kronig were concerned 
in explaining the anomalous intensity distribution of the x-ray emission lines and 
their satellites, and explained the anomalies on the basis of the high probabilities 
of the two transitions L;>Ly;Myy and Ly>LyMy. The Coster-Kronig 
phenomenon was first taken account of in relation to the Auger effect in radio- 
active decay by Kinsey (1948). 
Obviously Coster-Kronig transitions can take place only when 
EL >F +E 


lat: he ey ey, 
It is interesting to note that energetically this condition holds only for elements 
of atomic number less than 53 and greater than 74. For Z=81, 83 this condition 
most certainly obtains. 

Coster and Kronig estimated that as much as 80°, of the initial L; ionization 
is transferred to the Lj, shell by the transitions at present in question. Although 
the precise figure is rather uncertain, we take this value to obtain ultimate 


intensities of ionization per 100 disintegrations in the three L shells as follows. 


Table 2 
a Li Ly Lor 
83 12S 9-09 22°61 
81 3°13 1:57 12:69 


It may be remarked here that the Coster-Kronig transitions should give rise 
to very strong Auger lines in the region under 1 kev. So far, due to the difficulty 
of accurate neutralization of the earth’s magnetic field for these low energy 
electrons, this energy region of the spectrum has not been investigated. 

Using the ionization figures of Table 2 as a guide, an attempt may now be. 
made to fit the observed L-Auger lines (Table 1) intoascheme. For this purpose, 
values of the L, M and N ionization energies used were taken fron Lindh’s (1927) 
table. 

As indicated by Ellis (1933), the energy of the Auger lines cannot be obtained 
accurately from such x-ray data, for it must be remembered that the atom left by 
an Auger electron is ionized doubly, not singly as in the case of the x-ray ionization, 
and no systematic information regarding x-ray levels in deeply ionized atoms is 
available. For the K-Auger lines of 'ThB Ellis estimated the additional binding 
energy for the second L ionization between 100 and 300ev. In the present case, 
in some disintegrations, the second ionization occurs either in an M or an N level, 
and consequently the additional binding energy will be considerably less (probably 
well under 100ev.), but in others the atom is already doubly ionized, after a 
Coster—Kronig transition, and becomes triply ionized after the Auger transition 
which is observed. In these disintegrations obviously a slightly greater additional 
binding energy will apply. 

From Table 2 we naturally attempt to explain the observed lines as originating 
in transitions from.the most frequently ionized states, namely 83 L,,,, 81 Lin 
and 83L,,. Table 3 gives possible transitions on this basis. 

It is satisfactory that on this basis we find, for any strong line originating from 
the Ly; shell of element 83, a corresponding line for element 81. Such pairs 
are ab, cd and ef. Since the strongest x-ray originating from the L,,; shell is 
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the LyyMy(La,), it is further satisfactory that the most prominent Auger lines 
involve the same transition. 

Lines g and h, on energetic grounds, are not likely to originate from the Lj, 
shell of either element. They probably originate from the L,, shell of element 


Table 3 
Line (kev.) Element 83 Energy (kev.) Element 81 Energy (kev.) 
h 12-62 LiMiyNyr 12:58 
ey 12-26 Lo MirNiy: ib Pori| 
f 10-39 La ivNy. 10-36 
e 9-82 LiMyNy 9-86 
d ~8-10+0-05  LyMyNiy 8-12 
Cc 7°55+0:05 LiMyNiv 7:74 
b 7:10 LinMyMrr 7-11 
a or7l LinMyMi 6°81 


83. The main x-ray originating from this shell is the L,,M,y(Lf,), and they 
have therefore been attributed to the internal conversion of this x-ray. It is 
satisfactory that the relative intensities roughly deducible from Figure 4 give 
general support to these identifications on the basis of the ionizations given in 


Mable 2. 
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Photoelectric Dissociation of the Deuteron 


By C. H. COLLIE, H. HALBAN anp R. WILSON 
Clarendon Laboratory, Oxford 


MS. received 5th December 1949 


ABSTRACT. The cross section for the photo-disintegration of the deuteron has been 
measured by detecting the photo-protons formed in a high pressure ionization chamber. 
The values obtained are 13-91-+1 x 10-28 cm? for 2:62 Mev. y-rays from radio-thorium and 
15-6+1 10-28 cm? for the 2°76 Mev. y-rays from radio-sodium. 


§1. INTRODUCTION 


OR the investigation of nuclear forces two-body problems are of primary 
interest. This is as much due to the possibility of solving the wave 
equation exactly as to the relative simplicity of the experiments 

involved. 

It is well known that the assumption of a strong nuclear attraction (about 
20 Mev.) acting over a short range (2:8 x 10-"%cm.) is adequate to reproduce 
the main properties of the deuteron and the interaction of slow neutrons and 
protons. The direct production of 7-mesons by the interaction of high energy 
nuclear particles leaves little doubt that nuclear forces are due to a meson field 
whose nature is to be determined by a combined appeal to direct experiment 
and the restricted principle of relativity. The latter is of course an appeal to 
imaginary experiments of which oneis not supposed to doubt the issue. As is well 
known (Bethe 1940), the latter principle severely restricts the interaction potentials 
which come into consideration to terms of the type e-“’/r, in which the range of 
the forces is unambiguously given in terms of the 7-meson mass p by « =ye/h. 
Thus on a strict meson theory the range of the nuclear forces is not a disposable 
constant unless one imagines that mesons other than z-mesons play a part. 

Once the nature of the interaction between the two nucleons has been 
determined, the wave functions of the deuteron can be obtained by solution, 
usually numerical, of the wave equation. A knowledge of the wave equation 
enables the interaction of the deuteron with an electromagnetic wave to be 
calculated by conventional methods and thus a direct comparison to be 
made between theory and experiment. ‘The experimental quantities which 
can be measured are : (a) the absolute value of the cross section at a given value 
of incident quantum energy ; (0b) the relative variation of cross section with 
energy ; (c) the angular distribution of the emitted particles, which is in effect 
a measure of the relative proportions of the disintegration due to the electric 
and magnetic vectors of the incident wave. 


The present measurements are concerned with the first two of these 
possibilities. 
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§2 PREVIOUS WORK 

Although many measurements of these quantities have been made 
previously, owing to the technical difficulties of the measurement and the 
exigencies of the times they have lacked the precision necessary to distinguish 
between the competing theories which have been put forward, and have done 
little more than lend support to the general wave mechanical treatment of the 
deuteron. ‘The most important measurements which have been made in this 
field are tabulated in Table 1. The more precise of them had not been made 
when our work was started. 


able i — = 


Cross Section 


y-rayenergy oX1028cm? 


Author Date Method 
(Mev.) 

Chadwick and Goldhaber 1935 —_ Tonization chamber 2-620 6+6 
Halban 1938 Neutron standard 2-620 10+1 
Richardson and Emo 1938 Cloud chamber 2:758 11+11 
van Allen and Smith 1941 Ionization chamber 6:13 11:6+1:5 
Lassen 1948 Proportional counter MUS) 25+8 
Russell 1948 Neutron standard 2:758 1644 
Snell, Barker and Sternberg 1949 Neutron standard { se poe 5; : 
Wilson, Collie and Halban aeeae Ionization chamber { oe i are 
Waffler and Younis 1949 Photographic plate 17-6 8+3 
Barnes, Stafford and Wilkinson 1950 Ionization chamber af ee see 

Z ah 17 6 8 5 Sf 1 2 


Angular Distribution 
H=counts 


be Sei at 90°/counts 
‘ at 0° 
Chadwick, Feather and Bretscher 1938 Cloud chamber 2-620 0-25 
Halban 1938 Activation detector 2620 0:07+0:06 
Myers and van Atta 1942 BF, counter ee £0 \ 0-87 -+-0:07 
Graham and Halban 1945 BF, counter 2:620 0-20+0-05 


Jf Proportional counter AS 2-758 0:15-++0:03 


Lassen 1948 

¥ for protons 
Woodward and Halpern 1949 ain Aa ee ie 
Hamermesh and Wattenberg 1949 BF, counter 2°758 0-19 
Gibson, Green and Livesey 1947 Photographic plate 6:13 0-2 
Meiners, Smith and Slack 1949 BF, counter 2°758 0:16-+0-02 


§3 METHOD OF MEASUREMENT 
The measurement of the absolute cross section of the photo-distintegration 


of deuterium, 
7H + hy = {H+ jn+Q, 


is in principle simple, all that is required being a known source of y-rays, a 
specimen of a deuterium compound whose deuterium content is known, and a 
detector capable of detecting either the proton or the neutron. In practice, 
however, it is found that there is some difficulty in evolving an experimental 
technique capable of the overall accuracy to within 5°% which was aimed at as 
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being the least which could be of use in a comparison between theory and 
experiment. After careful consideration it was decided to detect the proton 
rather than the neutron emitted in the disintegration. It is true that the great 
penetrating power of neutrons, together with the ease with which they can be 
detected, enables the difficulties due to the small value of the photo cross section 
(~10-27 cm?) to be in part off-set by the use of a large amount of deuterium in 
the target. However, none of the devices (e.g. BF; chamber) for detecting 
neutrons are absolute and the lack of a reliable neutron standard* as yet prevents 
satisfactory absolute measurements being made in this way. This objection 
does not apply with equal force to measurements of the variation in the relative 
cross section with incident quantum energy. A further disadvantage of methods 
of measurement based upon the detection of the neutron, when using the very 
convenient y-rays from ThC, lies in the relatively large neutron background 
of this source. The unwanted neutrons are produced by the intense «-particle 
bombardment of the light elements unavoidably included in the carrier used 
in the chemical separation of the radio-thorium from the mesothorium parent. 


§4. THE IONIZATION CHAMBER 
The protons were detected in a spherical ionization chamber, filled with 
deuterium under pressure, connected to a linear amplifier. A cylindrical 
chamber was used in some preliminary experiments. The design of both 
chamber and amplifier are dictated by the experimental difficulties which have 
to be overcome, and are as near the optimum as was possible by simple means. 
Details of the ionization chambers are shown in Figure 1. 


il || 


(a) (0) 
Figure 1. Deuterium ionization chambers. 


(a) Shows the cylindrical chamber used in early measurements. 
(b) Shows the spherical chamber used in later measurements. 
A=collecting electrode of nickel. 
B=special (low capacity) seal. 
C =cylindrical copper chamber 4 cm. X 4 cm., volume=50 cm. 
D=Kovar-to-glass seal type Cathodeon, GMSE. 363. 
E =outer copper pressure case and guard ring. 
F =final soft solder joint. 
I =spherical copper chamber of volume 30 cm3. 


The metal-glass seals are soldered with silver tin solder softening at 300° c. 


Using the y-rays of ThC (2-62Mev.), the energy of the photo-proton is 
210 ky., so that its range will be 1-0 cm. in deuteriumat N.T.P., and it will produce 
about 6,000 ion pairs. This small ionization, even under the most favourable 
circumstances, will not produce a change of potential of the first grid of the 


* The position in this respect has recently improved, but no standard reliable to better than 5% 


enone, and the neutron flux is spread over a very large energy range which makes comparison 
ifficult. 


‘ 
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amplifier greater than 70 microvolts. The pulses due to the protons can be 
distinguished from random noise and other pulses by their uniform size, and 
advantage must be taken of this fact to obtain sufficient discrimination for 
absolute measurements in the presence of recoil protons from unavoidable 
neutrons and random fluctuations in the ionization current produced by the 
y-rays traversing the chamber. 

The pulses from the photo-proton are not quite uniform in size owing to 
the momentum of the incident quantum which must appear in the combined 
momentum of the proton and neutron. ‘This effect amounts to a spread of 
+ 26 key. in the energy of the recoil protons. Recombination and shielding of 
the inner electrode by the positive ions will also cause a spread in the measured 
pulse height. The detection and measurement of this spread affords a 
satisfactory indication of the correct functioning of the whole system. 

Some of the photo-protons will arrive at the walls before they have expended 
all their energy in producing ions in the gas, and this ‘ wall effect’ will be an 
additional cause of variation in the recorded ionization produced by the 
photo-disintegration. The magnitude of this effect, which can be fairly simply 
calculated provided the range-energy relationship is known, will become 
proportionally less as the gas pressure in the ionization chamber is increased, 
and will be least and most easily evaluated for a spherical chamber. Another 
advantage of a spherical chamber is that it has no corners in which the field 
driving the ions to the central electrode is very weak, thus increasing the time 
available for recombination. ‘These considerations determined the design of 
the chamber : a pressure of about 6 atmospheres in a 4cm. diameter sphere 
reduces the wall effect to manageable proportions, but brings other difficulties 
in its train. 

Owing to the small absolute value of the photo cross section, the ionization 
in the gas due to y-rays greatly exceeds that due to the protons and increases 
with the pressure. ‘The statistical fluctuations in this constant ionization 
current ultimately approach, and even exceed, the pulse of ionization due to the 
proton, and this sets an upper limit to the gas pressure which can usefully be 
used. If the collection time of the ions is a few milliseconds, as is the case for 
molecular ions, this ‘ piling up’ of the y-ray pulses in the time required for complete 
collection makes measurements at 6 atmospheres impossible, but if the negative 
ions remain as free electrons, so that the collection time is a fraction of a 
microsecond, the statistical fluctuation in this short time is reduced sufficiently 
for one to be able to choose a deuterium pressure (about 6 atmospheres) at 
which the wall effect is small and the y-ray noise sufficiently low to make 
measurements possible. 

The achievement of electron conduction can only be attained if the gas is 
extremely pure and free from molecules such as oxygen, which have large 
electron attachment coefficients. It was widely believed that even in pure 
hydrogen electron conduction ceased at about 3 atmospheres pressure but, 
as was first shown by Stafford (1948), this is not so, electron conduction being 
maintained up to 80 atmospheres in hydrogen of extreme purity admitted 
through a palladium leak. For our work, using pressures below 30 atmospheres, 
we found that a less rigorous but more convenient system of purification was 
sufficient. After careful cleaning and final assembly, the chamber was pressure- 
tested and usually x-rayed to see that the central electrode was symmetrically 
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placed. This avoided wasting time and deuterium on a chamber which would 
ultimately prove useless. It also provided a record of the exact position of the 
inner chamber with respect to the outer wall. The chamber was then connected 
to a vacuum system through an ordinary union rendered vacuum tight with 
hard Cavendish wax (Everett’s No. 3) and baked out at as high a temperature 
as the last soldered joint allowed. The baking was usually for two days at 160° c. 
When no more occluded gas was given off the chamber was flushed out with 
deuterium and then filled under pressure from a cylinder of gas which had been 
purified by distillation from charcoal cooled with liquid oxygen. The whole 
of these operations were greatly facilitated by the use of the small diaphragm 
taps described by two of us (Halban and Wilson 1949) which have the merit of 
acting both as vacuum and pressure taps. The pressure was indicated by a 
miniature Bourdon gauge, but the actual pressure in the ionization chamber 
was determined at the end of a series of measurements by expanding the gas 
into a known volume and measuring the pressure on a wide-bore mercury 
manometer. Specimens of the gas were retained for analysis. On occasions 
chambers which did not show electron conduction immediately after filling 
developed it in the course of a few days, showing that the clean copper walls 
of the chamber were capable of absorbing the impurity, presumably oxygen, 
responsible for the unwanted electron attachment. 


$5. THE AMPUIPIER 

In order to make full use of the small collection time (found experimentally 
to be about 0-2 psec.) so desirable from the point of view of reducing y-ray 
noise from the ionization chamber, the bandwidth of the amplifier must be 
increased, thus introducing more noise from the amplifier. ‘The necessary 
compromises were on the whole found experimentally by measuring the 
signal/noise ratio with the pulse amplitude (Figure 3, see below) discriminator. 
The 2-usec. pulse duration found in this way agreed very well with the value 
obtained by a careful analysis of the circuit carried out subsequently (Wilson1950). 

The first, or pre-amplifier, stage avoided input circuit noise by using a high 
value (about 1,000 megohms) of the grid leak so that the current pulse from the 
photo-proton produces a rapid quasi-permanent increase in the grid potential. 
Since it is important to avoid the scattering of y-rays by apparatus near the 
ionization chamber the first valve was connected to the main amplifier through 
a transmission line which was matched into a grounded grid triode. The 
resulting low gain of the first valve, which would in a conventional circuit allow 
appreciable noise from the second valve, is off-set by the degenerative effect 
of the high impedance of the first valve in the cathode circuit of the second. 
This circuit was adapted from circuits used in radar by Wallman, MacNee, 
and Gadsden (1948). The output from this pre-amplifier stage consists of a 
sharp voltage rise due to the arrival of the electrons at the collecting electrode, 
and before passing through a conventional amplifier this was converted into a 
true voltage pulse by means of the network shown in Figure 2. The advantage 
of an independent pulse-shaping network is the very much greater control one 
has over the pulse shape. When working near the noise limit the relatively 
small improvement in signal/noise ratio obtained by using a delay line as first 
suggested by Frisch results in a substantial gain in discrimination against 
unwanted signal. Working in this way, a factor of 1:5 in the signal/noise ratio 


’ 
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over that given by a simple differentiating circuit can be obtained. Exact analysis 
by Wilson (1950) shows that the best shape of pulse occurring at a mean time 
t, Should be exp {—|(¢—#,,)/7|}, in which the parameter 7 is determined from 
the ratio of the valve noise to the y-ray noise. In this way one could obtain a 
further 20° increase in signal/noise ratio, but so far no sufficiently stable circuit 
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Figure 2. Pulse-shaping circuit. 
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Figure 3(a). Single channel kicksorter. 
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Figure 3 (6). Discriminator and pulse shaper. 


for producing this pulse shape has yet been devised. After the pulse is shaped 
it is further amplified by a wide-band amplifier of conventional design and 
applied to the discriminator unit shown in Figure 3(q), whose purpose is to 
enable only counts above a certain threshold, or between two thresholds, to be 
recorded. 

The whole of the power supplies were stabilized against fluctuations and 
the overall gain was constant to better than 1%. 
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§6. THE y-RAY SOURCE 


Since the photo cross section varies rapidly with quantum energy, the 
monochromatic sources from excited nuclei are the most useful. ‘Three ways 
of determining the number of y-rays emitted from such a source are available, 
namely : (a) Direct ionization measurements in which the theoretical cross 
sections for Compton scattering and pair production are assumed (Gray 1936). 
(b) Absolute calibration based on the coincidence method first proposed in 
this connection by Rotblat (private communication) and Dunworth (1940). 
Extensive measurements on the y-rays of radio-sodium, "4Na, have been made 
by Putman and Willment (1950). To apply this method the decay scheme of 
the nucleus must be well established with a B-ray spectrometer (Siegbahn 1946). 


(c) Measurement. of an accompanying «disintegration, which can be made — 


much more easily. This method is applicable to the y-rays from a radio-thorium 
source. 


The number of radioactive substances of reasonable half-life emitting y-rays ~ 


of greater quantum energy than the binding energy of deuterium (2-220 Mev.) 
is limited, only RdTh (2-62), 24Na (2-758), Ga (2-51), #°La (2-56), *®Mn (2-7) 
and Ra (2-42) being readily available. Only the first two of these have been 
used in the work described here. 


§7. THE RADIO-THORIUM SOURCE 


The radio-thorium source used consisted of about 400 Mc. of radio-thorium, 
carefully purified and mixed with platinum black before sealing in a stainless 
steel cylinder. The preparation of this source was undertaken by Thorium 
Limited. It was hoped that the use of platinum black might reduce the number 
of fast neutrons emitted by the source since these were thought to come from 
the «-particle bombardment of light elements in the thorium oxide carrier used 
in the separation of the radio-thorium from its parent mesothorium. This 
device proved unsuccessful and the source emitted approximately 2,000 neutrons 
per second. 

Radio-thorium disintegrates according to the scheme 


PS 


The 2-62 Mev. y-rays are emitted from the ThD nucleus, which is left in an excited 
state after the ThC” disintegration. Higher energy y-rays (3-2 Mev.) were looked 
for using the high-pressure ionization chamber but were not found. If they 
exist, their intensity, less than one per thousand disintegrations, is low enough 
to be neglected in these measurements (Bishop, Wilson and Halban 1950). 
‘Thus the absence of any high-energy (5 Mev.) B-rays in the ThC” B-ray spectrum 
shows that every ThC” disintegration is accompanied by one y-ray, so that the 
absolute number of y-rays can be obtained by counting the number of «-particles 
emitted by a known amount of ThC. 

__ This measurement has been made by Ricoux (1937), who measured the total 
lonization current produced by the «-particles from a known amount of ThC. 
To recover Ricoux’s calibration a similar ionization chamber (the well-known 


poe : 
RdTh-> ThX-—> ThEm-> ThA-> ThB> ThC ‘Sm. 
The’ 2 


4 
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Paris pattern first described by Madame Curie (1912) for radium measurements) 
was constructed and calibrated with a known radium source (249-5 mc., N.P.L. 
Certificate No. RG.467). Ricoux’s measurement is not based upon direct 
counting of the «-particles, and therefore depends upon a knowledge of the 
relative ionizations of «-particles of different ranges. He assumed that the 
energy required to form an ion pair was the same for all «-particles. This 
assumption is not quite correct, as was shown by Gray (1944), but a recalculation 
using the new data leaves Ricoux’s result unchanged because the weighted 
mean range of the two «-particles from ThC is the same as that of the «-particles 
from RaC. The radio-thorium source was also measured with a standard 
N.P.L. chamber. Each intercomparison could be made to about 1%, but the 
calculation of the absolute y-ray intensity introduces two further uncertainties. 
Since Ricoux measured the total «-particle ionization of ThC the branching 
ratio eee) at enters directly into the calculation. This ratio is usually 
quoted as 35°, based upon the work of Albrecht (1919), but Kovarik and Adams 
(1938) obtained 33-7°%, and we have adopted 34°%% as a weighted mean. This 
agrees well with the results of Meitner and Freitag (1926), who used a cloud 
chamber. 

Since the whole method of measurement is based upon comparison with 
a radium standard, the absolute disintegration rate of radium must be known. 
This is often taken on Rutherford’s authority (Rutherford, Chadwick and 
Ellis 1932) to be 3-70 x 10!°/sec/curie, but the recent measurements of Kohman, 
Ames and Sedlet (1947) suggest that 3-61 x 101°/sec/curie is nearer the truth. 
Adopting the latter figure, we find the effective disintegration rate of our source 
(No. RdTh. 139) to be 4-54 x 109 quanta/sec. calculated back to 9th June 1948. 
The greatest uncertainty in this figure is the branching ratio in ThC, which 
could easily be in error by 3%. The calculations, however, do include a change 
in the accepted wall and salt absorption correction (Kaye, Aston and Perry 1934) 
to take into account the different geometries of the N.P.L. and Paris ionization 
chambers. 

The disintegration rate of the source may also be based upon the data of 
Shenstone and Schlundt (1922), who compared the number of «-particles from 
a ThC’-ThD source with those from a RaC’-RaD source of equal y-ray intensity. 
The result is 4% lower than that based on Ricoux’s measurements, a figure 
which also agrees with an estimate based on ionization chamber measurements 
made by Gray. On the whole we have given greater weight to the very careful 
and consistent measurements made by Ricoux. It is unlikely that the value 
we have adopted is more than 5% in error. A further comparison of radium 
and radio-thorium sources on the Paris ionization chamber has been made 
by Winand (1939). In these experiments the total «-particle energy was 
measured with a micro-calorimeter and the results agree with Ricoux’s 
measurements to within 1%. 


§8. THE RADIO-SODIUM SOURCE 
The radio-sodium source was prepared by exposing 30 grammes of pure 
sodium carbonate in the A.E.R.E., Harwell, Pile in an aluminium container. 
After shaking well (to minimize small effects due to non-uniform activation 
of the sodium carbonate), the main bulk of the material was transferred to a 
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flat tinned steel box (4 mm. wall thickness) and weighed. An aliquot was reserved 
and its absolute disintegration rate determined at A.E.R.E., using the coincidence 
method. This method of calibration depends upon the validity of the decay 
scheme for radio-sodium given by Siegbahn (1946) and Wiedenbeck (1947). 
Two y-rays of quantum energy 2-76 Mev. and 1-38 Mev. are emitted in cascade. 
The f-ray spectrum is a simple one with maximum energy 1-38 Mev. At least 
99°/, of the disintegrations follow this simple scheme. A further point is that 
the calculation of the intensity from the observed B~y coincident rate involves 
the assumption that there is no angular correlation between the B- and y-radiation. 
This point has been checked by Grace et al. (1949), who showed experimentally 
that the departures from a random distribution are less than 1% of the observed 
B-y coincidences. This method of calibration has been checked by the absolute 
ionization method (Gray 1949). The two methods of calibration agree within 
the experimental error. | 
The half-life of radio-sodium enters directly into this method of calibration ; __ 


the value 14-90 + 0-02 hours found by Wilson and Bishop (1949) in an investiga- 


tion undertaken for this reason was used in reducing the experimental values. 


The radio-sodium sources used emitted some 10° quanta/second, which enabled — | 


the source to be placed about 20cm. from the counter. 


Experimental Results 


Typical differential bias curves using radio-thorium and radio-sodium are 
shown in Figures 4 and 5. The hump due to the photo-protons is clearly defined 
and measurable, but it is clear that it is only just separated from the y-ray 
background when the lower energy protons from radio-thorium are being 
counted. 

In order to determine the true number of photodisintegrations produced 
in this case one has to decide what the differential bias curve would have been 
had no photo-protons been produced. This was done in two ways. A radium 
source was brought up to the counter and its position adjusted until it gave the 
same y-ionization current in the counter as had been given by the ThC source. 
The differential bias curve obtained under these conditions, which represented 
the noise due to the ionization current from the radium y-rays, was assumed 
to be the same as the contribution of the ThC y-rays to the composite differential 
bias curve in which the noise due to y-rays and the tail of the proton peak are 
superimposed. ‘The true photo-proton count in the region of overlap was 
obtained by subtraction. The assumption that the shape of the differential 
bias curve for y-ray noise depends primarily on the ionization current and not 
on the y-ray energy is reasonable since most of the electrons crossing the chamber 
are Compton electrons produced in the chamber wall. Alternatively the profile 
of the well-resolved proton peak obtained with radio-sodium y-rays was 
carefully measured and corrected for wall effect. This measured photo-proton 
peak from 2-62 Mey. y-rays was then used as a basis for extrapolating the tail 
of the measured proton peak from ThC y-rays into the region where it overlapped 
with the y-ray noise. 

The agreement between the two methods was satisfactory and the correction 
itself was often zero and always less than 2°/, of the total counts. The total 
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Figure 4. Differential bias curve for photo-protons from radio-sodium y-rays. 
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observed number of disintegrations obtained in this way still has to be corrected, 
as follows : 


(a) Wall effect. 

Some 8° of the protons liberated in the chamber are not counted because 
they hit the walls. The number was estimated in the following way. ‘The 
number of protons of range R produced within range of the wall of the chamber 
of area A is AR, of which }AR will strike the wall if it is plane. This remains 
true to a second approximation for a spherical wall, and the main difficulty in 
calculating the wall effect is due to other causes. 

Some protons will strike the wall towards the end of their range and will 
produce sufficient ionization to be counted. Thus, to give a concrete example, 
if the counting level is 150 kev. and the mean pulse height is 290kev., the wall 
effect is approximately 1A(Rygy— Ryo), Where Ry and Ry) are the ranges of 
290 and 140kev. protons respectively. The ranges used are given below and 
were calculated from the ranges of protons in air and the stopping power of 
hydrogen given by Gray (1944). 


Estimated Proton Ranges in Deuterium at 760 mm. and 15° c. 


Residual energy (kev.) 50 100 150 200 250 300 
Range (mm.) 1-92 3°91 6°33 O21 12-41 16-15 


To some extent the correctness of the estimated wall effect can be checked 
by observing the minimum of the differential bias curve for radio-sodium 
(Figure 4) in the region where it should have fallen to zero. If the whole of 
the residual count is attributed to wall effect, the agreement between calculated 
end observed wall effect is good. A similar check could not be made when 
radio-thorium y-rays were used owing to the y-ray background overlapping 
the tail of the proton curve. 


(b) Absorption of the y-rays in the walls of the counter. 

The absorption of the y-rays in the walls of the counter is appreciable and 
was estimated by extending the calculations of Evans and Evans (1948). The 
calculations are to a certain extent simplified because quanta scattered through 
more than 10° by Compton collision can no longer disintegrate deuterium. 
Thus scattering from the ends of the ionization chamber can be neglected. 
For the various electrode systems used the mean effective thickness @ of the 
walls was related to the measured thickness ¢ by the relation #@=1-77¢ for the 
cylindrical counter and #@=1-57t for the spherical counter. From these values 
the following correction factors to be applied to the crude value of o were 
calculated : cylindrical counter, radio-sodium 1-:0596, radio-thorium 1-0605 : 
spherical counter, radio-sodium 1-053, radio-thorium 1-0538. , 

The absorption was also investigated experimentally by surrounding the 
counter with a copper cylinder of equal wall thickness. Agreement within 
the experimental error was obtained. 

A further correction due to Compton scattering in the source must be made. 
. Quanta scattered in the source will contribute to the observed y-ray 
intensity as observed with an ionization chamber, but nearly all of them will 
have a quantum energy below the deuterium threshold because of the energy 


lost to the recoil electron. This effect was estimated roughly at 1:45°% for 


; 
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measurements on the Paris chamber and.1-9°% for measurements on the N.P.L. 
chamber. ‘This particular type of correction is not necessary in the case of the 
radio-sodium source, which was calibrated absolutely through an aliquot: an 
ordinary correction for absorption in the source must be applied, but was in fact 
small (1-8°%) owing to the shallow container in which it was placed. Because 
of the low specific intensity of the radio-sodium source the container was rather 
large and its average distance 7 from the ionization chamber had to be estimated 


from the formula 
1 1 a x+y 
Amato ap 


in which p is the corrected measured distance to the centre of the counter, a is 
the counter electrode radius, and 2 and 2y are the breadth and width of the 
Square source. 


(c) Purity of the deuterium. 


The deutertum was made by electrolysing 99-5°% heavy water, so that, 
allowing for an enrichment factor of 6, the gas produced electrolytically could 
have had up to 3% of ordinary hydrogen in it. After the plant had been running 
for some time the composition of the gas should have slowly approached that 
of the water. Since the possibility of exchange with ordinary hydrogen occurring 
in tap grease, etc., could hardly be avoided, samples of the deuterium gas from 
the ionization chambers were analysed with a mass spectrograph by Messrs. 
K. Wood and J. Blears at the Metropolitan-Vickers Research Laboratory, 
‘Trafford Park, Manchester. ‘The deuterium concentration did not fall below 
98-7°%, in any of the samples analysed. ‘The measurements were made by 
assuming that the measured icnization current in the spectrometer was 
proportional to the concentration of the corresponding nucleus in the gas. 
The peaks measured were due to D,+ and DHt+ under conditions in which it 
had been verified that the discrimination between mass number 4 and mass 
number 2 was less than 10°%. 


(d) Neutron background. 


The radio-thorium source emitted an appreciable number of fast neutrons 
which produced recoil protons in the ionization chamber. By bringing the 
source up to the side of the counter a countable number of high energy (up to 
500 kv.) recoil protons was observed having a continuous energy distribution. 
The number to be expected in the energy range of the photo-protons could not 
be observed directly under these conditions because of the y-ray fluctuations, 
but was estimated to be 2:2°/, by extrapolating the observed distribution into the 
y-ray background. The extrapolation could be carried out with some confidence 
since it was based upon a blank experiment using a counter filled with ordinary 
hydrogen. No recoil protons were observed when a hydrogen-filled chamber 
was exposed to the radio-sodium source. 


Results 


Complete data for a typical run are shown in Table 2, and the collected 


results in Table 3. 
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Table 2. Radio-Sodium y-Rays 


Counter B, Run 2, 25th October 1948. Clock time running 0-3% fast. 
Observed counts are already corrected for natural background of the ionization chamber. 


Corrected to 0000 G.M.T. 


25.10.48. 

(1) (2) (3) (4) (5) (6) (7) 
30 16:12 196-7 1-2 416-9 25 27 
20:26 165-1 0-9 425-9 DS 205 

22-14 139-5 0-4 396°5 el DY 

22:59 15327 0-75 446-2 22 29; 

2328 133-0 0:77 Bo 7c D233 297, 

62 16:58 188-3 0:8 413-5 ile7/ 4% 
16-21 188-7 1&3 412-8 2°8 3% 

34 20-11 161-7 0:9 412-3 2°3 295 
Mean counting rate at bias of 30 is 417-0 19 


(1) Discriminator bias; (2) mean time; (3) integral counting rate per minute ; 
(4) differential counting rate per minute; (5) integral counting rate; (6) differential 
counting rate; (7) statistical accuracy in integral count. 


The above table shows the experimental data from which was constructed the 
approximately horizontal portion of the differential bias curve shown in Figure 4. 
The remainder of the curve showing the pronounced peak at a bias of 47 is obtained in 
a separate experiment and is only used for calculating the subsequent corrections. 

The total number of photo-protons, 417 per minute, is given by the integral counting 
rate at a bias 30, i.e. all those particles whose amplitude exceeds the threshold bias of 30. 
The choice of 30 is arbitrary, any point on the flat part of the differential bias curve could 
have been chosen. If a lower voltage had been chosen, a higher integral count would have 
been obtained, but the correction added for wall effect would have been less; if a higher 
bias voltage had been chosen, the integral count would have been smaller but the wall- 
effect correction to be added would have been greater. 

If C is the mean counting rate at a bias of 30 (=417/minute), N the source strength 
at zero time=1-:632 x10" disintegrations/minute, p the distance between the centre of 
the sodium source and the centre of the counter=21-:91 cm., 7 the distance when corrected 
for the finite size of the source.and counter=22:03 cm., p the counter pressure =6°865: 
atmospheres at 0° c. or 7:23 atmospheres at 15° c. (for use in calculating the wall-effect 
correction), V the volume of the counter = 30-14 cm*, L Loschmidt’s number =2:687 x 10, 
then the uncorrected value of the cross section o is given by 


4nCr? 
o=>= 
2LVpN 


The further corrections which have to be applied are given in Table 3, which also gives. 
the results of other runs. 


=13-94 x 10-28 cm?. 


Discussion 
Theoretical values of the cross section have been given by many authors, 
but those of Hansson and Hulthén (1949) correspond well with the conditions 
of our experiment. These authors have based their calculations upon the 


Moller-Rosenfeld theory and its variants, in which the static interaction in 
various states is of the form given below. 


State 
Type of theory aS gS Be 
Moller-Rosenfeld (MR) —A —B B35) ee 
Symmetrical (O) —A —B Ome 
Neutral (S) aay) | 7 =P y 
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Table 3. Consolidated Results 


o Wall Source Wall 
(uncor- absorp- _absorp- effect D,+H,/D, o Cine Gores 
rected) tion tion (%) 
A, Na 13-60 1:0596 1-015 TTP? 1:00835 15-98 
RdTh 12-01 1-0605 1-015 4-7 13-68 1-17, 
B, Na 13-07 1-053 1-018 5:88 1:0082 15-01 
RdTh 12E3S 1-0538 1-013 4-1 13-84 1-085 
A, Na 12-97 1:0596 1-018 6:58 1:0161 15222 
RdTh 1219 1-0605 1-015 43 14-27 1-066 
Be Na 13-94 1-053 1-018 6°87 1-0080 16-17 1:156 
RdTh 12-43 1-0538 1-015 4-2 13-99 
A,;* Na 16-33 1-053 1-018 5:88 1:0136 ——leak—— 
RdTh 12-98 1-0538 1-015 5-78 
B,* Na 14-39 1-053 1-018 5:67 1-0087 15°59 OU’, 
RdTh 12-33 1:0538 LE OS 4-18 13-94 
Na 15-6 
Meand Rip 13-9 ios) 


* The counters had both begun to leak after these fillings. With filling B, two radio- 
thorium counts were taken which showed that the leak was small. The pressure was 
measured immediately after one of the counts so that there is less than $°% error introduced 
in the cross section quoted. 

The estimated accuracy of these measurements is about 6% for the radio-thorium 
measurements and 5°% for those using radio-sodium y-rays. Since both measurements 
are made under very similar conditions, the observed ratio 


Gore Oorga— tel Da 0-08 
should be fairly immune from the effect of unsuspected systematic errors. The radio- 
thorium cross section is different from the value (14:8+1)x10-*8cm? given in a 
preliminary note (Wilson 1948). The latter figure is an error owing to the inadvertent 


use of a wrong estimate of the effective number of quanta emitted from the radio-thorium 
source. 


The values of A and B are adjusted to fit the observed neutron—proton scattering 
cross section and the known binding energy of the deuteron. ‘The absolute 
cross section obtained is sensitive both to the binding energy of the deuteron 
and to the meson mass assumed. 

As pointed out by Hansson and Hulthén, a good general fit with our results 
is obtained for a meson mass of 300 electron masses and a high (2-237 Mev.) value 
of the deuteron binding energy in agreement with the work of Bell and Elliot 
(1948). The degree of agreement obtained is shown by the results given below, 
taken from Hansson and Hulthén’s calculations for a 2-237 Mev. binding energy. 


Meson mass = 200 Meson mass = 300 
hv MR O N MR O N 
2:62 133 13-6 14:6 11-8 11-9 12-1 
2°76 16:5 16°9 18°5 14-3 14-4 14-8 
6:20 25:0 25:6 28:0 21-2 PLS) 22-9 


68-2 
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There still remains a substantial discrepancy with the measurements of 
van Alien and Smith, who found o =(11-6 + 1-5) x 10-*§ cm? for 6-2 Mev. quanta, 
and it is in this experimentally difficult energy range that the discrepancy 
between the cross section calculated from a Yukawa potential and the simple 
short range theory of Bethe and Peierls (1935) is greatest. Recent measure- 
ments (Barnes et al. 1950) are, however, in much better agreement with the 
calculated cross sections for a Yukawa potential. A certain discrepancy still 
remains between our results and the theoretical cross sections. This can well be 
expressed in terms of the ratio 09.76/G9.62, which we find to be 1-12 + 0-08: ~The 
value derived from the theory is not very dependent upon meson mass or binding 
energy, and has a value of about 1:22. At this stage it is not possible to say that 
this discrepancy is well established. Owing to the very involved nature of the 
calibration of the radio-thorium source we are inclined to attach greater weight 
to the radio-sodium results, which are in fair agreement with the theoretical 
cross section based upon a high meson mass. 

Further experimental work is being undertaken in collaboration with 
Dr. Siegbahn to establish the relative y-ray intensities of the available sources 
upon a sounder basis. 
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Gn the Lateral Spread of Extensive Air Showers 
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ABSTRACT. ‘The lateral and angular spread of extensive showers is calculated in an 
approximation where ionization loss and incomplete screening are neglected (approximation 
A). The atmosphere is taken to be homogeneous. It is shown that the spread at the 
cascade maximum is equal to the average spread, but the spread increases considerably 
throughout the shower; the above results are in qualitative agreement with results given 
by Borsellino. The quantitative discrepancy between the results of Borsellino and other 
authors is commented on. 


§1 
N a recent paper Borsellino (1949) has calculated the lateral extension of air 
if showers. His results differ from those obtained by other authors, see for 
example Nordheim and Roberg (1948); they also differ from results 


previously given by the author (1948). I wish to take this opportunity to elaborate 
my previous calculations and also to show how the discrepancy with Borsellino’s 


results arises. 
S22) 
Considering a cascade in a homogeneous medium, the mean square angle 
of scatter or the lateral spread can be expressed according to Nordheim as follows : 


(92) =A,? (22) (?) =A? (22) he caer? (1 
where 
5 Wo 2 
Lee ja f q(wo, w’, £—C')q(w’, w, £)C'?% (=) dw’ |q(wo, w, ©), 
Fee Vel Meuron: er, (2) 
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Here q(w, w, ©) dw is the average number of electrons with energies in an interval 
w,w+dw at a depth ¢ produced by a primary of energy wp. 

The integral (2) can be evaluated analytically provided we use the 
‘approximation A’, i.e. neglect ionization loss and incomplete screening. In 
this approximation we have (e.g. Janossy 1948, § 379) 


1 Yo ttc Wo y 
= —) K OV 2) Ee 3) 
ueiy, 20) = = |" (8) KO. dy, (3) 
DAG as - = fs 
K(y, ) = a ees pare exp ( "30 eee (4) 


where @,, @» stands for a,(y), a9(y). Inserting (3) into (2) we find 


2 a5 © day’ Ts tic y! 1 ! / 
Agape, t= | at [SE |" (Et) KO E—e) dy 
@ « j w 


(at) "5. Gal ae 


Yo" .: : cS w ‘ Y 3 Mu" / 12, Mu" 
| (=) KE (5) 
J yp" ic 

Interchanging the order of integrations we can carry out first the integration over ~ 
w’. Collecting the factors containing energies we have 


rwo N=l=3 if W y"—2 1 
Wy we-9" | we day eat (6) 


i eae jay See 
provided Re(y"—y’)>2. We introduce (6) into (5); this is legitimate for 
Vo >No +2. The path of integration into y’ can now be deformed into a loop 


round the singularity at y’ =,” — 2 and can be evaluated by the method of residua. 
This gives 


1 g rYyo tio y 
Aj?q(Wo, &, C) = \ dl | (22) K(y, 6-0) K(y +2, ()C2* dy, 


2T1Wy . Oss 


where we have written y for y”—2. ‘Taking the Laplace transform of (7) with 
respect to ¢ we find 


ne Lee one OL ,(y +2) 
4c A 2 as 0 A Sy 
F € i G(Wo, w, C) dl it ies (22) L(y) a2 CA Reme etey Ve? (8) 
where L,y)= F eK (ysthdl 2) a ee (9) 
and with the help of (4) 
Dex 
iG hae ae SU ee 
ECD RS a ee (2a) 


We evaluate (8) forA=0. We note that L,(y) has a pole for y=2 as a,(2) =0 and 
further poles to the left. If wy Sw, then the pole with the largest real part gives 
the greatest contribution to the integral (8) and we may write approximately 


[Artal Cae = (22) 15 (SEO) (10) 


Ww Ove 


In the same approximation one can obtain the well-known formula 


ip q(Wo,2%, 6) dl = (22) LiG) = 5° Sr eee (11) 
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Thus, defining an average A,? as 


= |. Aj2q(ty, w, £) dt / | a(e09,%0, 2)@l vases (12) 
. “0 
we find 
Se ee es 
Aj? = - 
k ( a2 i Silhts (13) 
In particular we have 
A= ( p ) =():5/0 
85} y=4 


— 2 fed oe 
LS (a tae ) rel 
y=A 


Q>d,—a, a e~a,—a,; 
(Janossy 1948, §§ 556-557). 


Bue 
The integral (8) can be approximated for arbitrary values of A in a similar 


way. ‘The most important singularity in y is attained for y=y, where y, is the 
largest root of 


a DO ea Sa) Ss ed eae (15) 
We have thus 
eae 1 /wy\¥2 D-a OL AY +2) 
eA 2q(Wo, w, 6) dl = — 1) aa Cane Gamer IE 
iF 1K Y(Wo, , C) dl 1, (= a,'(y)(a,— a) OA 


For simplicity we replace = by the equal sign, and write y instead of y,; taking 
the inverse transform of the above equation we have thus 


1 
Ay?4(@o 2,0) = 5—— | 


2m11W, 


tren Wo YD —a,o*Lj(y+2) dx 
Iy—t0 Wy ag—a, on ay'(y) 


We change the integration variable from to y with the help of (15) and find 


Lo rte Wy\4 D — a, 0*L,(y + 2) 
= — —- — d 
Imi, ieee exp { ay (©) @,—a, or Be 


A,?q(o; w, ¢) =; 


The above integral can be evaluated by the saddle-point method. We note 
however, that, apart from the last factor of the integrand, the integral is one 
which is often used to represent g(w ,w,¢). Regarding therefore the last factor 
under the integrand as a slowly varying function, we can replace it by its value 
at the saddle point. In this approximation we can divide both sides of (16) by 


g(Wo, w, 6) and find 
: 2k TL ”, 
Te (aa) tora eee (17) 
A= — (Yo) 


The value of y) has to be determined in the usual way from the saddle-point 
condition, namely 


—a(y)$+In (22) + (in ies 


ay — ay 


) =0 ere eee ve (18) 
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Equation (17) is written explicitly as 
gis woe =a ee 
° [ai +2) —ay(y)] Laaly +2) — (9)] 
= z DONG ATE ee eee (19) 
*  [a(y +2)—a,(y)]}* ay +2) — ay(y + 2) 
De ay +2)—D 
Fe 1 ah VEE Glee eee a 
[a(y+2)—a,(y)Pa(y+2)-a(yt+2) J 
Equations (18) and (19) give a parameter representation of the lateral spread. 
Comparing (19) and (14) we see that (A4;,”),-2=A,?. Thus the spread of the 
shower at the maximum is equal to the spread averaged over depth. 


§ 4 
In the Table we have collected relevant numerical data. ‘The results are 
also represented in Figure 1. The abscissa is log w,/w, ordinate . The straight 
lines represent regions of constant spread; these are labelled with the corre- 
sponding values of A,. Orbits of constant shower size are also plotted; these 


£=a(y) log (w/w) —B(y) log (w/w) 

yy, A,(y) A,(y) oy) By) C=14 C=25 
ont 0-433 0-117 0-092 — 0-009 == = 
ity 0-510 0-195 0-244 —0-005 = = 
113) 0-563 0:268 0-426 0-012 —= = 
1-4 0-604 0-389 0-630 0:042 = == 
ios) 0-638 0-412 0-856 0-084 — = 
1-6 0-668 0-488 1-101 0-140 = = 
127) 0-694 0:568 1-368 0-213 = = 
1:8 0:717 0-654 1-658 0-305 8-628 sss 
ley) 0-737 0-748 1-976 0-418 C207 — 
DO 0-755 0-851 D325 0-554 6-260 = 
ie Oe77Al 0-968 2o7 0 0-718 5-429 9-487 
2) 0-785 1-100 3-141 0-913 4-748 8-250 
2S 0:798 1-250 SAL 1-140 4-180 UPD 
2°4 0-809 1-422 4-163 1-403 3-700 6-342 
Mes 0-819 1-619 4:772 1-704 3-294 5:596 
Dao 0-828 1-846 5-460 2-045 DES) 4-953 
Di, 0-836 2:104 6-236 2-428 = 4-398 
2°8 0-844 MSY) PAN 2°853 = 3-916 
D9 0-850 2°733 8-099 3-321 == 3-497 
3-0 0-856 3-110 9-208 SE853 = =a 
Sieh 0-861 SI°R)S)5) 10-451 4-389 — = 
Sw 0-866 4-010 IAL eeSy 3X8) 4-989 = = 
3308) 0-870 4-539 13-384 5-632 = = 
3-4 0-874 5-126 15-097 6°319 = = 
BE5 0-878 5:774 16-991 7:048 — = 
3-6 0-881 6-487 ID Oya BESO = one 
SO7/ 0-885 7:269 21:367 8°637 == — 
3°8 0-888 8-118 23-873 9-495 = =z 
3°9 0-891 9-052 26-605 10-394 med =o 
4-0 0-895 10-055 IDO IRS hd 337 a = 


* Average. 
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Figure 1. 


The lateral spread of air showers. 
Straight lines give the location of constant spread in the (log w,/w, ¢) plane. 
Lines of constant shower size (OQ=1, 100, 10,000, 1,000,000 are also shown). 


3-0 


Log Q(H,,H, ¢) 
» 


ii 4 (E=25) 
2 } 10 
=— A(é=14)— 
05 
i t 
oo 4 5 6 7 8 oa 0 
Log (w,/w) 


Figure 2. Lateral spread of showers for ¢=25, 14, as function of energy. 
Corresponding values of log Q are also shown. 
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show the average numbers Q of particles with energy above w. Figure 2 shows 
the lateral spread of showers for ¢=25 (sea level) and also for Gra 4 

As can be seen both from inspection of the Table or of Figure 1 the lateral 
spread of a shower increases with ¢ and the increase continues quite considerably 
above the average value; this result is in qualitative agreement with the results 
found by Borsellino. Indeed from equation (19) it can be seen that both A, and 
A, tend to infinity for yoo. This does not mean, however, that a shower 
can develop an arbitrarily big spread, as can be seen in the following way. For 
very large depth ¢ the average number of particles O(w, w, ¢) in a shower becomes 
very much smaller than unity. This means that at such depths the shower 
has usually no particle left, but in exceptional cases one may be left and only 
very rarely more than one. 


The one surviving particle can have arrived there only as a result of a — 


fluctuation of energy loss; thus it must have travelled without serious collision 


for a considerable length. Along this length this one particle may have moved ~ 


at an angle to the original shower axis and thus it may have moved away from 
the axis. Thus the large spread is simply caused by the moving away of such a 
single surviving particle. 


The question arises, what is the largest spread which can be reached before 


the number of particles has been reduced? ‘To answer this question we note 
that the number of particles in a shower can be represented (Janossy 1948) as 


In Ooms 0G) in =? (» — =) + terms depending only weakly on w, 
1 
i 


and ¢= —lIn 0 + terms depending on y. 
eee 


For values of y for which y—a,/a,;’>0, Q can be made large by choosing 
In (z/w) sufficiently large. On the other hand for y—a,/a,’<0, Q is always 
small. From the tables of a, and a,’ we find y—a,/a,'=0 for y=3-27. Thus for 
y>3-27 all showers are small. Therefore 4,=4:3 for y=3-27 is the upper 
limit for genuine spread, i.e. for spread at a stage where the shower may still 
contain many particles. 


The strong increase of shower spread with depth shows that the spread 


develops over large distances. This makes it likely that showers developing in - 


the inhomogeneous atmosphere show a spread rather larger than that calculated 
for a homogeneous medium. We hope to return to this question in the future. 


§5 
_In deriving expressions for Aj? we have inserted (3) into (2). Instead we 
might have integrated (2) in ¢ thus 


| : Aj24( Wy, », 0) dt = | k (=) foe ie g( yy 2, £) dt : g(w',, w, CL? at. 
ee (20) 


Introducing into (20) the following approximate expression derived from (3) and 
(8) (see also (11)), viz. 


oO 5 ke 
i (tg ww", 0) dl = (sam 14(2)) ae Gh) 


A=2 


aetna | : Ai2q(t0o, , 0) al = L (2) (sR L,()) Sees (22) 
A=2 


¢ 
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‘The equation (22) is, however, very inaccurate, the inaccuracy arising from the 
fact that equation (21) is valid only for w, Sw’ whereas in inserting (21) into (20) 
it is used essentially for values w)~w’, as this region gives the most important 
contribution to the integral in (20). This inaccuracy accounts for the discrepancies 
mentioned at the beginning of this article. 
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ABSTRACT. Following a method due to Harish-Chandra, a correspondence between the 
Proca and the $-formalisms of the meson field theory is developed. An important feature 
of the method is that one can at any stage of calculations pass over from one formalism to 
the other. No explicit representation of the B-matrices is required nor need one know the 
exact relationship between the elements of the wave function ¢ and the Proca field quantities 
U,, and G,,. This introduces considerable simplification in the introduction of nuclear 
interaction as compared with earlier treatments where one has to choose a special explicit 
representation of the f,’s and identify the various elements of the matrix 7% with the com- 
ponents of U,, and G,,,._ The nuclear interaction is, in fact, here introduced very simply in 
an elegant way. The quantization is performed for the vector, scalar and pseudoscalar 
cases assuming in all cases both charge and dipole types of couplings. As an illustration of 
the method we calculate the scattering of positively charged mesons by protons. ‘The spur 
calculations, so tedious in the earlier works, are now replaced by simple relations which are 
easily handled, and the results are obtained with about as much labour as is used in Proca’s 
case with the additional advantage that both the longitudinal and the transverse parts of the 
meson field can here be treated together. 


INTRODUCTION 


OLLOWING a method developed by Harish-Chandra (1947, 1948), we estab- 

lish here a correspondence between the wave and the particle aspects of 

the meson theory. One of the important features of this method is that 
one can at any stage pass over from the f-formalism to Proca’s equations or vice 
versa. No explicit representation of the 8 matrices is required, nor need one 
know the exact relation between the components of the wave function % and the 
Proca’s field quantities U,, and G,,.__ This introduces considerable simplification 
in the introduction of nuclear interaction as compared with Wilson’s (1941) 
treatment where one has to choose a special representation of the fs and identify 
the various elements of the matrix y% with the components of U, and G,,. ‘The 
nuclear interaction can, in fact, be now introduced very simply in an elegant way. 
The spur calculations also, so tedious in the former calculations, are now replaced 
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by simple relations which are easily handled. These equations are quantized 
here for vector, scalar and pseudoscalar meson fields. As an illustration of the 
method, we calculate the scattering of a positively charged meson by a proton. 
The calculations in this case have already been done by Wilson (1941) and our 
result agrees completely with his. This method has been applied by the Delhi 
group of physicists working with Professor R. C. Majumdar to a number of 
other problems. 


$i. DHE FREE FIELD 
(i) Some Results of Harish-Chandra’s 
It is well known that the wave equations of the meson field can be written 


in the form 
enOU 
1p! aa + yb O0N a Oe I Se er ee (1) 


where # is acolumn matrix with 5 elements in the case of scalar mesons and 10 in t 
the case of vector mesons; f’s are four 5 x 5 or 10 x 10 matrices according asthe 
wave field is scalar or vector. The 8,’s obey the commutation rules 


sd edcnetled Matyi cS eS omen: (2) 


We follow here the convention that the Greek suffixes can take the values 0, 1, 2, 3 
while the Latin suffix can take only 1, 2, 3. The fundamental metric tensor 
£2» 18 given by 


£00 = —£u = — 822 = —$33= 1, 8 =9 (uA). 


The equivalence of the equation (1) with the wave equations of the Proca 
field can be shown by taking a suitable representation of the 6’’s. Harish-Chandra 
(1947, 1948) has, however, shown that the equivalence can be established without 
making use of any such explicit representation. We first give a summary of his. 
results with slight changes of notation. 

The Hermitian conjugate of (2) is 


pr pret + pepe pe = ginpet + grape 
B,"’s thus follow the same commutation rules as the 8,’s. There thus exists a 
matrix A, which can always be chosen to be Hermitian, such that 


IN BANS Se BH. a a ee een ee 3 
Writing S 
Aes Dam Wesel yes rt) ae el ner cc (4) 
we get the equation conjugate to (1): } 
OF ay 
tae BYE sell oP ae eters (5) 


It can be shown that for spin 0 and 1 there exist one column matrices ‘|, such 
that %"], transforms as a vector. Also, we write 


FG re MS ee eee (6) 
BSA eee ee eee ee (7) 
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The following relations are then easily established 


p*=B, (La) 
BB, + B.B =B,. (1d) 
le (Ic) 
he) (Id) 


The relations (I) are independent of the commutation rules (2). In the case 
of spin | the following relations hold: 


BiB Be =u ly + 80 bs (Va) 
TB Be =Vy'810— Te Bues (V6) 
Pak =P 1,, (Ve) 
ee ate. (V d) 
B+ 38,68" =1. (Ve) 
For spin 0 the following relations hold: 
BLBy | = |.8r0 (Sa) 
pee sale" (Sd) 
BL = 28h |, (Sc) 
LP, oil > (Sd) 
B+28,86"=1, (Se) 


BBB. =8uwP B+ PBS (Sf) 


With the help of ], andT_,’, equations (1) and (5) can be transcribed to Proca’s 
form if we write: 


For spin 1 [ipa a z ee Se tee (8) 
c 
hese es 
[| B,xs Ve Gree eMail eet (9) 
and their complex conjugates. Multiplying (1) by I7,’, we get 
0G 


Es eu), =0; 
on =p a 


and, multiplying by I,'8., 


0U,> ou, 
i Ox! = Ox” . 
: 1 
For spin 0 ry = V/ (hex) 1 O Red eters (10) 
+ 5 Xx 
,' Bx =e hic Ogi eersienene (11) 
ae, 0U 
giving am Ani 


QF 
and ax" +y7U=0. 
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(ii) Lagrangian, Energy-momentum Tensor, and Current Vector Density 


The wave equations (1) and (5) can be derived by the variation of a Lagrangian 


M _ [ Late, 
oa; as a he 
where Bee mie pee a pv | +5 re (12) 


L” vanishes if the wave equations are satisfied. For energy-momentum tensor 
we obtain 


the Os xt 
T= o. | By a Aye Bt | _ Le 


= FEV ee — BoB BxB M+ Se aes fH BAyB.~ BAB .BodYh 


T,,, is not symmetrical in p and v. We define therefore 


Dy = wo on = [yt (B sPaP v — £,B,,.B,)6] 


hex ,. 
= FEY Gn BiBr— BB nates (13) 
We note | 6, d3x =| Ty dx. 


When translated in terms of G,, and U,,, (13) agrees completely with the energy- 
momentum tensor usually employed for Proca’s vector meson field. For scalar 
field, however, one finds that the two expressions differ by a minus sign. This 
situation can be remedied by taking for the Lagrangian for the scalar meson field 
an expression (—1) times equation (12). As, however, in actual problems this 
is not of much importance, we shall retain for the Lagrangian the same expression 
(12) for both scalar and vector meson fields. 
From (13) we get the Hamiltonian density 


CN NON Mis 
H=6)= a b nob = zat Aggie <a eee eee (14) 


This is positive definite if we choose 


N= ip. <a iy ee eee (15) 

where Ho = iS Zhi ee ee (16) 
We notice that 7,*7= land. A 7, 
The current density four-vector 

Gone 
hime ree Bibs aera ele 6 le (17) 
so that, with the help of (15), 
Skew 
ee? Box. SOGO.6-0 (18) 


The expression (17) also agrees with the corresponding expression in the case of 
Proca’s field, except for a minus sign in the scalar case. 
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This establishes complete correspondence between the two aspects of the 
meson field. We proceed next to transform separately the free vector and the 
scalar meson fields to momentum space. The treatment follows closely that of 


Booth and Wilson (1940). 


(11) Vector Meson Field in Momentum Space 
is easily seen to satisfy the second order wave equation 


a) 


ee 2 = 
Ox,0x” gs oe 
We can therefore introduce a plane wave solution 
Up eSD Ch iRpty PI KX) eo. 6 2 ie. (19’) 
Beers 6. LO) Ry (ke ie") 2, 
For a given k there exist six independent solutions of (1). We label these by 
the suffixes /, m,n,.... ‘The most general solution of (1) with a given k is therefore 


a linear combination of all these six solutions: 
6 
4 =1/ (4m) 2 aq exp (—tRoxptikx), ...... (19) 
1=1 
Ro, = Roz = Roz = Ro; | 
Rog = hos =Rog = — ko. J 
We have here introduced 4/(47) for reasons of normalization and have supposed 


the eigenvibrations to be confined in a unit volume. 
We postulate the relations 


where 


10 is 
* 0 
p> U4) usm ym Eis 1 2a ere 0. ond Od on (21) 
10 Ne 
and Pea Weg) ast pas ome | A 5 oes (22) 
A,B 
10 6 . ; 
so that = CR iG, Ay thay, EXP Ror — one or 
A=1 im=1 
6 
y * 
=)CRy ys a, a aol eis) sie (23) 
i=l 
10 10 % . s 
and Jo=e & Da, U4, (Bo) 4p%mUBm EXP (Roi — Rom)Xos 
A=) Lav=1 
6 
* 
=e & q (Bo)im@m sees. (24) 
(=a 
ae one 
Sly Oy Gi LOY Gi Pee ae ey rman (25) 
Tn =i 


where we have taken a representation of the B’’s in which f, is Hermitian and 
diagonal, the first three diagonal elements having the values 1, the next three 5 ile 
and the rest are zero. It then follows that B’’s (k=1,2,3) are anti-Hermitian. 
We may also derive an orthogonality relation inverse to (22): 

10 


6 * 
x pa {(Bo?) 42¥p1}( Po) ttm (Bo")cn} = (Bo) an» A, D=1,2,...10, 


B,O=1 Lm=1 
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which can be written simply as 
6 L 
XZ (Bo?) a(Bolin(m Bo )e=(Bo)ap vee (27) 
l,m=1 


One can easily derive the initial conditions analogous to those employed in 
Proca’s method. Multiplying (1) by 1—£,? from the left and noting that 


Bo? =o, Bo By + BaPor =Bin 


J fe) 
we get b = Byrb— sf 5, Pao" a ce, 2 ee eee (28) 
Similarly, from (5), Wawpt =SE B eh, eas (29) 


(28) and (29) do not involve time and serve as initial conditions. 
Applying (28),"(29) to (14) and noting that 


BoP ro =9, 
Dee — EX] wea 3 So Baebes Se | 
- a wpe {i SEP ba, eee (30) 
where (kB) = —R'6,. 


(iv) Scalar Meson Field 


The most general solution corresponding to (19) in the vector case is here 


by=V (47) X amy exp (— tox) +7kx), Ae. ee (31) 
1=1,2 
where Rote Ree Kea = hoes = ee eee (32) 
The orthogonality relations corresponding to (21) and (22) are 
u * ko 
x Uy Lam — Sims (/, m— 1, 2) MSS er ek (33) 
A=1 x 
and ; ee 4, (Bo)ide Ym = (Bolin: 0 ea ee (34) 
The Hamiltonian density H=hek5 2G; Gd; ae a (35) 
1=1,2 
and charge density Jo=Oat (a, (Be) 0e ee” ee (36) 
1=1,2 
=e(Gp Oi=— ds Gs) eo A = Sree (37) 


We have here, too, taken the representation of Bs in which By is diagonal, the 
first element being 1, the next —1, and the remaining three zero. 
The inverse orthogonality relation is 


p> x {(Bo *) BY B1s(Bo)im Min (Bo Jeo} = (Bo) AD? ASD= V2... 5 


Cb Milly A 
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which can be written simply as 


> (Bor )ia(Boln (tm Bo')s=(Bo)ape vanes (39) 


l,m=1 
The initial condicene are the same as (28) and (29) and the Hamiltonian 
density is given by (30). 


NTE RACTION IW IT He NUCLEAR PAR TICLES AND 
QUANTIZATION 


With the help of the matrices |, and [,,' introduced in the previous section, 
we can now introduce nuclear interaction in a simple way. We treat the vector 
and the scalar cases separately. 


(i) Vector Meson Field 
We assume that in the presence of nuclear Particles: the Lagrangian density is 


given by 
Fe | + ee (A) iT RH 


_ the = ui ee 
+ Seav(Rex)Pw TYPE V] os... (40) 


(g, and g, are coupling constants of dimensions of charge e and of a length times’e 
respectively). 


Ox" Axl 


a = bY uTXP d; Ie = ip'y VT xp . 
Here ¢ is the compound wave function of the neutron—proton 
$= ($8)  =#'0) 
$x 
obeying the absence of mesons in the wave equation 
ee eR: 
ty a + % (Mptp+ Myry)¢ =9, 
(M,=mass of proton; My=mass of neutron). 


The »’s satisfy the commutation relations y“y’+y’y“=2g” and the matrices 
Tp, Ty, TNP» Tpn are given by 


10 00 ND VE goes 
Sea CO eT NO 1 eee Od TPN = O Ou 


We shall also write 
yt yo = 162,13); (eae Vie ot: 


From (40) we obtain the wave equations 


Pew a 4g . x a 
Dee paar i Gen) R, V+ i2mg¢ J) so ae WSL We (41) 
0 4 Bi oN ah ee ee 
and ae BY + ybt= hex i TR“ —i2ng, J(é) Balla ees eae (42) 
Translated to Proca’s form these equations reduce to 
OG” p 2U"=4ng,R : 
Ox" 
CHES MICH Ui 
and G= Oxt = ax = + 47gP yy 
PROC. PHYS. SOC. LXIII, 9—A 69 
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From (41) and (42) we can derive the initial conditions corresponding to (28) 
and (29) for the free field. We obtain 


= porb— "ByBet SE + IB R(t —BayY-+i THe Pll ~ BEY 


* xv (hex) 
ad OR i ee (43) 
ee Ome z Of" 2Rk 4783 *Tet(] — B2 ; ; 2782 *T et = 
p' ='Bo tuesoae Oe SOE a [' (1 Poe) = (hex) x) ee [ pr Bo?)- 
vet Sos (44) 
They will be used to eliminate the non-independent components of ¢. 
From (40) we derive the energy-momentum tensor 
the Ob = ast 
18 a eS = 
T= (ee v ax! Ox, Bt | Lge 
We define again 
t6 -Ow a, 
6 = ti o= Ox. [¥"(6,8,,8,— BBB], 
@ 
so that the Hamiltonian density 
hex RO if A 
H == gan oh — ii /(~) b no PRo- 381 il “\R a Tne 
p . 1 ; eo 
= gS2v (hex) 08" PP os + Z82V (hex) Po TF Bn ou 
4 he Tut 
+321 5 [R, s+ Rb] 
+ Foo (hex) [P,btB’Y — PTB] 
Ney,» WON os ae 
= Fe vrnob tan |) tee PR, + RT 
+ Sav hex )[b BoB" BPs = Pog Te B05 ye eres (45) 


From yo we eliminate the non-independent components of % and ot by 
means of the relations (43) and (44). We then get after some simplification i in 
which we use the relations (I), (V) and the commutation rules (2) 


HH oe 0) eee (46) 
where Fae ee xy ee 5 - BoB. BiBo? x | cf tie (46a) 
and 
HMN 4 (~) E ms Bolo * TR oT Bose — URE RT| 

= aay (hed| — Page T+ 5 Put TGS iy PO — Pog THB gh |. 


see (46 b) 


4 
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We have here replaced ' by "np. It is to be noticed that (46) is Hermitian. 
When translated to Proca’s form (46) agrees completely with the usual Hamiltonian. 
To pass over to the quantum theory, we postulate the commutation rules 


[(Bovx)a» ("(%')Bo?) 5] = — 47(Bo*)and(X-¥’) eee (47) 


and for the Fourier’s expansions of % and %* we take 


3 6 
by =V/(47)z, {3 Ayy U4) + Zi bait Paes aad tar phe (48 a) 
3 6 
ba = a/ (477) X;, 12 Cee + bua} Cae Raeicae rene (48 b) 


a,, and b,, are g-numbers obeying the commutation rules 
[4x Qi 1=1, [b4, 634°] =1 (no summation over &, 1), 


all other combinations commuting. a;," and a, are emission and absorption 
operators for positive mesons, and b,," and b,, for negative mesons. 
Making for ¢ the expansion 


p = UC V ze, 


we now easily get (bar over H denotes that integration over 3-dimensional volume 
has been performed) 


rs. Anh g 1s ape 
BS=—g,, |( =) pS: (KK +k) {> Ay Uy + UY Oxy un} 
NX K, K’,k iL (4 


x Vy Cy {Bok + x0} LC eve’ 


wv 


h - pe he eee Sl 
+0,/(= ‘) L (KK +k) > Ajigy + UX Oxy Mid f 
= l=4 


AG K,K’,k 


x Oy CR {Yl %n&m 7 Cine — XVo%ns se Cyvp 


Ant ere ee ee : 
+a/( | Ss o(K—K —k)vp Cy, Tin {BoR™ + xa" Cex 


x K, K’,k 
3 6 
x {= Ay Upt & baat} 
Pi 1=4 
Anhc\ , Pentel ae 
Se: ee x i(K—-K —k)vp Cx [ LEV 0(%n%n a Chri %a) 


xX K,K’,k 


3 6 
+ XV %njCRUp. { DY ag” yt busta} 


2 
44a 38 8(K —K’' +K"—-K”) | & Oy ae Crop.) 


K, Ky Kz, K’”’ 


x (Oy "Cg Candp'") +8 2(vy Cx (yos)CKep’) 


x (Vy""C Re" (Yo8) Cures”) | soseee (49) 
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(ii) Scalar Meson Field 


For the Lagrangian density in ea of nuclear particles we assume 


pa the the [weg 5 _ Opt re |+ 5 hex ae, 


ax! 
) he ae 

+ po) (S) RT RV 

+ gon/ (hey) LP ee re (50) 
where R=¢'txpd, P,=¢$'y,tNp?- 
From (50) one obtains the wave equations 

oy img J (4 ‘) 

7p" aye Od a WAGE pea An2> PASS ae (51) 
7 fee: _ 118, cod mw Le 2S + Du 

and ee OF" + yp" = Sure rtp tng, |(£) [Pee tonee (52) 


From these equations one can easily pass over to the usual equations of the Proca 
field. 
From (51) and (52) it is easy to derive the initial conditions 


ONO We se 4 Img ly 4785 te” 2\ qu 
bBo = AaB + EL RU BBP SB P= By) 
and Se ee (53) 
Po 4ip2 1 Ob" 9 a0% mS ‘ _ 8 2)\pR* ang or + _ 2 2) pe* 
Pati + © 5 Bop — — ET — BR" — ETA a 

REA a 


From (50) we get the ee for the energy-momentum tensor 


wm Be [ HBSS — Sout | ~ La, 


Defining again Oy =Tyy— 5 5 LH"(B,B,P.— BB, 


we get 


= FEW Cie BaP BBM + 581] (2) (RYE EAT RTA Bl 
— sev (hex BEVP,+ PTB. ees (55) 


Eliminating the non-independent components of % and %' by means of (51) 
and (52), we get 


H=H™ +H, (56) 


#16110) eel e 


M _ he 2 i; Ost wy 
= — FX) wie 5 Fe POP ABPe se | oan. (57a) 


and 


ne he F pre 6 eo 
z =ie.,/ &) [Ry"Bo ly — RT o' Bort] — gav/ (hex) [bt yPy + Po Tote] 


he rs) il : 02 F 
— 185 JAP a Boly —TelBos Pa’ | euepeke tore (57 b) 
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Making now the Fourier expansions 
P= V/ (Am) E (ag py + Olle, 
: Pei Oman Odi, Jeo ty 1 Alene (585) 
we easily get 


TIMN G 4irhe 9 * 
AS = —ig, ei ) =H [SK-K’' +k)oy "Cg yoC ge Dp (Attar +57, Ux2)Bo by 


x KR’ kk 


+8(K—K! —k)Po'B (ag Ui + btie)op Cr VoC Ken’ | 


Anhe a 
+e./(™) x [8(K-K’'+k)(ayuyy” +5, Uy") 
x K, K’,k 
x vy Cx"(x — Bok &») bCrvp’ 
= d(K AS k) Ttup"C;,"(x — BoR™&n) hCx ee’ |] 


tang’ = UK-K’+K" —K")(0y"C "aC 0’) 


,K’,K’,K’” 


Kop Cer OC gy) nnn es (59) 


(iii) Pseudoscalar Meson Field 


Proceeding as in the scalar case it is easy to derive the following expression 
for the interaction Hamiltonian in the case of pseudoscalar mesons : 


aes Aniic\ . ae 
AS sig, a ( - “) E [8(K-K’ + k)oy*Cg*yoxta?a5vp' 


xX K, K’,k 


«(Aug +, Ups )Bo ly 


+8(K—K’ —k)['Bo(ay Ups + Oita) Op Ce “youra?aP Cp dy’) 


4nhe + , * * 


X K, K’,k 
X Uy "Cx" (ixata2a® — Boo,,k™) hCevp’ 
—6(K—K’—k)[fy'up Cg" (txata%a3 — Byo,,k™)C xy’ 
X (Ay Uy + Bj2tt x2) 


Anes a. (KK + KK), Cy oC ze) 
K, K’, 


oye 


x (vp CR OC gn Dy”). in: (60) 


SiGe TAPPLICATION TO THE SCATTERING OF POSITIVE MESONS: BY 
PROTONS 


As an illustration of the method developed above, we shall here calculate 
the cross section for the scattering of positive vector mesons by protons. ‘The 
scattering is described by the scheme 


Y+(po) + P(— po) > Y*(Po) + Y+(p) + N(— Po— Pp) Y*(p) + P(—p), 


where by. Y+(p) we denote a positive meson, by P(p) we denote a proton and by 
N(p) a neutron each with momentum p. We have to make a second order 
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perturbation theory and the interaction matrix elements for the transition are 
easily seen to be 


A, =i 2 (Fs) v'( ? Po) Tin (PD Bo PE pica” )o( iP ae p)u(P) 


hx? 


Amc 


Hne= — 8 a (Fs) uy; (Po)? (—Po— P)(BoPo” + pc”) 1,,0( — p) 
(pp =tk, p=tk, w=hx). ...--: (61) 
Here the suffix i denotes the initial state, f the final, and n intermediate. The 


suffix I’ is used to denote the polarization of the initial meson and / for that of the 


final meson. Hence 
Hy, ne 
te E; aa Es 
4ng °C l 5 T T( -_m m 
= ee ooo oom (= + pe 
hys (E,— EP — WP —2(p, +p) (= Po) in'(P Bot pen’) 


x {Ey — E—c(a, po + p) + Me? }4y(Po)uy (Po)(P"Bo + vex”) 1,0( — P), 


and 
4ar*g ihc? 3 +¥(Am mf 2 
Ai, A = REO?” (= po) Fin'(p™ Bo + ca”) {Ey — E—c(a, pp +p) + Mc?yo} 


X Uy(P)u;"(Po)(BoPo” + vex”) 1,,v(— p) 
—v"(—p)lP."(Bopo” + wca” )ue( po) uy“(p){Ey — E—c(a, py + p)+ Mc?yo} 


x (p™ Bo + ca”) In?(— Po), 
where 


Q=(Ey—E)?— M*ct—c(pot+p), E=c(p?+pyc?)¥?, Ey =c(p?+ M%c?)1? 
Ip|=|Pol=2, (PPo) =p? cos 0. (M=mass of nucleon). 


Summing over the final states of the nucleon and averaging over the initial ones, 
we get (iy denotes initial nucleon states, fy denotes final nucleon states, and by 
iy, and fy we denote the corresponding meson states; spy denotes spur of the 
nucleon matrices, i.e. of yo, a etc.) 


4 
BALM Hy = pap a sal (Po) BoPs" + He2" Fy {By —clap) + Met} 
xT "(Bopo” + pcx” )uy(Po)ey"(p){Ex — E—c(a, po +p) + M Cryo} 
x (p™Bo+ pcx”) 1, {Ex — capo) + Mc? yo} 
Ti (D™ Bo + pcx”) {Ex — E—c(a, py + p)+ Mc?yo}u;(p)]. 


We now write Sr ereleres 


3 
ea ee Uy (Po){BoPo” aie pca”} |, {Boy = c(ap) a Mc*yo} La Bobe” a pce.” }uh( Po) 
Snthes (64) 


and to make this summation over the meson states we use (19’) and the initial 
conditions (28), (29), so that 


ae Bota(Po)t4r"(Po) Bo = = {Ebo + a ( p™B,,)?— perl Rie Sehie (65) 


( 
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We then get 
== PuPat 7pm! {BB+ PEE _ yes) 8, t0,— (ap) + Mey} 
+ Filn'Bo{ BB+ PHBL” — tl ye, pots clap) + Metye}am 
+ Fo" Ex —clap)+ Meryain!{ BB, + PHB _ yet} 8.1, 


+ He aE —c(ap) + Mc?yo}xT,,'Bo {BB 4 ee m Pam pet} Bale. 


To simplify further we note the relations 8 ©. ©) 
fire 3 De Boo leis 
Pn'Bol, =F" BoB Bs hy =Tin'BBsBo 1, = 9, 
lin ly = Pras | eres (67) 


ep eton b =§mn8rs— &msE rn» 
Ve Boke aoe = Ss mr ns: 


Hence 
ral pit(c*p? + spc) Ey + Ec*(pop)}— $urct Mc?y, — EE xc(apo) 


+ 3(uict— cXp?)c(ap)—se(apo)(ap)(apo)]- seve (68) 
Denoting by x’ the matrix obtained from (68) by interchanging p and py we get 
from (63) 


S42 AyHy= 


iv, fu in, fn 


2n2g,4c2 
3h yS Eye? 
x x" [Ey — E—c(a, po +p) + Mc*yo]}, 


and we obtain the differential scattering cross section 


Spy (x[Ey — E—c(a, po + p) + Me*yo] 


ol 


2 A dQ 
ot Ee) pee 
where 
A =8Ey*c*p* + 8EEx3c*p{ — E* cos 6 + c2p?(1 +4 cos 6)} 
+ Ey?{18% — 20E*c?p?(2 + cos 0) + 2E%c4p4(23 + 18 cos 6 + 4 cos 76) 
+ c6p8( — 24 + 8 cos 6 + 16 cos? 6} 
+ EExc*p*{E* — 10 +22 cos 6) + E2c2p?(5 — cos 6 — 8 cos” 8) 
+ c4p*( — 10 — 18 cos @ + 28 cos? 6 + 4 cos? 6) } 
+ E®c*p?( —9 + cos 6) + E4c*p4(35 — 13 cos 6 + 12 cos? @) 
+ Ec&p%( — 43 + 23 cos 6 —4 cos? 6 — 12 cos? 0) 
+ c8p8(17 — 11 cos 8— 8 cos? 6 + 20 cos? @). 
rene: (70) 


This agrees with the cross section obtained by Wilson. ‘The great simplifi- 
cation of the procedure must, however, be emphasized. 
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A New Method for Determining the Radial Distribution 
Function 


By KAI-CHIA CHENG 
Department of Mathematical Physics, The University, Edinburgh 


Communicated by M. Born; MS. received 23rd January 1950 


ABSTRACT. It is shown in this paper that Green’s integral equation resulting from the 
Kirkwood approximation for the radial distribution function can be transformed into a 
differential equation of infinite order, to which the method of successive approximation can 
be applied without ambiguity. The solution at moderately high temperatures is discussed. 
The result is a correction to the usual Boltzmann distribution function; the radial distri- 
bution function oscillates with decreasing amplitudes at small distances and tends to unity at 
large distances. The validity of the approximations is discussed in two appendices. 


§1. INTRODUCTION 


HE present paper proposes a new method of treating Green’s integral equa- 
tion (Green 1947) resulting from the Kirkwood approximation for the radial 
distribution function of molecules. This method consists in transforming 
the integral equation into a differential equation of infinite order, which can be 
solved by successive approximations. However, we only present it as a method, 
and no numerical solution of the equation is attempted. The solution at 
moderately high temperatures is discussed. The result is a correction to the 
ordinary crude approximation exp{—¢/RT} and shows that at relatively low 
temperatures and for high densities in the gas region, the curve of the distribution 
has a minimum, in addition to the usual prominent maximum. This is fully 
confirmed by experiments on the scattering of x-rays at low temperatures, and 
is also in agreement with the results obtained by Kirkwood (1942). 
The validity of the approximation methods is discussed in two appendices. 
Some criteria are given. In an example for nitrogen, use is made of these criteria 
to determine the vapour densities and pressures at several temperatures near 


saturation. _ They are in fair agreement with the observed ones even in the first 
approximation. 


¢ 
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§2. TRANSFORMATION OF GREEN’S INTEGRAL EQUATION 


Green has shown that the distribution function m,(r) in equilibrium states is 
determined by the solution of the following integral equation: 


BENS 2 aan ee t+r ‘(s 
l es) alas 2 42 r i. ’(s) 
at? BRT. al J_¢ £°) —— ital + t) — my? jna(s) RT Ud 


where , is the average density of molecules, r is the relative distance between a 
pair of molecules and (7) is the potential function between them. Equation (2.1) 
can be transformed into a differential equation and can be solved by the method 
of successive approximation. Define first 


g(r) = 2 


; and LUCAS: (8 are (2.2) 


Then (2.1) becomes 


Ing(r)+ = i [_@-2 f(r+2) (9) BO at EO) 


It can be verified by the expansion of f(r) into power series of 7 that the following 
relation holds: 


SOSA) | Cee oe hl ae (2.4) 
where D stands for d/ér. Hence (2.3) becomes, with the help of (2.4), 
Ing(r) ent = Zi |” | "(st #)eP f(r) (5) 4'(s) dt ds. v0... (2.5) 
RTY) 


Treating D as an operator which operates on the variable r and which commutes 
with the variables s and t, one can perform the integration 


ic (s?—2?) ee? dt= = (sDcoshsD—sinhsD). ...... (2.6) 


With the help of (2.6), (2.5) can be reduced to 


47, 2 


“TA | pa(eD coshsD—sinhsD)f(0) $(9)E(9) 48. +++... (2.7) 


Expanding (2.7) in powers of D, one obtains 
8am, = 
OU SES p= USS OUD pear (a) ae Bek We aio 2.8 
Ing(r)e REO D (r), (2.8) 
where M are the moments defined by the expressions 


ure II, g()d'(NPttdr,  v=1,2,3,.0.. «one (2-9) 


In the following section the solution at moderately high temperatures will 
be discussed. In Appendix 1 the validity of such an expansion will be considered. 
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§3. SOLUTION AT MODERATELY HIGH TEMPERATURES 


We shall approach the solution of the equation by the method of iteration. : 
For the sake of convenience in illustration, assumption will be made that only | 
M® and M®) will be needed, while the remaining terms can be ignored. Then 
(2.8) becomes 


Ing(7) eMOME = ST (14 MOD YG). sae (3.1) 
Vf 


The iteration method consists in the evaluation of the distribution function of 
the (+ 1)th order g,, ,,(7) from the left-hand side of (3.1) in terms of the pth order 
quantities M, M,® and f,(r) from the right-hand side of (3.1). For instance, : 
let g,(r) be the solution of (3.1) in the jth order, and define =| 


M,0= =| giryp'nrrar, | 


CO SS ee elt” Rh eg a os Bee 
M2 = = | gte)$' yr dr, | C4) 


Sidr) =718,(7) — 1. 


Then the (w+ 1)th order distribution function £,41(7) 18 given by 


87, 


ln 2)43(7) er? = RT (VL ® 2 MODE) ee teen (323) 
870 7 
or | Zia? N= exp ee CR ar Pa) OL) Ce a a : 
Taking, as a trial function, 
Era e a (3.4) 


and substituting into (3.2) for M,, M,® and f,(r), one obtains to the first order 
of approximation, 


1 
&1(r) =exp kT { —(r)+ =e (M+ mpeDyfse)| sae see (3:5) 
with My = : [ * e POIRT by) 78 dr, 
/0 
M,” = = a eT POMEL (yy dr, ecceve (3.6) 


Sot) = 4 Cmaddland ae 1). 


With the help of (3.5) and (3.6), (3.1) can be solved to any desired degree of accuracy. 

In Appendix 2 the conditions for the validity of the iteration method will be 
discussed. An application of this method will be given in the following example. 
Take the potential function 4(r) as follows (Lennard-Jones 1924): 


b(r) =Ar-? — ur-8, eeas oa) 
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Retaining M® and M® in (3.5), one obtains by a straightforward computation : 


£;(r) =exp — {sl +2n, pesor= 1)-B” (e 2 a zs (f) } | is 


ee egg) rg ee: (3.8) 
where B’ is the usual second virial coefficient, being equal to 
, 2a fr? |, eed | 2m—1 
iB (— eal ie ee =e? m 
Sp) MOS ami?) 
oo ol ee ‘ 
pee = a (CO SAS ee el ee te et ree (3.9) 


_mfA\Ar 77) S21 fi 5 
(aD) {r(5) - BE ar? (Ge)? } | 


and where y =(?/ART)#2. 


B70 


gl) 


Figure 1. g(r) at high temperatures. 


0 r 


Figure 2. g,(r) at low temperatures for large values of m. 


Figures 1 and 2 show the general behaviour of the curve £,(r) given by (3.8). 
At high temperatures, it nearly coincides with that of e kT For large values 
of m, and at low temperatures, at which B’ and B" become negative, the curve 
begins to show another minimum before reaching the maximum. This effect 
becomes more and more pronounced as the density 7, becomes larger, in agree- 
ment with the experimental evidence of Eisenstein and Gingrich (1940) on 
liquid argon at 90°K., which shows a similar curve. Our theory here indicates 
that these phenomena exist quite generally in liquid states. 

However, it is expected that at very low temperatures this approximation 
will not be valid and other methods are required. 
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NI DVAS INGE MOS I 
THE VALIDITY OF THE EXPANSION (2.8) 


In order to discuss the validity of the expansion in (2.8), let us consider the 
integral, 


ary | atl D ps | (sD cosh sD — sinh sD) ¢'(s)g(s) f(r) 4s, 
which, on integration by parts, is reduced to 
T= 5 ssinhsDG)SOas, dao’ (A 1) 
0) 
1 
=—a| e(t)h(t)dt. ———saeaee A2 
where G0)=- gp] OK OA (A2) 


G(s) becomes of the order of ¢(s) when s tends to infinity. Since * 


[- Walar= |" Ofratey}ar 


converges absolutely, one can express f(r) as a Fourier integral, 
fi) = | h(a) ei dr, 


so that the operator D can be regarded as a purely imaginary quantity zx. Let 


, sin as 


i= es GO = Se (A 3) 


sin as 


Then ie ce 


So 


| G(s) | ds = Ofs9°S(50)}; 


which can be made as small as one likes by choosing sy sufficiently large and 
independent of « Therefore (A 3) converges absolutely and uniformly. In 
consequence an upper limit sy can be chosen such that for any small quantity «, 
independent of «, 


"SiN 5p 


—— s°G(s) ds| <e, $555)» aes (A 4) 


I(«)— | 


On the other hand one can expand the function 


0 4S 


°sin as 


N= {= = 52G(s) ds 


in powers of «: 
T(a1, 89) = 5 M(sy)( = tly 2-2 
it 
1 
where ue j 
(= Bal Di}, Cds tees (A 5) 


* r 7 
Here the symbol O is employed to denote ‘ of the order of ’ when r tends to infinity. 
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Replacing ix by D, one has 
Lisi) = 2 Mons.) D2, 
1 


and with the help of (A 4) one obtains 
hea S Ie gy Ole r)et 8 aes (A 6) 
1 


It is clear from (A5) and from the definition of G(s) in (A2), that the ratio 
M°*(s9)D”f(r)/M(s9)D”—?f(r) will increase when s, increases. Consequently 
the number of terms of the series in (A 6) which are required to ensure the accuracy 
one may desire in the sum, will also increase with sy. One observes that when 
Sq tends to infinity the higher moments actually diverge and that the number of 
terms required is infinite. The expression (2.8) is therefore the limiting form 
and has to be interpreted in the above sense. 

For a given value of so, for which the relation (A 4) holds, a practical criterion 
for the determination of the number of terms N required in (A 6) is that 


M¢+0(5_)D™f(7) 


M%s,)D®f(r) <p=<1, epee ate (A 7) 


where p is a measure of accuracy in the approximation. Equation (A 7) also 
defines the relative importance of individual terms in the series of expansion. 

In the example of (3.7) amd (3.8) etc. the condition that M®) can be neglected 
in comparison with M is 


MOD*f(r) 


Mar) <p<l: 


oo 


By the substitution of B’, B” etc. from (3.8) and (3.9) one has 


3¢| L185 —3n®y%(8—1)?| | Fay) 
SUS! 3 gee aap arya eee ese wd hints tt esa A8 
; 2" | exp dy'nX(y>—2)—1 |" | Fay) (A8) 
r T 
“ee REY AN ee ge oe 1 ca el eee rE (A9) 
meee er US Sie ae a (A 10) 
Fy =1-53 —(4-81y + 0-318y? + 0-0667y?+...). cer ees 


¥o =(2A/u)"% is the distance at which ¢(r) assumes a minimum value and 


1 = — $(7)/R. 

If either 7)><r, or T,<T, then 7 will be small and B” can be neglected. At 
moderately high temperatures, namely T7=7), this condition becomes 

2,3 | Lin? —5—3y*8(yP—1)?| _ [5Fi 

exp a(9°— 2)y°n?— 1 Py 


which states that there are regions of 7, for instance, at the position r =7, in which 
B’ is by no means negligible. On account of the high power of 7 in the expression 
(A 11), this condition will be fulfilled at relatively large distances. 


< 


= 


ME ee eee (A 11) 
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APPENDIX 2 
THE VALIDITY OF THE ITERATION METHOD 


The validity of the iteration method is intimately connected with the choice 
of the trial function g(r), which has to satisfy certain boundary conditions. 
For instance, if an exact solution of the equation (3.1) is chosen as trial function, 
then there would be no question of validity at all, provided of course the solution 


is stable against small perturbations. 
First of all, we shall examine the validity of the method at large distances. 


A necessary condition for validity is that lim g,(r) exists and is equal to 2(r); 


ee) 
which satisfies the boundary condition g(r) +1, as r+ 0. From (3.3) one has 
for u—> © 


82, M 820, M® Z 
o(r) =exp— ar Oe ee i ee oO at (B 1) 
At large distances the terms corresponding to the higher moments 


+ MeD*-I(p) =O(G(ry?—*) 


are negligible compared with the first moment 7-1Mf(r) =Oj{¢(r)}. Neglecting 
in the exponent of (B1) terms of higher order in 71 than ¢(r), one obtains the 
asymptotic expression of (B 1), 


e(r) exp {2n, Be) =1)}ae PO, ee (B 2) 
where B is a function of m, and T given by 
2a [° ; 
B=- 3RT | EHEORrRA: Ale eae (B 3) 


which equals approximately the second virial coefficient at high temperatures. 
Now the solution of (B 2) is simply 


o(r) =exp \- (1+ 2n,By"\ te tah cae eee (B4) 


When B>0, namely at high temperatures, the factor (1 + 27,B)~1 of the exponent 
(B4) is always smaller than unity, and it is obvious that if (B4) is adopted as a 
trial function, the iteration method is valid. 

When B<0, namely at lower temperatures, the exponent will become infinitely 
large if m,n, where my is a solution of the equation 


2n, Bin), Ty=— ee ae eee (B5) 


The solution (B4) is singular when 2 =n, in the neighbourhood of which the 
iteration breaks down. The limiting values m) must be in some way connected 
with the densities of an assembly of molecules undergoing a change of phase. 

According to Green (1947), the pressure is equal to 


2 
“a 


—n, kT ~— 
Dee 


| : glr)$'(r)4ar dr =n, RT(L +B)... (B6) 
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By the substitution of (B5) into (B6) the corresponding limiting pressure 
is obtained 


Puasa) saiser neh aes BAe ten as (B7) 


which states that the pressure produced by each molecule, e.g. the molecule of 
a saturated vapour,is about }RT. 

Putting, in the case of nitrogen (Fowler 1936), 4 = 140 x 10-®, A=37 x 10-14, 
so that y? =386-7/T, and assuming that at moderately high temperatures 


2 NN, Sips | 2m—1 
Beeb | 1/2 eS ae m 
3 (;) ‘ iG) a 4 ea abs (B 8) 


approximately, one obtains from (B 5) the limiting density and the limiting pressure 
My = 0-464 x 10?%y4?1 1-22 —0-906y — 0-153 y? — 0-037y3 —... }-1, | 
Po=in kT. J 
It is clear from (3.1) that on letting n, tend to zero, one has 


Ze MORE pes B’, 


so that (B 8) represents the first term in the expansion of B in powers of 2,._ The 
assumption thet B= B’ is equivalent to the assertion that terms of higher powers 
than the zeroth in 7, in the expansion are small. Hence m, in (B 9) is the smallest 
root of the equation (B5) expanded in powers of m,, and represents the limiting 
value of m, in the gas region. ‘Thus m, and py given by (B9) should be near the 
values of 1, and p of a saturated vapour. ‘The Table gives the computed values 
of m) and py at several temperatures for nitrogen and also the observed densities 
and vapour pressure of the same element in saturated vapour states. The 
agreement between the two sets of values is close enough in this rough approxi- 
mation. 


Saturated Vapour Pressures and Densities of Nitrogen 


Temp. No (X10-*4/cem’) my (gm/cm*) — my (gm/cm?*) Po (atm.) Lean 
(° K.) (obs.) (obs.) 
108: 8 2-39 0-110 DRG) 13-60 
116°8 Boil 0-125 — 27-24 20°47 
119-4 2-76 0-128 0-118 = — 
1DAlkeS 2°95 0-136 0-164 a — 
£2322 2:98 0-138 _ 30-83 26°75 
124-9 Sas 0-141 0-200 — — 


Finally, the validity of the approximation method at small distances will be 
considered. In the following, a tentative criterion will be adopted, as a rigorous 
treatment is not available in the present stage. 

Let g,(r) be varied to g,(r)(1+¢), « being a small quantity; g,,,,(7) is thereby 
changed to g,.,(7)+«W(r). According to (3.3) and (3.2) one has 


Y 


SMO =M,%,; SM =eM,®, 


16 9 
so that Wr) =2,41(7) { : oT (M+ M,D*)r(g(r) — y} . 
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‘The ratio 


as «W(r) = 8usil”) Ger (M+ M©D®)r(g(r)-4 \ 
PA IEO IOS batt 
gives a measure of the relative deviations of g,(r)’s of two consecutive orders 
produced by sucha variation. A small r, means that higher order approximations 
are not sensitive to small deviations of the lower order approximations and thus 
ensures great stability and rapid convergence of the sequence g,(7). “Then each 
higher order iteration will bring about closer approximation to the exact solution. 
Hence a condition for validity of the iteration is 
|7.|= 24 
Su 
Equation (B10) also gives an estimation of the closeness of the approximation 
in the pth order. It is interesting to note that at large distances (B 10) reduces to 


{ar (Mo fia M,D*)r(g(r) — p} | Ce 1 aL BS oic 0.4 (B 10) 


Y 


San, uo 
f Bre 


=) .0F 2 bliee 


in agreement with the condition given by (B5). 
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Statistical Equilibrium and Radiation Damping 
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9th May 1950 


ABSTRACT. The reason why the usual form of radiation damping theory does not 
satisfy the principle of detailed balancing is discussed. A form of damping theory satisfying 
statistical equilibrium is developed. 


mesons during the transition of a system having a continuous energy 


spectrum from a state q to a state p lead to a transition probability p,, | U 


where U,,,, satisfies 


[ine damping equations for the scattering, creation or annihilation of 


pd FE 


Uyq=Hyg—in | Heid, Uo 2 a ae (1) 


* Now at Christ’s College, Cambridge. 
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r runs over all final states (i.e. states of the same energy as p and q) and p,dQ, is the 
density of such states in the solid angle dQ... Hf, is the compound matrix element 
connecting r and p. 

The damping equation (1) was derived by Heitler (1941) and Wilson (1941) by 
assuming an exponential time decay of the occupation amplitude a, of the initial 
state gq. The usual perturbation theory treatment assumes a,=1 for all time, and 
gives the well-known expression p, | H,,,|? for the transition probability. 

It was shown by the author (Hamilton 1947, to be referred to as I) that an 
equivalent formulation in terms of approximate eigenfunctions of the coupled 
system of matter and radiation could be given, and that it led to equation (1). In 
this treatment the sum of the occupation probabilities of the various states is 
normalized, so that a decrease in the occupation amplitude of the initial state g 
accompanies any increase in the total occupation amplitudes of the final states r. 
(Furthermore, the problem is regarded as being the calculation of the transition 
probabilities arising from given matrix elements. The latter are to be determined 
by the methods of Schwinger and Tomonaga.) 

U,,,1s not a Hermetian matrix, and several examples have been given (Hamilton 
and Peng 1944) of processes involving elementary particles in which|U,,,|? 4|U,,,|?. 
Consequently it appears that the principle of detailed balancing is no longer valid 
when the damping correction is included. However, a closer examination of the 
damping solutions (according to either of the methods mentioned above) shows 
that they are, in principle, not compatible with detailed balancing. The transition 
probability p,|U,,,|? is derived by considering the system initially to be in the 
state q only, and as time proceeds there will be a finite probability of the system 
being in any one of a set of final states 7, one of whichis p. Similarly, p,| U,,, |? is 
derived by considering the system initially to be in state p only. ‘Thus U,,, and 
U,, do not refer to processes which are in any sense inverse. For detailed 
balancing to be valid in the usual sense it is necessary that probability for 
transitions from gq to p be independent of the occupation probabilities of other final 
states ry. As it is an essential part of damping theory that these probabilities 
are related, detailed balancing becomes meaningless in a damping problem. 

It is also important to notice that the normalization of the occupation probabili- 
ties over the initial and final states means that the process described by U,,, is not 
related to statistical equilibrium. In statistical equilibrium the initial state will be 
anixcure- of states. q, Go, .-..- and the normalization condition has to include 
their occupation probabilities. It will certainly not be true that a solution, in 
which the occupation probability of a single state decreases exponentially from 
unity, while those of the other states increase from zero, can in any way represent 
statistical equilibrium. It is necessary to choose a solution which satisfies initial 
conditions more appropriate to equilibrium. 

In the remainder of the present paper it is shown that the damping approxi- 
mation as developed in I can be adapted to give results which do not violate 
statistical equilibrium. It is necessary to consider an ensemble of systems such 
that their initial states are consistent with statistical equilibrium. ‘The 
approximations which have to be made at various stages are of the same types as the 
approximations which were necessary in I. 

The notation and methods of I will be used. Any state ‘I’ of the system can be 
expanded in terms of the wave functions },. of the free electrons, photons, mesons, 


etc. 
P= Dab, = 2 a,d,exp (—1E,t). 
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, is the time independent wave function, and £, is the corresponding energy 
(i=1 throughout). The amplitudes a, vary with time according to the equation 


a,exp(—1E,t)=2c,ALexp (= 14,1) es eee (2) 
where the A, are the energy values of the interacting system, while the A’ are a 
normal orthogonal set of amplitudes satisfying 
NARS AS a ee (3) 
The H,, are the usual (compound) matrix Sets connecting the final states s and 


r. The c, are determined from the initial values of the amplitudes a,(0) by the 
relation 


c,==A’a,(0), 


where A is the complex conjugate of A. 
The inital conditions for the ensemble of systems in statistical equilibrium can 
be written in the general form 


@,(0).4,(0)=0, ifrAp;  a,(0).a,(0)=|f(Z,) ,.- + (4) 


where f(E,,) is some ‘slowly varying’ function of E, (i.e. f(£,) varies little over the 
line breadths involved in the problem). ‘The double bar denotes the average over 
the ensemble of systems. Equation (4) can be written in the form 


a,(0) =f(Z,) exp (29,), 


where exp 7(p,— ©, )} —0 af hoy 

The roots A, are repeated in each energy interval AE, (cf. 1§6), and the roots. 
inside one interval AF will be denoted by the index «. ‘The solution of (3) is 
obtained by substituting Af =d,V,,/(A,,— E,), and then the V,, satisfy the integral 


equation 


V.=7 cot {(Aq—E,)n/AE} | H,p,dQ,V yy 
the integration going over all final states. Normalizing the V,,,,so that 
| p,dQiV.,.\Vig oe ee ee (5) 
the condition X| Af |? =1 gives |d, |? =(AE/7?) sin?{(A,— E,)m/AE}. 
The variation of the a, with time is given by 


a, = |d,P=U,.f(E.) exp (263) VsaV ral (Ay,— E,\(Ay— E,) exp f(E,—A,,)t} 
7 s : 


Summing the A, over the different energy intervals gives 


ee ~ 2a UrflEs) exp (2¢,) [exp (WE, =e E.)t} os 1]/(E, rx E.) 


“8 (WAB)AE,)E exp (id,) E VscV ppt Un | ee ae ee (7) 
[s a 


) 
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where X(n,) denotes the summation over all E, values except E,, while © denotes. 
Es [s] 


4 the sum over the s values corresponding to the energy value E,. Also, 


ily Ue Ve Ke 


rea Taga. os TS (8) 


U,.. satisfies the integral equation (1). In(7)«=fA£. Using the relation 


Sexi AN = Lith = By 


=0 otherwise 
it is easy to see that 


ad, = | Fe,, E =U, Ural exp {(E,- E,)t} aa 1] 
x [exp (WE, =“) E,)t} 7 1)/(E, eg, eae? = E,.) 


+ (AE) 8(E,, Ey)E(2V pgVegt 18 ,)(2V paV sa t¥ Uy) 
[s] « 


AE) 13 Ei v1 Up EV pV sg + 120 ps) — Uy EV pV sq 18 )F 


In deriving (9) and in what follows, it is assumed that it is sufficiently accurate to 
take the values of U,,, V,,. and {(£,) as those on the energy shell E,=E,. Also 
Bo(L,, E,)=1, if £, aa » and zero otherwise. ‘The summation over EF, splits up 
into the two cases em Hyde =e Ee, giving 

a,a,,=(1/AE){1 —8(E,, E,)} | - 2ir—x)JaU,.U,.— 2ixXU Ups 
[s] 
ca | 
= AK CE cay Ury) | [exp (W(E, —E. eel )ty = 1) | E,) PAGES oa E,) 
fe = =. 
+ (1/AE)*8(E,,E,,) | 2nx— MejeC Umer tee (Us Cle) 
[s] 
* 250, Up | PCOS pas pe ot es ee ae ak (10) 
[s] 
Here the orthogonal property of the V,,, is used to give 


SOV Va) Upp AEDV 5 | 


[s] « 


(ye MON He EB Cpv a og OR 


sa ps pr 


From (8) 


Upp — Uy = — (21/1) ZV pV pq: 8in?{(Ag— E,) 7 AE}, 


while the orthogonality relation (5) gives 
ps es U »,=(AE/n*)= iZeVes sin?{(A,— E,)m/AE}. 
[s] rod 


Thus (10) reduces to 
Gd, =8(Eyy Ey)(1SE)EV rgV pal AE) P field) 
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Taking the s values as a discrete set on the shell £,, the orthogonality relations (5) 
can be written = 
(L/AEB)ZV .V og = Sap» 
Is] 
so that the other set of orthogonality relations become 


(ARS 0 Mg eee UR ] 
; Stes z 


Hence equation (11) shows that the initial equilibrium conditions (4) are invariant 
with respect to time. 

Although the mutual interaction of the final states, which is necessary for 
radiation damping, makes it impossible, in general, to use the concept of detailed 
balancing, it is clear that radiation damping does not lead to any effects which are in 
contradiction with the general ideas of quantum statistical mechanics. 
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LETLEERS LO--VoEe Ee Dito; 
A Photographic Method for determining Half-Lives 


This letter describes a photographic method of determining the half-lives of radioactive 
materials. 

The method is similar to that of Zuber (1948), who compared the intensities of 
blackening produced by placing a radioactive metal on a photographic film at different 
times during its decay, the exposure times being so chosen from a previous knowledge 
of the approximate half-life that the images were of roughly the same degree of blackness. 
The same procedure was adopted for the present experiment; however, owing to the 
stronger source available, the decay could be observed over a large number of half-lives, 
and thus greater accuracy obtained. 

*4Cu was chosen as the isotope whose half-life was to be determined. A copper specimen 
consisting of 1 cm? of spectroscopically pure foil of thickness 1-2 x 10-* cm. was irradiated 
with slow neutrons in the B.E.P.O. at Harwell, and the experiment was started 36 hours 
after the foil’s removal from the pile. Ilfex x-ray film was used and the first exposure was 
8 seconds. During subsequent days, further exposures were made on the same piece of 
film, the exposure times being increased so as to’maintain approximately the same degree 
of blackening. Thus after 39 hours an exposure of 68 seconds was given, after 120 hours 
an exposure of about 100 minutes. At the end of the experiment, the film was developed 
for 10 minutes in ID 19, thorough rocking of the dish being maintained to prevent uneven 
development. 

Another piece of Ilfex film was exposed for 60,.63, 66, 69, 72 and 75 seconds, 39 hours 
after the start of the main experiment. The intensities thus obtained covered the range 
of intensities expected during the main exposure experiment, and provided a convenient 
“ step-wedge ’ for evaluating these exposures. "To compare the various images, positives 
of these two films were obtained by contact printing on ultra contrasty bromide paper. 
This showed up d fferences of intensity well, and intensities could be compared to within 
about 2%. Thus the activity of the source at various times could be calculated. 
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The Figure shows the logarithm of the intensity of the source plotted against time, 
the intensity at T=O being taken as unity. It will be seen that after 120 hours the effect 
of the background due to other activities in the specimen became noticeable, but, on the 
assumption that these have half-lives considerably longer than that of ®*Cu, this background 
can be regarded as constant and a simple calculation then enables us to eliminate it. Thus 
useful observations could be extended over a further two days. 

The half-life, calculated by least squares from the readings plotted in the Figure is 
12-82+0-03 hours. This value, however, must be corrected for the effect of latent image 
fading, which causes the earlier images to appear weaker than if they had been developed 
soon after exposure. The correction was evaluated concurrently with the main experiment 
using a long half-life isotope (250-day ®°Zn), and this leads to an adjustment of the half-life 
of ®*Cu to 12:80-+0-04 hours. This value agrees well with that of van Voorhis (1936), 


12-8+0-1 hours. 


3° 
O corrected tor background 
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The fact that a straight line was obtained in the Figure indicates that in this case the 
blackening depended only on the total amount of radiation received (reciprocity law). 
Since the blackening of the photographic film was mainly due to electrons and positrons, 
this confirms the view that f-rays do not produce reciprocity failure. 

The photographic method of determining half-lives will not in general be more accurate 
than counting methods. It would, however, appear to be useful for determining half-lives 
of the order of seconds, since manipulation at the time of the experiment is extremely 
simple. However, in cases where the blackening is produced entirely or mainly by y-rays, 
it may be necessary to apply corrections for reciprocity failure of the photographic film. 


Research Laboratory on the Physics and E. RaBINOWICZ. 
Chemistry of Rubbing Solids, 
Department of Physical Chemistry, Cambridge. 
23rd June 1950. 
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An Order—Disorder Transition Curve in Cobalt Ammonium 
Sulphate below 0°1° K. 


Measurements have been made of the magnetic properties of a spherical single crystal 
of cobalt ammonium sulphate along its three magnetic axes at temperatures’ below 0:1°K., 
using the technique of magnetic cooling. Along the K, axis * an interesting new effect 
was observed. “ae 

Measurements were made of the dependence of the differential adiabatic susceptibility (as 
found by applying a 40 c/s. alternating field of amplitude about 1 gauss) on steady parallel field. 
Usually such a steady field results in a decrease in differential susceptibility (see, for example, 
Garrett 1950). This is due partly to the magneto-caloric effect and the related difference 
between the adiabatic and isothermal susceptibilities—this effect must always result in a 
decrease in differential susceptibility—and partly to ordinary paramagnetic saturation. “a 
Along the K, axis at temperatures below the Curie point, however, the application of a a 
steady field was found to give an increase in the differential susceptibility. Further 
investigation showed that the susceptibility continued to increase to a fairly sharp maximum 
and then began to decrease again. The value of the steady field required to reach this 


0 100 200 300 400 500 0 700 200 300 A, 
H (gauss) Hf (gauss) 
Figure 1. Figure 2. 


maximum increased from zero at the Curie point to about 310 gauss at the lowest 
temperature reached. The effect is shown in Figure 1. On integration the magnetization 
M is found as a function of the field H along an isentropic. From these magnetization 
curves Figure 2 has been drawn, showing lines of constant magnetization and lines of 
constant entropy on an (H, T) diagram. In constructing this diagram : (i) the entropy S 
was found from statistical thermodynamics, knowing the principal ionic g-values and the 
initial field and temperature from which, on demagnetization, each isentropic is reached; 
(ii) the absolute temperature T in zero field was found as a function of S from relaxation 
heating measurements (method of de Klerk, Steenland and Gorter 1948) made along the 
K; axis (y” was found to be inapplicable at all temperatures when the alternating field was 
applied along the K, axis); and (iii) the magneto-caloric rise or fall in temperatures on 
applying a steady field was found by evaluating the increase in enthalpy 


7H 
IN = MdH and deriving AT= (0AE/OS)q. 
JQ isentropic 
* There are two sets of ions in this crystal. At low temperatures the ions behave as though 


there were only one spin per ion, and are also highly anisotropic. ‘The K, axis is inclined at angle 
of 33° to the strong axis (g=6-2) of the ions in both sets. 
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It will be seen from Figure 2 that the region of anomalous behaviour, where a steady 
field produces an increase in susceptibility and a fall in temperature, is bounded by a 
transition curve of crudely parabolic shape. This curve also is shown on Figure 2, being 
the locus of the minima in the curves of constant magnetization, and must also be the locus 
of points of lowest temperature on the isentropics, since (OT/0H)s= —(O0M/0S) yz. 

Many features of these experiments are in agreement with the predictions of the model 
of Sauer and Temperley (1940), designed in the first place to apply to the paramagnetic 
alums. The authors show that, as a result of dipole-dipole interaction between the ions, 
the state of lowest energy at the absolute zero and in zero field is that for which alternate 
zows of ions are aligned in opposite directions. At temperatures other than the absolute 
rero Sauer and Temperley make use of the Bragg—Williams approximation, and show that 
the appearance of the ordered phase should be confined to a region on the (H, 7) diagram 
bounded by a transition curve very similar in shape to that shown in Figure 2. This curve 
meets the temperature axis at right angles at the Curie temperature 7, and at absolute 
zero meets the field axis at a field He. The values calculated for T> and He; hold only 
for a cubic crystal, so that no quantitative agreement would be expected: but the ratio 
kT¢/uHe should be 2 irrespective of crystal structure. In cobalt ammonium sulphate 
the alignment presumably takes place along the strong magnetic axes of the ions in each 
set, at an angle of 33° to the applied field. Allowing for this, Hj=300 gauss: then, 
taking T~=0-084° kK. and w=3-18, the ratio RT¢/wH¢ is found to be 1-35, in rough agree- 
ment with theory. A feature of Sauer and Temperley’s model is that the transition should 
be second-order between Ty and 37, (discontinuity in (0M/0H)r), whereas below $7, 
there should be a first-order transition with a discontinuity in MM. In practice the transition 
is found to be smooth throughout the accessible temperature range (Tq to 0:-47,). This 
discrepancy, and the departure of the observed kT¢/uHc from the theoretical value 2, 
are probably due to the fact that Sauer and Temperley assume the Lorentz form for the 
local field, that is, the average field acting on the ions is taken to be the same (for a spherical 
specimen) as the externally applied field. In higher approximations this is not so. Both 
the Onsager and Van Vleck theories indicate that the internal field in a spherical specimen 
should be less than the applied field by an amount proportional to (M/A xt)”. ‘It can easily 
be seen that this would prevent any discontinuous change in the magnetization, and would 
result in the appearance of an intermediate state, where the magnetization might be expected 
to adjust itself to such a value that the internal field is precisely the critical field at that 
temperature. Calculation, using either the Onsager or the Van Vleck form of the local 
field, shows that the width of this intermediate state should be about 100 gauss, which 1s 
of the order observed. In this region there should be a marked Barkhausen effect. 

From the point of view of Sauer and Temperley’s model, a more direct test would be 
provided by a crystal containing only one ion in unit cell, such as cobalt fluosilicate, since 
in this case the principal ionic axes must necessarily coincide with the principal magnetic 
axes of the crystal. 

A. fuller account of the experiments, giving the results obtained along all three magnetic 
axes, will appear later. 


Royal Society Mond Laboratory, C. G. B. GaRRETT. 
Cambridge. 
12th June 1950. 
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The Half-Life of the Isomeric State in 077T-m 


The half-life of the short-lived 1°’I'm isomer has been measured using the differential 
type delayed coincidence recorder described by Bunyan et al. (1949). De Benedetti and 
McGowan (1948) found the half-life of this isomer to be 1:0 microsecond; the present 
investigation has yielded a value which is considerably lower, namely 


T,=0-658+0 -003 microsecond. 


The error given was derived from a least squares analysis of the observations shown 
in the Figure based on the assumption that all errors are confined to the counting rate 
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values, the pulse delay times being assumed correct. There may, in addition, be an error 
in the calibration of the delay lines in the recorder. The calibration was carried out both 
before and after the measurement using a crystal-controlled megacycle oscillator. An 
independent calibration had also been made by Mr. D. E. Bunyan of this department 
(cf. footnote to Bunyan et al. 1949). Our measurements agreed within 1% and gave values 
of the pulse delay times within 2°% of those computed from the constants of the delay lines. 


Physics Department, E. W. FuLuer. 
University of Birmingham. 
19th June 1950. 


Bunyan, D. E., et al., 1949, Proc. Phys. Soc. A, 62, 253. 
DE BENEDETTI, S., and McGowan, F. K., 1948, Phys. Rev., 74, 728. 


Scintillations from Naphthalene-Anthracene Crystals 


Anthracene is one of the most suitable organic crystals for the scintillation counting 
of ionizing radiations. It is transparent to its own fluorescent radiation, and gives large 
scintillation pulses in the region of 4000-5000 a., which match the spectral response of 
photo-multipliers with Sb—Cs cathodes. The growth of very large crystals, which are 
desirable for the absorption and measurement of high energy y-rays, is however technically 
difficult with anthracene, due to its relatively high melting-point (216° c.) and the 
deterioration of the liquid when in contact with oxygén. 

: Naphthalene (M.P.—80° c.) can be readily grown in“large crystals several inches in 
diameter, by controlled slow freezing from the liquid state. The scintillation properties 
of naphthalene—transparency, decay time, energy conversion efficiency—are comparable 
with those of anthracene, but its fluorescent emission is in the region of 3100-3700 a. 
This is on the edge of the photo-tube spectral response, and the pulses obtained are only 
about 20%, (dependent on the photo-tube characteristic) of those from anthracene. 
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Large mixed crystals of naphthalene, containing up to 9% by weight of anthracene 
(M.P. 77-80°c.) can be grown in a similar manner to those of pure naphthalene. 
A range of crystals, about 2 cm. diameter x 1 cm. thick, have been prepared, containing 
various concentrations of anthracene, and the scintillations excited in these by Po «-particles 
have been investigated. Measurements were made on several small transparent samples 
cut from each crystal to allow for variations from its mean composition. ‘These gave 
satisfactorily consistent results, and a mean value was obtained for each concentration. 
The a-source was mounted 2 mm. in air from the crystal sample, which was placed in 
direct contact with the face of a 5311 photo-multiplier tube. Integral pulse size 
distributions were plotted in the usual manner, using a cathode follower, linear amplifier, 
discriminator and scaler. A broad flat ‘ plateau’ was obtained in each case, giving the 
total «-particle flux through the specimen, and the average pulse size was determined from 
the discriminator setting at which the counting rate was reduced to 50%. The mean pulse 
sizes for each concentration, expressed relative to those from a pure anthracene crystal, 
are shown in Figure 1. 

The fluorescent emission from the mixed crystals contains two spectral components, 
N (3100-3700 a.) characteristic of naphthalene, and A (4000-5000 a.) characteristic of 
anthracene. ‘The intensities of N and A have been measured independently by introducing 
suitable filters between the crystal and the photo-tube. The average pulse sizes under 


0 


10 


= 
co 


08 
wo 
aN | EN 
& e 
2 06 SHG 
— » 
= s 
S S 
= 0-4 04 
< Ey 
ce co 


02% 


i=) 
N 


0 zi | x | 
01 | 10 100 05) | 10 100 
% ( by weight) of Anthracene % (by weight) of Anthracene 
Figure 1. Figure 2. 


a-excitation were observed as above, and compared with those from pure naphthalene 
and pure anthracene under the same conditions. The total intensities obtained by adding 
the relative proportions of each component measured thus, agreed within 2°% with the values 
obtained from the previous direct observations. This shows the absence of any further 
fluorescent component within the spectral range of the photo-tube, and also indicates the 
consistency of the experimental method. The intensities of N and A, relative to the pure 
crystals, are plotted against composition in Figure 2. 

As the concentration is increased N is progressively suppressed, but all the suppressed 
radiation is emitted as A. This is consistent with the ‘exciton’ theory proposed by Bowen 
et al. (1949) to account for similar phenomena observed under ultra-violet excitation. 
The electronic energy excited by the incident radiation is transferred by resonance from 
molecule to molecule within the crystal, until it is captured by a single molecule, which 
fluoresces. The probability of capture of the ‘exciton’ is greater for an anthracene 
molecule than for a naphthalene molecule because of the lower energy of its excited state. 
On the simplest model, if & is the ‘ exciton ’ capture probability of the anthracene molecule 
relative to that of the naphthalene molecule, and x is the anthracene concentration in 
mols/mol. of naphthalene, then the proportion of N radiation will be 


re 
~ 1+kx 


and the proportion of A radiation will be 1—n. The solid curves in Figure 2, computed 
from (1) for k=43, are in excellent agreement with the observations. 
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This value of k—43 for Po «-excitation is very much smaller than that of 5 x 10* found 
by Bowen et al. (1949) for naphthalene-anthracene mixtures excited by 2900-3200 a. 
radiation at 20° c. The difference is significant and is probably related to the low energy- 
conversion coefficient for a-excitation (>100 ev. per emitted photon) compared with the 
approximately 100% ultra-violet quantum efficiency of fluorescence. 

The ‘ excitons’ produced by the a-particle will be of higher energy (approximately 
ionization energy) than in the ultra-violet case, i.e. the inter-molecular resonance process 
may be one of ionization transfer, rather than of excitation transfer. This will considerably 
reduce the ratio of the differences between the ‘ exciton’ energy and the excited energy 
levels of the anthracene and naphthalene molecules. The dissipation of the incident 
a-energy in alternative non-radiative processes may also contribute to the low efficiency 
and reduction in R. 

The dependence of k on the mode of fluorescence excitation in naphthalene—anthracene 
and other mixed organic crystals is being further investigated, since it represents a new 
fundamental approach to the study of excitation, ionization and fluorescence in solid crystals. 

I am indebted to Mr. F. A. Black, who prepared the crystal specimens. 


Department of Natural Philosophy, J. B. Birks. 
The University, Glasgow. 
10ch June 1950. 


Bowsn, E. J., Mixiewicz, E., and Smiru, F. W., 1949, Proc. Phys. Soc. A, 62, 26. 


REVIEWS OF BOOKS 


Calcul Opérationel, by E. Lasin. Pp. 145. (Paris: Masson et Cie, 1949.) 
780 fr. 


This work may be regarded as a list of general rules of the operational calculus. ‘These 
are given without proof, but the conditions under which they apply are stated rigorously. 
The subject is treated from the viewpoint of Laplace transform and Complex Variable, and 
there is now no doubt—as experience shows—that for any but simple problems operational 
procedure must be based upon both of these subjects. There is brief mention of Abelian 
and 'Tl'auberian theorems, while procedure for the solution of linear integral, integro- 
differential, and differential equations by the operational method is outlined. 

There are three lists of Laplace transforms, two inverse and one direct. The first two 
have the transforms classified as functions of p, i.e. 6(p), the original f(t) being on the right- 
hand side. The total number of entries is 114. In the third table, f(z) is on the left of 4(p), 
and there are 116 entries. It is very desirable that those who use operational calculus to 
any extent should have a dual list, i.e. direct and inverse. Unfortunately publishers in this 
country are averse to issuing lists of formulae of any kind ! 

Although it contains only one worked example, this book is recommended to those who 
have an adequate knowledge of the underlying mathematical theory. The printing is clear, 
and the treatment concise. N. W. M. 


Dynamics of Real Fluids, by E. G. RicHarpson. Pp. vii+144. 1st Edition. 
(London: Edward Arnold and Co., 1950.) 21s. 


Dr. Richardson’s book is unique : it has no Preface and the reader is left to guess the 
purpose for which it was written. ‘The subjects selected for mention range over a very 
wide field; they include a short introduction to the classical theory of a perfect fluid, the 
behaviour of fluids of small viscosity, laminar and turbulent boundary-layer flows, experi- 
mental and statistical aspects of turbulence, some features of compressible-fluid flow, 
fluids with a temperature gradient, hot-wire anemometry, surface waves, break-up of liquid 


jets, liquid drops, the behaviour of particles in suspension, emulsions, some aspects of 
erosion, elastic liquids, thixotropy. 
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Clearly, it is not to be expected that a small book of 144 pages could deal comprehensively 
with these various aspects of fluid dynamics; indeed, the entire space would be needed for 
adequate treatment of any one of them. It is also doubtful whether the book would fulfil 
the needs of those who seek an introduction to fluid dynamics, for some previous knowledge 
on the part of the reader is required. Quite a number of the subjects selected for treatment 
are those on which the author has made his own contribution to knowledge and, undoubtedly, 
the attraction of the book is its esoteric appeal, for it gives glimpses of the trend of recent 
progress in a fascinating subject. A. F. 


Lonization Chambers and Counters : Experimental Techniques, by B. B. Rossi and 
H. H. Straus. Pp. xviii+243. 1st Edition. (New York, Toronto, 
London: McGraw-Hill Book Co. Inc., 1949.) 19s. 6d. 


This volume is the second of the National Nuclear Energy Series (Los Alamos Project, 
Division V) and it has an obviously close association with the first—Electronics, by W. C. 
Elmore and M. L. Sands. The first four chapters deal with the fundamental features of 
ionization and the general properties of detectors based upon the ionization process. Many 
readers will find Chapter 1 of special interest since it collects much useful data on the drift 
velocity and agitation energy of electrons; a good deal of material due to the authors 
themselves appears here for the first time and it refers closely to counting techniques. 
Chapters 2 and 3 are concerned with the general design considerations affecting the perform- 
ance of chambers, and the latter deals in very thorough fashion with the relatively new 
problem of the measurement of variable ionization. In Chapter 4 the basic ideas of the gas 
multiplication process are rather briefly presented. A series of curves gives the values of 
multiplication observed in a variety of gases in typical operating conditions. ‘The studies of 
pulse shape, end-effect of tubes and multiple-wire arrangements follow. While the alloca- 
tion of space to the different sections seems somewhat strange in a textbook, all of it is of 
interest to the specialist in this field. 

The last five chapters consider in some detail the design of detectors for beta-, gamma- 
and x-radiations (Chapter 5), alpha-particles, neutrons and fission products. Chapter 5 is 
comparatively modest in its scope and rather unbalanced. On the other hand, the compact 
discussion of aw-ray spectroscopy is very meritorious. "he two long chapters devoted to the 
neutron form the core of this section of the book; the whole subject is treated with origin- 
ality and the authors are to be congratulated on their thorough treatment of the many 
difficulties encountered in this field. Again, a considerable fraction of the material is 
original. Finally, Chapter 10 discusses fission detectors. 

The book suffers in several directions due to its origin and it is particularly disappointing 
that there is no bibliography. It is a notable contribution to the very limited number 
of texts on a subject of major importance in atomic and nuclear physics. ‘The number of 
engineering drawings is unusually generous. All closely interested in the subject will be 
grateful to the authors for this excellent book. S. C. CURRAN. 


Numerical Analysis of Heat Flow, by G. M. Dusmnperre. Pp. xiv+227. 
1st Edition. (New York and London: McGraw-Hill Book Co. Inc., 1949.) 
38s. 6d. 


In discussing the analytical solution of heat conduction problems by setting up the 
differential equation and solving with specified boundary conditions, the author of this book 
makes rather too sweeping a statement when he says that “‘. .. . this direct approach is of 
little use to the engineer’”’. Many problems of interest to physicists and engineers have 
been solved by the methods described in the well-known work of Carslaw and Jaeger. It is 
true, nevertheless, that occasions arise when it is necessary and desirable to use numerical 
methods. A formal solution may be impracticable, for instance, because of the complicated 
shape of a bounding surface. On the other hand, a solution in the form of a slowly con- 
verging series can well involve computation far more laborious than that required in the use 
of a numerical method ab initio. 
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The book starts with a discussion of the fundamental concepts of heat transfer in a some- 
what non-fundamental manner. Heat capacity, for instance, is defined by referring to 
readings on a watt-hour meter. The British system of units is used throughout. Some 
elementary steady-state problems are solved by simple algebra, and there is a brief account 
of analytical solutions in simple cases. Then, the method of relaxation, extensively 
developed by Southwell and his fellow workers, is used for the numerical solution of the heat 
conduction equation in the steady state, in one and two dimensions. Relaxation, however, is 
not a suitable method for transient heat flow problems, the equation being of the parabolic 
type where the boundary conditions are not specified over a closed boundary. ‘There are, 
of course, other methods available. The one developed in this book is a useful numerical 
adaptation and extension of the Schmidt graphical method. The book contains some pro- 
blems involving convection and radiation. There are numerous worked examples to show 
the details of the numerical technique, and also problems for the reader at the end of each 
chapter. These problems are of a practical nature, dealing mainly with the heating of 
rooms, heat transfer through walls and so on. 

Anyone whose work entails the consideration of problems of this kind would benefit 
by reading the book and working through the examples. He would probably be stimulated 
to learn the full power of the relaxation method by further reading elsewhere. A physicist 
would be surprised by some curious remarks on the theory of relativity and the quantum 
theory, which are misleading and not relevant to the main subject-matter. M. R. HOPKINS. 


Artificial Radioactivity by P. B. Moon. Pp. vit+102. Cambridge Monographs 
on Physics. First Edition. (Cambridge: University Press, 1949). 12s. 6d. 


It is well known that the scientist in the general sense has nowadays nearly become 
extinct, as there is nobody who masters botany and physics, zoology, chemistry and 
astronomy equally well. Unfortunately this process of sub-division has not stopped; from 
the physicist comes the nuclear physicist, and even he has undergone fission into the beta- 
ray spectroscopist, neutron physicist, cosmic-ray expert and accelerator engineer. But 
the subconscious urge of the old days to attempt to attain a universal grasp of nature and the 
burning desire to be up-to-date in neighbouring fields still exists in all of us specialists, and 
we are therefore grateful if we are helped in this endeavour by the appearance in the series 
of Cambridge Monographs on Physics of so neat a condensation as the book on Artificial 
Radioactivity by Professor Moon. 

Its first 30 pages cover the experimental methods used to measure the energy and intensity 
of the beta particles and gamma rays emitted by radioactive bodies. The description is very 
clear but necessarily (and rightly) short, since one could write a whole tome about beta 
spectrograph design alone. ‘The next 40 pages cover processes in which the atomic number 
Z changes. Here electron and positron energy distributions are discussed and compared 
with Fermi’s calculations. One cannot help admiring Fermi’s intuition in constructing at so 
early a date (1934) so lasting a theory of so elusive a phenomenon as the continuous beta 
spectrum. <A very satisfactory section on K-capture and on the experimental evidence, or 
lack of it, for the existence of the neutrino concludes this part of the book. It is followed by 
a section of 30 pages about processes which leave the nuclear charge unchanged. Under 
this paragraph fall such topics as life-times of excited states, gamma-ray selection rules, 
isomeric transitions, internal conversion and angular correlation between cascade gamma- 
rays; the experimental facts are generally illustrated by one or more concrete cases. This is 
followed by a descriptive exposition of the relevant theory unencumbered by actual calcu- 
lations, but the final formula is as a rule given in a form just usable by the experimentalist. 
ee ae cannot be covered ; references to the most important papers will 

¢ penetrate further into the subject. 

The book does not contain a complete list of the radioactive bodies with their emission 
data and modes of preparation. Such information, together with some additions on the 
practical side, might have served those readers who will be interested mainly in using beta 
emitters as tools. The title of the book is therefore slightly misleading as the one suon 
deals rather with the physics of radioactive decay than with the radioactive material. It is 
also not quite obvious why the book is restricted to nuclei of light and medium weight 
only; to include the heavy ones would not have increased its size very much. One 
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serious omission (which may be due to an oversight on the part of the reviewer) seems to 
exist : the discovery of the effect of the atomic electrons on the decay constant of Be(7) by 
Segre and Wiegand (Phys. Rev., 1947, 71, 274), and by R. Bouchez, P. and R. Daudel, and 
R. Murart, 7. Phys. Radium, 1947, 8, 336). That the rate of K-capture of Be(7) depends on 
the value of the electronic wave function at the nucleus may seem trivial, but the net result 
is certainly most astonishing and noteworthy, since it is the first time in history that the decay 
rate of a nucleus has been influenced by external conditions. 

One requirement for a book such as the one under review is speed of publication. As 
the footnote on page 80 indicates, part of the manuscript was written in Spring 1947, the 
book was presumably completed in 1948, and it carries the date of publication, 1949. As no 
reference is made to such important recent additions to our collection of nuclear physics tools 
as the crystal counter (publications started in July 1947, but the method really only became 
generally used during 1948), it seems that one has to reckon with a delay of one year at least 
between the completion of the manuscript and the moment of publication. This clearly 
demonstrates the difficulty of writing anything up-to-date in these days when scientists have 
to bear such a heavy burden of teaching, and when the time required for book production 
has been so greatly increased. E. BRETSCHER. 


Cosmic-Ray Physics, by D. J. X. Monrcomery. Pp. viii+370. First Edition. 
(Princeton: University Press; London: Geoffrey Cumberlege, 1949.) 
40s. 


Although the preface to this book states that it is not supposed to be a laboratory 
handbook, it nevertheless treats the subject with more detail than is easily digested by the 
newcomer and contains much information that 1s of use mainly to the active experimenter in 
cosmic-ray work. For instance, a valuable appendix is devoted to curves connecting 
momentum energy and range for electrons, mesons and protons in air, aluminium and lead 
which are based on more observations than any previously published, and the twenty or so 
pages on the cloud chamber deal with the instrument in considerably more detail than is 
usual. On the other hand, the photographic plate method of investigation to which major 
recent advances are due, receives only four pages, dealing mainly with the capabilities of the 
method, rather than details of technique. 

The book—again despite a disclaimer in the preface—is well provided with references, 
which are rather unfortunately sandwiched in the text, making reading difficult. Lengthy 
mathematics is usually avoided, but the Stoermer theory of the motion of a charged particle 
in a dipole field is included in an appendix. Another appendix is devoted to tables of the 
results from the statistical theory of showers. 

Each topic is brought well up to date, and the features which are still problematical, are 
made clear. A particularly valuable chapter is included on heavy particles tn cosmic radiation, 
a subject which is receiving much investigation at the present time. T. E. CRANSHAW. 


Radiative Transfer, by S. CHANDRASEKHAR. Pp. xiv +393. Istedition. (Oxford: 
Clarendon Press, 1950.) 35s. 


In this book Professor S. Chandrasekhar, F.R.S., who holds the chair of theoretical 
astrophysics in the University of Chicago, makes an outstanding addition to the systematic 
treatises on astrophysical problems. "Though the book is intended primarily for astronomers 
it also has great meteorological interest, because the methods of calculation of radiative 
transfer which it describes are applicable to the very difficult problem of multiple scattering 
of light in the earth’s atmosphere; moreover, essentially similar problems occur in con- 
nection with the diffusion of neutrons, hence the book should also attract the attention of 
theoretical physicists. The main topic, however, is the radiative equilibrium of stellar 
atmospheres, treated with great breadth and also in much mathematical detail. 

The subject goes back at least to Rayleigh’s discussion (1871) on the blue light of the 
sky and its polarization; as far as the earth’s atmosphere is concerned, further progress 
was made by L. V. King, and also by S. Chapman and A. Hammad, but they recoiled 
from the difficulties of a complete formulation of the problem and its solution; 
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Chandrasekhar, in recent years, in a long series of papers in the Astrophysical Journal, 
has boldly faced and overcome these difficulties with great mathematical skill; one element 
in his success is the use and extension of a principle of invariance introduced SS the Russian 
astrophysicist V. A. Ambarzumian (1943). Chandrasekhar, in one of a series of Biblio- 
graphical Notes which conclude each chapter (and form a valuable and distinctive feature 
of the book), points out that some of the basic ideas involved resemble those introduced 
by Stokes and Rayleigh in the discussion of the reflection and transmission of light by thin 
plates. Among other authors whose names occur in these Notes are Schuster, Schwarz- 
schild, Sampson, Eddington, Milne, Hopf, van de Hulst and Crum. 

The book has a good subject index and index of definitions, but lacks an authors’ index. 

Radiative transfer is a subject of considerable intrinsic mathematical complexity, but 
this brilliant and powerful book will smooth the path of knowledge for many readers and 
thereby give a new impulse to further development of the subject. SHCe 


Corrigendum 
“The Infra-red Spectrum of Magnesium Oxide”, by J. C. WiLtmort (Proc. 
Phys. Soc. A, 1950, 63, 380). 


The author has asked that the following alteration be noted : 


The vertical scale of Figure 5 on page 398 should run from 10~* to 10 instead of from 
1Om ston: 
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Voltage Gradients in Long Gaseous Spark Channels, by J. B. HtcHam and 
J. M. Meek. 


ABSTRACT. Oscillographic measurements have been made of the variaticn with time of 
the voltage drops in long spark channels conducting aperiodic impulse currents rising to a 
peak in about 1/4 sec. and decaying to half-value in from 10 to 85 usec., with peak currents 
ranging from 60 to 700 amp. ‘The lengths of the spark channels have been measured 
photographically, and the mean voltage gradients along the channels have been estimated. 
Results are given for sparks in air, nitrogen, oxygen and hydrogen at atmospheric and 
reduced pressures. Some preliminary measurements are described for sparks in water and 
oil, and for sparks in glass tubes. The voltage gradients in the various gases, with the 
exception of hydrogen, are independent of the peak current, but are influenced by the rate of 
current decay. ‘The results are correlated with the current densities in spark channels, as. 
recorded in a separate investigation. 


The Expansion of Gaseous Spark Channels, by J. B. HicgHaM and J. M. MEEK. 


ABSTRACT. A rotating mirror camera with a temporal resolution on the film of 
1:5 mm/sec. has been used to determine the rate of expansion of long gaseous spark 
channels during their first 10 usec. of growth. The spark channels conducted aperiodic 
impulse currents in the range of 60 to 500 amp. peak, which attained their peak values in 
about 1/4 usec. and decayed to half-value in 10 or 28 usec. Sparks have been investigated 
in air at atmospheric and reduced pressures, and in nitrogen, oxygen and hydrogen at 
atmospheric pressure. ‘The measured areas of cross section of these spark channels bear 
a linear relation to the peak currents, except in hydrogen, and are little influenced by the rate 
of current decay in the range studied. Measurements have also been made with the rotating 
mirror scanning an image of the spark across a photomultiplier connected to a high speed 
oscillograph; by this method a direct record is obtained of the radial light distribution 
across the spark channels. The mechanism of spark channel expansion is discussed, and 
estimates are made of the average ion densities in spark channels. 


The Electrical Resistance of Liquid Gallium in the Neighbourhood of its Melting 
Point, by C. Dopp. 


ABSTRACT. In order to determine whether any discontinuous change in electrical 
properties takes place when a liquid is supercooled, measurements have been made on the 
resistance of liquid gallium. No such discontinuity has been found, the resistance varying 
linearly with temperature over the range investigated. 

The results yield values for the resistivity and temperature coefficient of resistance of 


liquid gallium. 


Hydrogen-filled Geiger Counters, by B. CoLLINGE. 


ABSTRACT. Permanent gas counters are discussed and a new quenching circuit is 
described. Quenching times of about 45 usec. have been used and plateaux with slopes 
better than 0-05°%% per volt obtained. A formula for the correction of counting rates for 
losses cz ased by the counter dead-time has been confirmed within wide limits. It has been 
found possible to measure counting rates as high as 5 <x 10° counts per minute with an 


accuracy closer than 1%. 
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Cavitation Produced by Ultrasonics, by B. E. NottincK and E. A. NEPPIRas. 


ABSTRACT. The problem of cavitation produced by ultrasonic vibration is examined 
theoretically. Equations are developed which describe the motion of a gas-filled cavitation 
bubble in a liquid medium subjected to alternating pressure; the case of an empty cavity is 
also considered. Information is obtained concerning the distribution of fluid pressures and 
velocities in the medium near the bubble surface during the motion. It is suggested that 
these theoretical conclusions may be used to show how the intensity of the various effects of 
cavitation will depend on ultrasonic frequency and intensity. In particular, it is predicted 
that all cavitation phenomena will diminish and finally disappear as the frequency is 
raised. The important part played by the nuclei from which the cavitation bubbles grow is 
emphasized. 


The Colour Temperature of Light Sources, by H. G. W. Harpinc. 


ABSTRACT. The definition of colour temperature is explained and the approximations 
of the spectral distributions of energy of lamps to those of black-body radiators are indicated. 
The value of colour-temperature measurements to the colorimetrist is outlined and 
difficulties in making the measurements are dealt with. The use of photoelectric devices 
for colour-temperature measurements are mentioned. 


Certain Properties of Electrostatic Fields Encountered in Electron Lenses, by 
P. A. Linpsay. 


ABSTRACT. The equipotential lines in a bipotential electron lens exhibit a fine structure 
which is revealed by application of the relaxation method to the solution of the Laplace 
equation in cylindrical coordinates. "The asymmetry of the potential field between the two 
edges of the cylinders is confirmed by the same method and is established to be of the order 
of 2°% in one particular case. 


Reduction of the Spherical Aberration of Magnetic Electrun Lenses, by 
U. .-<GiaNnova: 


ABSTRACT. 'The dependence of the resolving power of the asymmetrical bell-shaped 
magnetic field on the lens parameters is examined for varying degrees of field asymmetry. It 
is snown that a high asymmetry index favours high resolution, and a preferable value of 
lens power is specified. 

By suitably reinforcing the lens field with that of a small air-cored coil, the theoretical 
resolving power of the lens is increased. On account of the small size of the coil, the coil 
current must be pulsed intermittently if an appreciable reinforcing field is to be produced. 

The basic factors governing the design and operation of a suitable coil are discussed. 


The Computation of Wave- Front Aberrations of Oblique Pencils in a Symmetrical 
Optical System, by W. WEINSTEIN. 


ABSTRACT. Aberration formulae given previously, together with new ones for astig- 
matism and secondary spherical aberration, are put into forms which are suitable for 
computation and simplify the treatment of field lenses. A method of comparing the 
formulae with ray tracing is given, and the results for a triplet objective are discussed. 
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